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LESTER DURAND GARDNER—Founder, Institute of 
the Aeronautical Sciences; Fellow, Institute of the Aero- 
nautical Sciences, and Honorary Fellow, The Royal Aero- 
nautical Society of Great Britain; Daniel Guggenheim 
Medalist; Publisher; Soldier; Public Servant; Delegate 
at Large in the promotion of aviation throughout the 


world 


Ln Memortam 


We of the Institute owe an incalculable debt to Lester Gardner 


His was 


the vision to see the need for a strong and active association of aeronautical 
scientists and engineers at a time when air power and air transportation were 
selling at a discount. He realized the potential value of such a society to the 
a concept that needed no argument a short decade later 


national defense 
Never one to “‘let George 


For L.D.G., vision was but the prelude to action 
do it,’ he plunged into any project that intrigued him with a fervor that left 
his colleagues and associates gasping. 

And so it was when he took the Institute under his wing 
ent of its founding in 1932 to the day he went on inactive duty in 1946, he 


From the mo 


poured out his boundless energy, day and night, in season and out, to make« 
the Institute the foremost organization of its kind and, through it, to keep 


\merica in the forefront of aviation development 


And so” to Lester Gardner Hail and Farewell! 


Now that he has gone, we can do no less than rededicate ourselves to the 


objectives to which he devoted all the years of his life 

Aviation is vastly more complicated, vastly more important today than 
was even imagined in the Institute's early years. All the more reason, there 
fore, that we should make the best possible use of this vital force that Lester 
(Gardner nurtured, brought to fruition, and left in our hands 

The Members, Officers, and Staff of the Institute can pay their best tribute 
to the Founder by keeping the organization strong and alive to every oppor 
tunity to advance the arts and science of aeronautics. Through its channels, 
we must do what we can to make and to keep this country strong 1n the air 


for peace or for war 








Lester Durand Gardner 


Lester Durand Gardner, who founded the Institute o} 
the Aeronautical Sciences in 1932, died November 23 in 
New York after an tllness of several weeks. Burial wa 
at Arlington National Cemetery on November 27. 


Lester Gardner was an indefatigable worker for 


aviation. One of his last efforts in its behalf was a 
campaign last year to raise funds to establish the Jerome 
C. Hunsaker Chair at Massachusetts Institute of Tech 
nology. The drive raised $500,000. In 1954, M.I.1 
had established the Lester Durand Gardner Scholarship 
Fund of $50,000. 

It was only last summer that the National Aero 
nautic Association conferred upon him the title ‘Elder 
Statesman of Aviation’ —-the last of many honors be 
stowed on ‘Major Gardner,” as he was known to 
thousands here and abroad. He was the third Ameri 
can to be made an Honorary Fellow of the 86-year-old 
Royal Aeronautical Society of Great Britain. 

In 1944, Polytechnic Institute of Brooklyn gave him 
its honorary degree of Doctor of Laws. He was the 
recipient of the Daniel Guggenheim Medal for Aero 
nautics in 1947. He also was a Fellow of the Institute 
of the Aeronautical Sciences, honorary member of the 
Wings Club and The Aviation Country Club, and mem 
ber of the University Club, New York, the Army and 
Navy Club in Washington, the Society of Automotive 
Engineers, American Society of Mechanical Engineers, 
American Legion, the Society of Colonial Wars, and 
Sons of the American Revolution. 

He served in the Army Air Service in World War I, 
leaving it after the Armistice with the rank of Major 

He was educated in the Detroit public schools and 
was graduated from M.I.T. in 1898 with a Bachelor of 


Science degree. Following a year of administrative 
law at Columbia University, he became interested in 
publishing and was successively on the staffs of Th: 


New York Times, The New York [erald Tribune, The 


Sun, and Everybody's and Colliers magazines. 


1876-1956 


In 1915, he became president of the Gardner- Moffat 
Company, which later became the Gardner Publishing 
Company, publishers of Aviation magazine and T] 
Rubber Age. 
who continued its publication for many years befor 
Aviation Week. 


lisher of Aviation that Major Gardner helped aviatio: 


The former was sold to McGraw-Hill 


changing its name to It was as pub 
through many of its growing pains. 

Always interested in scientific developments, Maja 
(Gardner entered the field of radioactivity and use « 
radium, uranium, and thorium. He took out severa 
patents. Following a visit to the Curies in Paris, hi 
imported the first radium to be used for therapeutics 1 
this country. 

In 1915 he attended the first reserve officers training 
camp at Plattsburg, after which he and Philip J. Roose 
velt published American Defense, a magazine that di 
much to arouse both official and public opinion in su 


While in activ 


service in 1917-1918, he organized S9 aero squadrons 


port of United States defense needs. 


for overseas service. 

Major Gardner served on numerous aviation con 
mittees and official delegations appointed to further 
the development of aviation. He was president of the 
Aeronautical Chamber of Commerce of America in 1925 
He also was a former president of the Technology Clu! 
of New York. 

In 1932, with the help of other aviation leaders, h 
founded the Institute of the Aeronautical Sciences, witl 
Rockefeller Center Th 


Institute soon became one of the greatest scientif 


original headquarters at 


organizations in the world, and its activities wer 
directed by Major Gardner until his retirement 
1946. 

Major Gardner also founded the Aeronautical At 
chives of the IAS, which contains one of the worlds 
principal aeronautical collections. He was the firs 
Archives president. 
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Recent Developments in the Theory of 


Wing-Bodyv Wave Drag 


HARVARD LOMAX* ann MAX. A. HEASLET* 
Ames Aeronautical Laboratory, NACA 





SUMMARY 


Phe theoretical study of wave drag on airplanes flying at super 
nie speeds is re-examined in the light of the recent discovery of 
transonic area rule \ supersonic area rule is presented to 
ther with a discussion of its range of applicability In regions 
vhere the supersonic area rule is seriously in error other methods 


yased on linearized theory, for calculating the wing-body inter 


ference effects are outlined In addition, some recent results 


sed on the nonlinear small-perturbation equation for transoni 
w are briefly discussed. Comparisons between theory and 


periment are given 


SYMBOLS 


minor and major axes, respectively, of elliptic 
wing plan form (See Fig. 11 
x { functions giving strengths of wth = orde 
1,2 \ multipole distributions [See Eqs. (380) and 
' 31) 
line coefficients in Fourier expansion of S’,(x, 4 
See Eqs. (8) and (9 
1 coefficients in Fourier expansion of S’(x, 6) — 
B(x, @)/2q| [See Eq. (38 
Cy, ¢ lift and drag coefficients, respectively 
D drag 
6) + streamwise distance between oblique planes 
6 bounding airplane for given 6 
Le lift on the wing and body, respectively 
x, 0 projection in @ direction of resultant force on 


oblique section 
lift on body integrated from nose to station x 
free-stream Mach Number 


mY) difference in pressure on Opposed sides of wing 
free-stream dynamic pressure, ('/2)pol 


Mx, 6 airplane area in oblique plane given by Eg 
3 
fuselage area in plane normal to free-stream 


direction 


6 wing area in oblique plane 
maximum wing thickness 
free-stream velocity 
Vien J volume of wing and body, respective ly 
= 1 ° 2 ° . . 
Presented at the Aerodynamics-I Session, Iwenty-Fourth 


innual Meeting, IAS, New York, January 25-26, 1956 
\eronautical Research Scientist 
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v, VY, 2 = Cartesian coordinates, parallel to free 


stream direction 


‘vlindrical coordinates, x parallel to free 
stream direction 
8 V MW? - 1 


ratio ol specific heats 


u(x, 9 slope of wing upper surface 

A(x, ¥ difference in slope on opposed sides of wing 
p free-stream density 

A(X, ¥, perturbation velocity potential 


} 
semlapex ingle of circular cone 


INTRODUCTION 


‘ X JHEN WHITCOMB demonstrated that a reasonably 
slender airplane and a body of revolution with 
the same longitudinal distribution of cross-sectional 
area (normal to the free stream) could have similar 
wave drags at transonic speeds, a new impetus was 
given to the study of interfering wings and bodies 
Such studies proceeded, of course, along both experi 
mental and theoretical lines. The purpose of this 
paper is to present some of the theoretical methods 
used and some results of their applications 
In particular, we shall limit ourselves to analyses re 
stricted by the assumption that the disturbance ve 
locities induced by the airplane are small relative to 
lL, the flight velocity. 
body, the x axis parallel to the free-stream, and ¢ is the 


If the axes are fixed in the 


perturbation velocity potential, the inviscid fluid theory 
for supersonic flow can then be based on the linear 


differential equation 
Ber — Sw — ¢ 0 | 


where 6 MJ? — 1, J being the free-stream Mach 


Number. Our discussion of the recent results based 
on this equation will, for the most part, be limited to 
analyses heretofore classified. 

As \/ approaches 1, Eq. (1) is, except for very special 
cases, an inadequate idealization of the flow field and 
a general analysis must be based on the nonlinear 


equation 
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o, —[MPe(y + 1)/Uolergr, (2a) 


B° rr =— Cay ~~ <€ 
or, for the particular case of sonic flight speed, 


Puy + vee = [(¥ + 1)/Volerer: (2b) 
ANALYSIS BASED ON THE WAVE EQUATION 


Preliminary Remarks 

To begin with, a certain familiarity with solutions to 
problems involving wings and bodies in flow fields 
governed by the wave equation is assumed. For ex- 
ample, we recall that sources, doublets, and higher 
order singularities, when properly distributed, can simu- 
late the potential flow about such wings and _ bodies. 
Equally pertinent to our subsequent discussion is the 
division commonly made of the total drag into the 
separate parts commonly referred to as the viscous, 
vortex, and wave drags. In practical applications, of 
course, all three types of drags must be considered 
simultaneously. Nevertheless, in the succeeding dis- 
cussion the viscous drag will be ignored, so that the 
expression “optimum combination” will refer to that 
combination having, under given restraints, the lowest 
value of vortex drag plus wave drag. Further, by 
imagining the airplane to be enclosed by an infinitely 
long circular cylinder parallel to the free stream, we 
adopt the usual definition of vortex and wave drags 
as the forces in the free-stream direction given by the 
the 
air- 


flux of momentum across the ends and sides of 
cylinder, respectively, as the distance from the 
plane to the surface of the cylinder tends to infinity. 
The condition for minimum vortex drag is the same for 
both subsonic and supersonic flow and, as shown by 
Munk,! depends upon the horizontal and vertical dis- 
tribution of vorticity streaming behind the airplane. 
The condition for a lower limit to the wave drag will be 
presented later. 

Finally, in the following, considerable use is made 
of a theorem presented by Hayes? in 1946. Let r be 
the radius of the cylinder, mentioned above, across 
which the momentum flux is to be measured. Then 
Hayes’ theorem states: 

(a) To the lowest order in 1/r, as r tends to infinity, 
and at a fixed azimuth angle 6, the magnitude of the 
perturbation velocity potential ¢ in the vicinity of the 
cylinder’s surface is invariant to a finite translation of 
sources (or any other singular solution to the wave 


equation) along the plane 


xy x + By, cos 6 + Bz, sin 0 € 


where x is any reference distance along the axis of the 


cylinder. 


The Supersonic Area Rule 

The very simplicity of the concept suggested by the 
words “‘‘area rule’ suggests a universality that is hard 
to resist. Hence, shortly after the transonic area rule 
was discovered, a supersonic area rule was presented 
by both Jones* and Whitcomb. The development of 
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this rule as given by Jones, and later as discussed }y 
Ward,‘ stems from a classical result (see reference ; 
p. 60) which, for our purposes, can be paraphrased 
follows: 

(b) In linearized theory the potential flow about g 
airplane can be simulated by a distribution of simp] 
sources on the airplane’s surface provided that none 
the thin wing-like members carry a local loading. 

(c) Within the group of airplanes given by theore 
(b) there is a smaller family whose local source strengtt 
is proportional (within the accuracy of 
to the local airplane surface slope in th 


linearize 
theory ) 
streamwise direction. * 

To construct the supersonic area rule for this limited 
family of airplanes is simple enough. First consider 
a line 6 = constant along our enclosing control surface 
Along this line the momentum flux caused by an air 
plane of the prescribed group is, by theorem (a), the 
same as that which would occur along the line if all 
the singularities simulating the airplane were moved to 





the axis of the control surface along the oblique planes 
given by Eq. (3). By theorem (b), these singularities 
are simple sources only. Using theorem (c), one can 
show the strength of the equivalent source on the cyl- 
inder’s axis is proportional to the gradient of the normal 
(to the free stream) projection of the area intersected 
by the oblique planes. But Karman? has derived the 
equation for the wave drag of a line of sources. 
his result and averaging the drag of the equivalent 
source lines for all 6, we get for drag D, divided by free- 


Using 


stream dynamic pressure g, the expression 


21.10 2110 
dé | dx) | dx» x 
J —Lo(0 J —Lo(a 


S" (x1, 0)S" (x2, 6) Inixy — x2) (4 


> 


D | [ 
qd ar / 0 


where S(x, @) is the normal projection of the oblique 
section area (see Fig. 1), S”(x, @) (07/0x7).S(x, 6), 
and the interval —L)(@) < x < L,(@) is the airplane 
length relative to the oblique planes determined by < 
fixed Mach Number and 6. Eq. (4) applies to all air 
planes which satisfy conditions (b) and (c) and which 
further satisfy the condition that both S(x, @) and 
S’(x, 0) be continuous. 

Eq. (4) provides an abundance of practical informa 
tion for those cases to which it doesapply. To illustrate 
this, let us consider a thin wing centrally mounted on a 
body of revolution at zero angle of attack, and let the 
x axis coincide with the centerline of the body. Since, 
according to condition (c), the wing must be thin and 
the body slender, in finding S(x, 6) the variation with 
z, can be neglected when the oblique plane |Eq. (3 
cuts the wing, and the variation with both z, and } 
can be neglected when the oblique plane cuts the body 
as illustrated in Fig. 2. This permits the term S(x, # 
to be separated into two parts, 

*It is the latter family that Ward and Jones consider. B) 
ordering the source strength in terms of thickness ratio, Ward 
argues that the physical meaning of condition (c) is to limit the 
subset to slender bodies and their nonlifting wings or their com- 


bination 
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THEORY OF WING-BODY WAVE 


Oblique plane | 
x, 2 x+fy,cos 8 \ 
+ 82,8in8 ~, \ 


= r Ni 





x, 


7-S(x,8) , normal 
projection of 
obliquely cut 








area 
Fic. | \irplane section in oblique plane 
S(x, 0) = Sfx) + S,(x, 0 (5 


where S,(x) is the cross-sectional fuselage area normal 
to the free stream and, therefore, independent of 6, 
and S,,(x, 6) is the normal projection of the wing area 
cut by the plane x; = x + By, cos 6. In terms of 
\, (x, y), the slope of the wing’s upper surface, S’,,(x, 6 
is given by the equation 

. 


2 Ay(X a. BV; cos 6, v1 dy, (6 
J wing 


where the limits of integration are determined by the 
intersection of the oblique plane with the exposed por- 
tions of the wing. 

The next step in simplifying our results is to expand 
S’,(x, 6) in a Fourier series (by symmetry only cos 


n§ terms are needed, where 7 is even); thus, if 


l 
S'ix. # S-( x wk 4 (x .% 
U, 
! . = 
“eo A, (x) cos né (7 
Uy 2 
then 
l ] Lens - 
Ao(x) S’,.(x, 0)d0 (g 
l , 2a J 0 
l 1 (27 
A(x : S’ (xv, 6) cos né dé (9 
l 7T/0 


ind by Eq. (4) the drag becomes 


dD l s/ al 
ep dx; dx2[UoS" (41) + 
rl ‘ I . u 
A’o(x1)] [UpS" x2) + A’o(xe)] In |x, — xe] + 
} 2 D, (10 
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where 


D, l vl e] 
_ dx, dx» x 
q trl)? I J —I 


Y) — Xe (11 


From Eqs. (S) through (11) the following conclusions 


can be drawn. The fuselage area S,(x) has no effect on 
the A, and D,, terms, so the drag of all the higher har 
monies is a function of the wing design only. The 
complete effect of wing-body interference, therefore, 
appears in the sum of A (x) and the fuselage area dis 
tribution determined by planes normal to the free 
stream. 

Next let us consider some simple examples of how 
low wave drag fuselage designs can be calculated on 
the basis of the above results when the wing and Mach 
Number are fixed. The simplest case is to find the 
optimum modification to a cylinder on which a given 
wing is to be mounted, as shown in Fig. 3. This 
modification is given by 


l fia 
S’.(x, 0)d0 (12 


S’ Ax) — 
2m. 


One can show that the volume removed from the cyl 








y, 





/ 
/ 


Lys 
xzxtycosO 





{ x; 


Simplified evaluation of S(x, 6) for planar wing anc 
slender fuselage 


Fic. 2 








oe. SF 
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Wing enclosing Mach 
cones with 
apexes on 
body axis 
VA % 
x DA 
es 4 
Fr 3. Wing and optimum modification to straight-sided 


evlinder 


inder by such a modification is identical to the volume 
of the exposed wing and that the drag of the combina 
tion is equal to the drag of the (exposed) wing alone 
minus the drag of the modified cylinder alone. 

If, instead of an infinite cylinder, a closed body is 
combined with a given wing, the optimum body shape 
is, in general, more difficult to find, although the prob- 
lem often reduces to a simple reinterpretation of the 
result for an optimum slender body of revolution 
having the same geometric restraints. For example, 
consider the following problem: 


Given the (nonlifting) wing, body length, and total 
volume (sufficiently large) of the combination, what 
is the area distribution of the body which yields a 
minimum wave drag if the apexes of the Mach fore 
and after cones enclosing the wing lie within the body? 


The solution in this case is obtained by mounting 
the given wing panels on a Sears-Haack body (by 
definition, the closed slender body which has a mini- 
mum drag for a fixed volume and length) and then 
modifying this body by taking out of it the wing vol- 
ume in accordance with Eq. (12). This follows from 
the above; first, since the wing is given, the term 
l 
q 
cation given by Eq. (12) the initial term in the right 


> D, in Eq. (10) is fixed, and with the modifi- 


member of Eq. (10) is simply the drag of the body 
alone; second, since the wing plus the modification 
have no net volume, the total volume remains un- 
changed. Notice that the location of the wing along 
the body is immaterial provided the required indenta- 
tion can be assimilated by the fuselage. 
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Returning to our theoretical considerations, we ea 


combine Eq. (6) with Eqs. (S) and (9) and, by mea; 


of a simple transformation of coordinates, write E 


I if A,(E, n)dédn 
TJJo V Bn? — (x — & 


\ 


12) in the form 


S’ Ax) 


and the value of A,,(x) in the form, for 7 > 0 


ot i 
A,,(x I 
Ax(&, 7) cos 1) cos ix — et bn } adésdn 
V Bn? — (x — & 


where oa is the region defined in Fig. 4. 


Eq. (15 It tells 


under th 


has an interesting interpretation 
us that the optimum fuselage area change 

prescribed assumptions) is determined by integrating 
sources over the portion of the wing in the “region | 
In other words, S’,(x) is determined by 
that 
flow field at the point x, 0 in either forward or revers 


silence.” 
planar source distribution cannot disturb th 
flow. 

Another point to be made concerns the problem \ 
finding the optimum wing thickness for a given pla 


form. This is resolved by determining the thickness 


z. D,, in Eq LO 


that is, by minimizing (under the pertinent restraints 


that gives the lowest value of 


the higher harmonics A,(v) independent of the valu 


of Ao(x) which is absorbed into the fuselage analysis 
In the above, some simple illustrations were give 
for finding minimum wave-drag combinations under 
certain geometric conditions and at a fixed free-strean 
Mach Number. 
of Mach Numbers at once. 


It is also possible to consider a rangé 
Define the wing area cuts 


made at a given Mach Number, by 


= l i 
; S i<. @ dé 


) a 
—T7 


Then, using Eqs. (8), (10), and (15), and simple van 


ational procedures, one can show that 





ao = shaded oarea 


Wing panels 





i 





Traces of 
Mach cones 


through (x,0) 





vg 


Fic. 4. Definition of area ¢ 
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rHEORY OF WIN 
— “1 
> As —-->)) Sx, M 16 
N i=) 
‘ves the minimum value to the sum of the drag co 
efficients at .1/ M,, Me, Ms,..., M,. If a certain 


weighting is assigned to each of the drags, as would, 
for example, be the case if the airplane were required 


a given Mach Number in a minimum 


, accele rate to 
time, the best fuselage area distribution can still be 


readily determined. In fact, if, for a fixed wing, 


C,(\DdJ is to be minimized, the solution 


ul 
‘ 


for the distortion of a cylinder is simply 


— , 
§ AM) S’(x, MW)d AM l7 
VV, — WM, Ju 
is an average of the obliquely cut 


Notice that S(x, M 


wing areas, and S’,(v), as determined by Eq. (16) or 


is obtained by further averaging with respect 


tothe Mach Number 


Regions in Which Supersonic Area Rule Fails 


The power of the supersonic area rule is evident from 
the simple and direct geometric interpretations that 
ould be given in each of the above cases. Untortu 
nately, however, there are many combinations of wings 
ind bodies to which the rule does not apply. Equally 
unfortunate fact that the 
validity or the errors incurred by its use in specific 
We are particu 


is the boundaries of its 
cases are not always easy to assess. 
larly concerned here with those situations for which 
the perturbation velocities in the flow field are small 
enough for linearized theory to apply but for which the 
requirements given in (b) and (c) above are not satisfied. 

Combinations of wings and bodies which cannot 
be simulated by simple sources alone, in accordance 
with (b), Any thin wing cam- 
bered with respect to the local stream so as to support 


are easy to construct. 


local loading is a member of such a class, and in ex- 
treme cases the wave drag can be zero although the 
oblique area distributions can be any of an infinite set. 
Combinations of this type were pointed out by Buse 
mann and Ferrari many years ago. For example, the 
wave drag, as given by linearized theory, of a particular 
lorm of the Busemann biplane is shown by the solid 
line in Fig. 5. Notice that the total drag is a mini- 
mum (zero) whenever $//c is a simple fraction of the 
torm 1/2, 1/4, ..., 1/2m,... and is a local maximum 
whenever 8h /¢ equals . iva. | 
the local 


straight line connecting (1, 1) 


Oy «a» Stee ee DT), 
Further, 
Com- 


pare these results with those given by the supersonic 


with the origin. 
trea rule and one sees the latter sometimes overesti- 
mates and sometimes underestimates the correct value 
Ol the drag 

Examples of objects that can be simulated by simple 
sources only but which are still of such a shape that 


t} ° 
their wave drags are not adequately predicted by the 


supersonic area rule because they violate condition 


¢) can also be given. Consider a rectangular or tri- 


BODY WAVE 


maximums all lie on the 
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angular wing mounted centrally on a flat sided fuselage 
for which 8 times the height / of the side is greater 
than the root chord, both wing and fuselage being at 
zero angle of attack. Assuming the fuselage area does 
not change in size or shape in the free-stream direction 
we find that the drag of the combination as given by 
the supersonic area rule is the same as the wave drag 
of the exposed half wing alone. But the correct linear 
ized value is half that of a complete wing formed by the 
exposed half panel and its reflection in the plane of the 
The 
shown in Fig. 6 and the magnitude of the difference is 


wing-body juncture. ratio of the cwo drags is 
clearly appreciable. 

Another example along the same lines is provided by 
some results published by Nielsen’ who studied the 
wing-body combination shown in Fig. 7. For the par 


ticular case shown a _ sonic-edged sweptback wing 


mounted on a circular cylinder the wave drag given 
by the supersonic area rule is about 1.55 times that given 
by linear theory. 

Finally an interesting inversion to the relation 
shown in Fig. 5 between the drags given by linearized 
theory and the supersonic area rule for the Busemann 
biplane is provided by a tube having its surface corru 
if the radius 


gated in the form of a sine wave. Thus 


is given by 


R ro + 


6 sin (27x; IS 


the obliquely cut area (assuming control-surface theory 
and neglecting terms of the order 6° compared to 6) can 
be shown to be independent of @ by reason of radial 
symmetry and can be written 


S(x, 0 2arr6J)(2rBro/l) sin (22x/1 (19 


kind. 
According to Eq. (4), the wave drag of one ‘“‘corruga- 


where ./,(z) 1s the Bessel function of the first 


tion length’’ (determined by finding the average of a 
large number of them) is 


D gq S 2°( 179767 / 17) Jo7(2 Bro /1 (20 


and the comparison between this and the correct linear 
ized value, as given by von Karman,’ is shown in Fig. 
S. Notice the 
those shown in Fig. 5. 


difference between these results and 
In the one case, the approxima- 
tion given by the supersonic area rule gives zero values 
for the wave drag when the correct values are not zero, 
and, in the other case, the reverse is true 


Drag in Terms of Obliquely Cut Areas and Forces 


The preceding examples illustrate the fact that many 
types of wings and bodies and their combinations lie 
outside the class restricted by (b) and (c) above. It is 
interesting to examine these types still using Hayes 
theorem but applying it to a field containing singu 
larities that can simulate any object for which linear 
ized theory is applicable. A simple way of carrying 
out this objective is to combine Hayes’ theorem with 
Green's theorem which relates surface distributions of 
sources and doublets to surface slopes and loadings 


Thus, one finds’ for the wave drag 
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D 2? 21310 2/6 
dé | dx, | dxo X 
r” JO J —Lo(0 J —Le 
r 3/'(x1, 0 i 31’(X2, 6 
S W - 37 (x 4g) — rd 
2q qq 
In x) — x 21 

where /(v, @) is the component of force pointing in the 
§ direction on the oblique section defining S(x, 6) (see 
Fig. 9). 


The similarity between Eq. (4)—-the equation giving 
the drag according to the supersonic area rule——and 
Eq. (21) is immediately apparent. In fact, the manner 
in which the term 8//2q couples with the term S’(x, 6) 
gives the error made in using the supersonic area rule. 
Clearly this error is small when the Mach Number is 
near | or when the force coefficient //g is small relative 
to the area gradient S’(x, 6). 

An immediate application of Eq. (21) is its use in the 
calculation of zero wave drag, nonlifting, shrouded 
bodies of revolution. For example, consider a body of 
revolution surrounded by a thin circularly cylindrical 
tube of negligible thickness, though this is not essen- 
tialsuch as is shown in Fig. 10. Specify any radially 
symmetric loading on the tube’s surface. This fixes 
a value of 6//2q which is independent of 6 due to sym- 
metry. Construct a body having an area variation in 
normal cuts (for simplicity, we are considering only 
those cases for which the difference in cross-sectional 
area between an oblique and normal cut is small) 
identical to this variation of B//2g, and place this body 
inside a straight sided cylinder. The loading induced 
on the cylinder will be just that value of 8//2gq orig- 
inally assumed and, by Eg. (21) the combination 


will have a negligible wave drag. The existence of 


such a body was first mentioned by Ferrari.'’ Gra- 
ham, Beane, and Licher'! have published a construc- 
tion similar to the one given here. 

Another application of Eq. (21), first pointed out by 
R. T. Jones, is its use in determining a lower bound to 
the wave drag for a given total lift. Let 2A(@) be the 
distance in the x direction between the two oblique 
planes bounding the airplane—i.e., Lyo(@) + L,(@) in 
our previous notation. Then, in terms of the Fourier 


expansion, 


Bl(x, 6 


5'(x; 0) = S* C,(0) sin ny 22) 
a a 
2q 1 
where cos y = 7/A(@ 93) 
he drag can be written 
D L a tis 
> n | C,,7(0) dé 24 
q os 3 J0 


Notice that C,(@) can be interpreted in terms of S), the 
base area (for simplicity in concept consider the body 
to be followed by a semi-infinite cylinder having the 


Same shape as the base), and L, the total lift. Thus, 


Ci(6 [2/mA(0)] [So — (BL sin 6/2q 25 
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and the C,,’s for the higher 7’s depend on the moments 
of area and lift. 

Now fix the total lift and assume the area and loading 
> C,7(4 0 


) 


is distributed in such a way that 
The drag of such an arrangement is given by 
D | Sas BL sin 6 dé 
‘ ss 26) 
qd 2xr- J0 2q A-(@ 
If 1/A°(@) is expanded in a Fourier series, so that 


| 4 
S~ (8, cos n0 + a, sin nO 24 
, 


A‘ 0 0 
one can show that the minimum wave drag occurs when 
S, BL 2q (a,/B 2S) 


and is given by 


(D/q), 1/2x) (8L/2q)? X 

By — (B./2 — (a,°/28 29 
Notice that a combination with vertical symmetry 
(ay Q) or with a value of A(@) that is the same for 


all 6 has the lowest wave drag for the given conditions 


if the base area is zero. 


Method of Multipoles 


One general solution to the wave equation (see, 
for instance, reference 5, p. 527) consists of an infinite 
set of basic singular solutions referred to as multipoles 
These basic singularities can be distributed along a 
line and, by appropriate weighting, can be made to 
simulate any potential field exterior to a certain cyl- 
inder along whose axis they lie. The expression for 
the perturbation velocity potential for a distribution 
of nth order multipoles, starting at — Ly) and continuing 
along the x axis, can be written (for a discussion related 
to problems of the type concerning us here see reference 


12) 
,” boy Ff 
g(x, r, 6) - 
2x \r Or J —I 
[a,(x1) cos nO + b,(x;) sin né 
ax; 30a) 
1 (¢ — v) - — 827? 
or, alternatively, 
| . 
Q(X, 7, 6 _ 
2m J =I 
[A ,,(x,) cosn6é + B,(x;) sin n@|cosh }ncosh~! {(x—.x,)/8r}! 
Vix — x — Br 
xX dx (30b) 
A, (x (—8)" (0”"/Ox" Ja, (x (—)"a v)) 
B, (x (—B8)"b,'\"" (x f 31) 


The drag caused by these multipoles is simply 


D = 2D,+ > (D,. + D; 32) 


— 


where D,, and Dg, are given by Eg (11), the A,,’(x) 


being replaced by A,,’(x) and B,,'(x), respectively. 
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Clearly, th 


existence of a stream surface representing a re 


Vertically Symmetric Combinations In this part let where 2L, is the total multipole length. 


us consider a wing-body combination symmetrical il shan 


about a horizontal plane (z Q) and at zero angle of will depend on its distance from the multipole line ag 
attack. 


of the body they are to represent, we see the drag of the 


Placing the line of multipoles along the axis the induced velocity field begins to take effect as wel 


as on the strength of that field, but a more specifi 
wing and multipole line is given by Eq. (32) wherein statement would require detailed information regarding 


Dp, 18 zero (by symmetry) and D,, is now given by a particular case. 


Eq. (11) except that the terms A,,’/(x) are replaced by An example in which slender-body theory ~i« 
A,/(x) + A,'(x). This leads at once to the possibility np ANd Gng ~ 1/7 was used to calculate the bod 
of zero wave drag it shape at _1/ 1.4 is given in reference 12. An ellipt; 
| | 22 wing having an aspect ratio equal to 5 and a bicony 
a A(x * | S(x, Odd (55a airfoil section with maximum thickness ratio of 5 per 


cent was combined with a body having a maximun 


diameter about 1/5 of the wing span. For this cas 
that the full strength of the a) ten 


could be used, but only a portion of the ad. term. A) 


| l f* 
=e? - S’,(x, 0) cos nédé@ (33b . . 
U r Jo it was found 


which defines the optimum strengths of the multipoles é : ; 
experiment details of which are described in the refer 


so placed. 


; ‘ was carried out to test two combinations: one 
In order to illustrate the amount of drag that one 


ence 


the wing described mounted on a basic body (comb; 


can associate with the various multipoles throughout 
a Mach Number and aspect ratio range, the nonlifting 
elliptic wing with a biconvex section shown in Fig. 


nation of a Sears-Haack body and a Karman ogive 
and the other, essentially the same wing mounted on 4 
= wD ; 3 body modified according to the stream surface induced 
1] was analyzed. The wave drag of this wing was first ; 


= ; age by the wing source sheet, the basic body source line 
calculated by R. T. Jones'* and is given by the equation 


ay)(x) and the best possible variation of d2(xv), under the 
The results 


Cp (4/8) (t/2a)? } [1 + 2(a/bB)?] condition that the fuselage remain real. 


[1 + (a/bB)?]/*; 9 (34) for the experimental drag rise are shown in Fig. 13 

a ia : : : Less wave drag was attained throughout the Mach 
Phe reduction in wave drag (details of the calculations : coat 

; : ~ . ; : Number range shown by modifying the body eve 

appear in reference 12) as the wing is combined with a : oe 

though the modified body had more volume. Notice 


optimum multipole distri- 
Notice that for values 


the first three (non-zero) 
butions is also shown in Fig. 11. 


of the Mach Number near unity most of the drag is 


that the agreement between experiment and theory is 
especially good as to the actual drag reduction at the 
: : design Mach Number. 

destroyed by a line of sources only, but as the Mach a ; ‘ ; : 
: Btescae ; : If, in the region affected by the multipoles, the body 
Number increases, the higher order multipoles play 5 
: ‘ ; does not distort greatly from a known surface, the dis 
an increasingly important role. : : 


at iki ca : ; tortions can be calculated by using a control-surface 
rhe principal difficulty in applying this method to 


: é ‘ i : method in which the induced velocities are evaluated 
physical problems is to find the shape of the fuselage 
We shall 


consider here only those cases for which the wing source 


: : it ; on a given surface and assumed to remain unchanged 
corresponding to the multipole distributions. ; an é‘ 

in the flow field on either side. The advantage ol 
control-surface for the 


shapes lies principally in the fact that the distortion 


method calculation of body 


sheet and the line of multipoles simulate the entire 


flow field exterior to the airplane surface. For such 


. ous : ‘ 
sai ; fields are linear so that body modifications induced by 
cases, the possibility of zero wave drag, mentioned ; ‘ % 
coer, ae cs the various orders of multipoles can simply be super- 
above, is a fiction if the body is to have finite length. ‘ ities : ; 
imposed. It is disadvantageous in that accuracy falls 


This follows, since any sources added to those given 5 : ae i 
‘ ian ie ‘ off rapidly as the modified body surface approaches 
by Eqs. (33) must increase the drag; yet sources must ‘ ; 

; om ; : i the body axis. 
be added if the body is to close and be real since the a aes ; : 
ee : : : : : Nielsen’ using a circularly cylindrical control surface 
initial magnitudes of Ao(x.) are negative for a real wing. hae: 

has presented calculations of a body shape for the 
Nielsen found that the 


axially symmetric body distortion reduced the wave 


Generalizations as to what the numbers and strengths : : i 
combination shown in Fig. 7. 


of the multipoles can be in the presence of some basic 
source distribution or starting streamtube to insure the ae ; : 
: : aps drag 22 per cent below that for the wings mounted ot 
simulation of real shapes are difficult to develop. As 


: i : rane ¥ : the undistorted cylinder, and the distortion caused b\ 
an illustration, consider the variation with 7 of the radial : 


the next higher harmonic caused a further reduction 0! 


velocity ¢, (the induced velocity principally affecting ; q 
; : p : : : 12 per cent. At the wing-body juncture (the region 0 
the body shape) induced by a multipole line. For as 


highest indentation) the maximum body indentatio! 


small values of r, ¢, varies as 1/r”*', whereas for large ; ; nr 
: é : 4 required for the full use of the higher harmonic was 
values of y (staying always behind the foremost Mach : er od ; ‘ 
i arts : igs ; found to be six times that required for the lower. 
cone), it varies as 1/Vr. For instance, the magni In applying the method of multipoles as described 





Sew sz 








tude of the first crest of the wave in ¢, (shown in Fig. 
12) is produced by the multipole distribution 
Lo)? |* (4/Le 


a(x) coll — (x 


above, one should always remember that only for 4 
limited number of bodies can the entire exterior flow 
field be simulated by a set of multipoles along a single 
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line and, in general, the optimum fuselage to go with a 
given wing cannot be expected to lie within this limited 
set. Hence, the optimum multipole distribution asso- 
ciated with a nonlifting wing gives, when used directly 
to calculate the body shape, a fuselage better than, 
say, a straight-sided cylinder, but it probably does not 
vield the solution to the physical problem of finding the 
fuselage which combines with a given wing to produce 
the lowest drag for the combination. 

Asymmetry and Lift-—-Next let us consider a wing- 
body combination for which the wing is thin and no- 
where inclined at a large angle with respect to the free 
stream but otherwise of arbitrary thickness camber, 
and disposition in x), ¥1, 2 space. For convenience, 
place the x; axis along the center line of the body and 
let AX(X1, V1, 21) be the difference between the stream- 
wise slope on the opposed sides of the wing and Ap(x, 
yi, 21) be the difference between the pressures across 
the same surface. Now cut the wing with the oblique 
plane, Eq. (5), and obtain the projections in a y), 2; 
plane of the oblique section so formed. One can then 
show that 

, 


L(x, 6) F l 
sin @ Ap(x + By cos 6 + 
qd J wi) q 


. 
; | 
2) sin 4, My 21 \dyy + cos @ 4 
e vu q 
Ap(x + By, cos 6 + Bz, sin 6, v1, 2)dz 5) 
and 
. 
S’(x, 8) | AX(x + By cos 6 + 
82, sin 8, y1, 21) ds; (36) 


where ds; is an element of arc in the y,, 2 plane. Fi- 


nally, let us generalize the expansion given by Eq. | 





to include the asymmetry and lift. Thus, 
Bl(x, 0) 
S’(x, 6) — 
2g 
pi A,(x) cos n6 + B,(x) sin n@] (37) 
0 
where 
a ell Bl(x, 8) 
Ap(x) = > (2, 0) = dl © 
U5 2% J0 2q J ( \9O) 
etc 


If multipoles are placed along the x axis and used to 
simulate the body shape, it follows, just as before, that, 
from the point of view of lowering the wave drag, their 
optimum distributions denoted by A, and B, are given 
by the negatives of A, and B,, respectively. In par- 
ticular, let us study Ay and B,, since, if we assume 
slender-body theory applies for the calculation of the 
body shape, these give the total body volume and lift. 
For simplicity, consider only closed bodies; then, if 
V’, and Lz are the volume and lift on the body, 


1 eh 
Vp = - > | xAy(x) dx (39) 
Leal ~<Es 
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BLp ie 
— B,(x)dx rT 
2q Ves =1 


Place Eqs. (35), (36), and the negatives of Eq. (3s 
in Eqs. (39) and (40), reverse the order of integration 
and one can show 





32 7 Ap Ap | 7 
lp —V,-+ , Zz dydx + y dsdx | 
Lp I E - 
where |", and L,. are the volume of the wing and th 


lift on the wing, respectively. 
Notice, from Eq. (41), that a lifting wing about th 


body creates an optimum multipole distribution whict 
will have, by itself, positive volume— that is, a real body 
can be constructed that will destroy all the equivalent 
source distribution caused by the wing so located 

Two examples illustrating some of the above con 
cepts have been worked out. In the first, a uniformh 
loaded elliptic wing having a biconvex section was 
moved above the body axis as shown in Fig. 14. Typi 
cal variations of body area and integrated loading 

(x) represents the lift integrated from the nose t 
station v| given by the optimum multipole distributions 
are also shown in the Figure. 

The effects on the wave drag of combining bodies 
having these area and load distributions with their 
wings are shown by the solid lines in Fig. 15. These 
results are simply a generalization of those presente 
earlier in Fig. 11. 

If the full value of B,(x) generated under the lifting 
wing were used to simulate the body shape, the total 
lift of the wing-body combination would be zero |set 
Eq. (42) ]. 
total lift, this distribution should be modified by the 
With such a 


restraint the wave drag can only be reduced to the 


Hence, if the airplane is to carry a fixed 


condition that the body carry no net lift. 


dashed line shown in the Figure, provided only A 
and B(x) are used to generate the body shape. The 
difference between the dashed and solid line is given bj 
Eq. (26) if So is zero. 

For the second example, an infinitely long semicircu 
lar cylinder, without thickness but with the loading 


Ap/q = 1/(Rk? + x 


was analyzed. In order to use this model to estimat 
the performance of a practical combination, the m 
finite wing was assumed to approximate a finite wing 
with the chord required to carry SO per cent of the total 
lift. Further, a body with a length equal to that of the 
optimum multipoles which the finite wing would gener 
ate was assumed to have the same relative effects 0 
the wave drag when combined with its wing as the 
infinite body would with the infinite wing. Fig. |! 
shows the drag of the wing alone and the wing plus th¢ 
various multipole distributions. For Cp, the total 


drag was used, for C,, SO per cent of the total lift 
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THEORY OF WI 


Both coetlicients were based on half the wing wetted 
irea. 

\vain the dashed line in the Figure represents the 
level to which the drag can be reduced by combining a 


vith a body generated by a(x) and the optimum 


wing W1 
distribution of 6;(x) under the condition that the body 
carry no total lift. 


{| no supporting strut was included, the study shows 


Although the wing has no thickness 


un 
that a very large portion of the wave drag due to both 
lift and thickness can be eliminated by the arrangement 
chosen 

A comparison of the partially shrouded body with 
two flat-plate wings is given in Fig. 17. Shown are 
the sum of the total vortex drag and total wave drag 
fa body and semicylindrical wing (without thickness 


together with drags of triangular and rectangular fiat 
plates. To actually achieve the results shown for the 
semicvlindrical wing and body, the wing should be 
cambered so that the vortex drag 1s a minimum and so 
the variation of lift required on the body is as small as 
possible. This accomplished, which does not seem un- 
reasonable in light of the results shown in Fig. 16, 
the comparisons are fair since all coefficients are based 
on the wing Note, the partially en 
shrouded wing-body combination has almost as low a 
BC,’ 
a considerable amount of volume. 


wetted areas. 


as the triangular plate even though 
For 


value of Cp 
it is carrving 
example, a 40,000-Ib. airplane flying at a Mach Num- 
ber equal to 2 at 50,000 ft., with an aspect ratio 2 wing, 
would carry an optimum volume equal to about 700 


cu. ft. This volume is sufficient to fill a Sears-Haack 
body 45 ft. long with a maximum diameter equal to 
(ft 


The Limiting Case M — 1 


After the original transonic area rule had been dis- 
covered, many pointed out that it was already con 
This 


was, in fact, true, but the significance of the linearized 


4 linearized theory. 


tained in the existing body 


results in this respect had been ignored or overlooked 


by everyone. However, this should not prevent us 
irom now reinspecting the linearized theory and seeking 
to find from it what further information we can about 
the flow in the transonic vicinity. 

Naturally, linearized theory as it applied to near 
sonic flow fields was not disregarded without reason, 
the reason being simply that the solutions to most 
boundary value problems contained whole regions 
filled with infinite perturbation velocities, a condition 
which clearly violates the assumptions basic to the 
theory. A few special cases were known which avoided 
this difficulty (as pointed out, for example, in refer- 
ence 14 Thus, lifting surfaces without thickness as 
well as infinite sweptback wings with constant airfoil 
section and chord have flow fields in which the per 
turbation velocities are small compared to a sonic free 
the 


existence of such flows was the condition that there be 


Stream velocity. Recognized as necessary for 


) net source strength in planes normal to the free 


NG-BODY WAVE 


DRAG lt 








—_"* 


16 





2.0 


12 
Bispan)®/ Plan form area 


Fic. 17. Sum of wave and vortex drags for shrouded body and 
two flat-plate wings 
stream. But an infinite streamwise cylinder having 


the same constant cross-sectional area before and be 
hind a wing and an indentation in the wing vicinity, 
such that the sum of the wing and body cross-sectional 
areas equals the same constant, also satisfies this cri 
terion, and perhaps in this sense the area rule is best 
described by linearized theory. To apply this result 
to practical flows the additional assumption (and a 
very strong one, but one that experiment now justifies 
must be made—that any area variation added to the 
stream surface comprising the wing and fore and aft 
cylinders to distort or close it yields a drag at near 
sonic speeds peculiar to that area variation only 

With these concepts in mind, linearized theory can 
the 


over the central portion of 


be used to estimate actual pressure distribution 
a thin wing and slender 
body combination in which the body has been modified 
according to the transonic area rule. Such calcula 
tions have been carried out by Byrd 
Further, results such as Eq. (21) can be examined for 
the case when 8 & 0 in an attempt to find simple ex 
the 


line have been carried out by Baldwin." 


tensions to area rule. Investigations along this 
Studying 
just the zero lift case, Baldwin expanded Eq. (4) in 
powers of 8 and found, at Mach Numbers near 1, the 
drag depends upon the longitudinal distribution of both 
the area and moments of area about the body axis 
Theoretically, this ‘‘moment of area’’ rule yields the 
same value of the wave drag as the area rule at .1/ l 
and, in addition, decreases the subsequent rise in drag 
as the Mach Number increases. 


Experimental verification of the usefulness cf the 











moment of area rule is shown in Fig. IS. The drag 
---- Experiment 
O16 — Theory 
“ Arec-rule 
_ 
) modification 
GB 
008. 
= Unmodified 
e 
Moment. of- area 
O rule modification 
8 
Fic. 18. Comparison between experiment and theory of wave 


drags for three different wing-body combinations 
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rise above .\/ 0.7 is shown* for three configurations, 
all having the same total volume. 
undistorted body (combination of Sears-Haack and 
Karman ogive shapes). The other two have been 
modified, but both have the same area distribution. 
By carrying some volume out on the wing panels in the 
form of pods, one of these has a better moment of area 
distribution and the resulting decrease in wave drag 


is evident. 
NONLINEAR TRANSONIC THEORY 


Although, as has been pointed out, linearized theory 
can be used to study a few very special types of bound- 
ary value problems at sonic speeds, solutions to tran- 
sonic small-perturbation problems in general should be 
be based on Eq. (2a) or (2b). The nonlinear char- 
acter of these equations makes the analysis of any but 
the simplest kind of problem exceedingly difficult. 
However, some progress has been made which has led 
to a few rather general conclusions. 

The following discussion of the analytical results 
based on the nonlinear Eq. (2b) is closely related to the 
well-developed theory of slender wing-body combina- 
tions furnished by Jones” and Ward.'* Therefore, we 
shall first review briefly the pertinent portions of 
slender-body theory and then describe its generaliza 
tion to the nonlinear case. 

First, we recall that the perturbation potential for the 
flow field in the vicinity of a slender configuration is 
given by the relation 

g(x, ¥, 3) g(x; y, 3) + f(x) 15) 
where ¢) is the two-dimensional harmonic potential 
satisfying the boundary conditions in transverse or 
crossflow planes. The function /(«) represents an addi- 
tive effect that is determined completely by the longi- 
tudinal distribution of cross-sectional area and _ this 
can be calculated from knowledge of the flow about an 
equivalent body of revolution at zero angle of attack. 
For supersonic flow, /(«) is given explicitly by the rela- 


tion 
P Oh Of 2. 2(x — 4%) 
f(x) = - S’(x,) In dx, (44) 
27 Ox Jo 6 
where the body nose is at x = 0 and S’(x) is the gra- 


dient of cross-sectional area S(x). Implicit in the 
theory is the fact that S’(x) has no discontinuities along 
the body and that its magnitude is small where body 
length ZL is the reference dimension. Furthermore, 
S’(0) = 0, but no further restrictions are placed on 
S’(L), the gradient of S(x) at the tail. 

If D denotes the drag of an arbitrary wing-body 
configuration and Dz denotes the drag of the equivalent 
body of revolution, the two values of drag may then be 


related as follows 


*In each case the theoretical or experimental value of the 
body alone wave drag has been subtracted out 


AERONAUTICAL SCIENCES DECEMBER, 


One has a simple 


19 5¢ 


p E O¢e ‘ 
D Dp = (tap Gs = 
2 be On 


[ O¢ 
aa ita On 


In Eq. (45), po is free-stream density and ¢» and 
are the crossflow potentials associated with the wing 
body configuration and the body of revolution, resp, 
tively. The curvilinear integrals extend around th 
contours of the two configurations in the plane of th 
tail. As written, Eq. (45) contains contributions t 
both the vortex drag and the wave drag, but it is not 
difficult to evaluate the vortex drag separately and t 
show that it is proportional to the square of the angk 
of attack and is represented by the same expression at 
all Mach Numbers. 

The wave drag of the equivalent body of revolutior 


is given by 


P= if) 2) pe Ui2(-( 24) S SS" (x1) S"(x2) In |x, - 
xe] dxydx2 + [S'(L)/r] SP S"(x1) In (L — x)dx; - 
1S(L)]2/2e} In [BV S(L)/2V w]) (4 


We observe here a dependence on Mach Number for 
bodies having finite angles at the base, and, in fact 
the theory predicts excessive values of drag as Mach 
Number approaches 1. On the other hand, when the 
class of bodies is such that S’(L) 
independent of Mach Number and remains finite 
throughout the Mach Number range. Although, as 


has been mentioned, the predictions at sonic speed ar 


0, the drag becomes 


certainly suspect by virtue of the infinite pressures 
the theory appears to be telling us that only for certain 
classes of shapes does a slender wing-body combina 
tion and its equivalent body of revolution have the 
same wave drag. We shall now see that this con- 
clusion is supported by the more exact transonic theory 

In a discussion as to what extent nonlinear theory 
modifies the results of the slender-body analysis just 
outlined, the papers of Oswatitsch,'® *’ and Harder and 
Klunker*' should be consulted. The following re 
marks are, however, drawn largely from a forthcoming 
paper by Heaslet and Spreiter. 

In order to achieve some success with the nonlinear 
equation of flow, we restrict our attention to a com 
parison of wings alone and bodies alone. For sufi 
ciently slender configurations it is found possible to 
derive from Eq. (2b) results directly analogous to those 
based on slender-body theory, provided the angle ol 
attack of the wing is kept small relative to the thick 
ness-chord ratio of the wing. Under these conditions 
it can be shown that the formal expression of Eq. (4° 
remains the same and that the perturbation potential 
in the vicinity of the wing is given by 


glY, ¥, 3) GAX, V, 3) T ZX 


where g(x) is determined by the longitudinal distribu- 
tion of cross-sectional area and, thus, can be calculated 
from a knowledge of the flow field about a body ° 
revolution at the same Mach Number and at zero angle 
of attack. In the case of slender-body theory the 
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Pheoretical wave drag 


could be written in terms of airplane 
(44) the 
no explicit result for g(v) has been determined. 


function f(A 
geometry —-e.g., Eq. but in nonlinear case 
How 
ever, one can show that Eq. (45) still applies in the 
nonlinear case where Dx is now the drag of the equiva 
lent body of revolution at sonic speed. This result 
orroborates slender-body theory in stating that in 
certain cases the drag of a configuration and its equiva- 
lent body of revolution are the same but that this is 
nly true when certain end conditions are satisfied. 
Although the function of g(v) in Eq. (47) cannot be 
given explicitly, it can be evaluated in certain specific 
ises. For example, the solution given for the circular 
one-cylinder by Yoshihara’ does furnish sufficient in- 
formation to calculate g(«) for cne particular variation 
f area. Such calculations have actually been carried 
The 


end conditions in these cases are such that the integrals 


ut for the case of slender elliptic cone-cylinders. 


i Eq. (45) do not vanish, and the results show that the 
drag of various cone-cylinders having the same longi- 
tudinal distribution of cross-section area, but different 
The 
magnitude of the effect is illustrated in Fig. 19, which 


cross-section shapes, may differ substantially. 


shows the variation of drag with ratio of major to minor 
axis of the elliptic cross section for two families of 
cone-cylinders. The cross-section area of all members 
of one family is the same as that of a circular cone with 
semapex angle w equal to 0.10 radian, whereas that of 
the other family is the same as a circular cone with 
# equal to 0.1225 radian. It can be seen that, within a 
sngle family, the drag is a maximum for the circular 
cone-cylinder and diminishes continuously as the cone 
and wider. For the two families 


becomes thinner 


shown, the drag of wing-like members (aspect ratio 
bout 2) is about 3/4 that of the circular cone-cylinder. 
Since solutions for sonic flow around three-dimen- 
sonal shapes other than cone-cylinders have not yet 
een obtained, the remainder of the discussions will 
9€ confined to cases in which experimental information 


Sknown for two or more wings or bodies having the 


at 


W l 


major axis 
minor OxIS 


ol elliptic cone-cvlinders 
same, or affinely related, longitudinal distribution of 
the 

) 


of data of this type is that given in references 25 and 


cross-section area. Probably most extensive set 
24 for a large family of affinely related wings of rec 
tangular plan form having NACA 63AOXX sections. 
Since the results for .\/ | can be presented most con 
cisely by using the variables suggested by the transonic 
similarity rules, the experimental results for the zero 
| are given in Fig. 20 by 
A is the 
As 


can be seen, the results define a single curve which ap- 


lift pressure drag at .l/ 


: ’ /3 - : ° 1/3 
Jotting Cp/r as a function of Ar where 
] £ 


wing aspect ratio and 7, the wing thickness ratio. 


proaches a straight line through the origin for low aspect 
ratio wings and a horizontal line for high aspect ratio 
wings. 

Now, if it is assumed that the effects of the violation 
of the theoretical requirements that occurs at the round 
nose of the unswept leading edge of each wing of the 
present family can be disregarded, and that Eqs. (47) 
and (45) can be applied, it can be shown that all wings 
having a given longitudinal distribution of cross-section 
area have the same drag as stated by the transonic area 
rule. Combination of this fact with the transonic simi 
larity rules* yields the result that Cp/7r’ 


5/3 


is linearly 
proportional to Ar’”. the data of 
Fig. 20 shows that the experimental results display 


2 


Examination of 


this trend for wings of Ar’® less than about unity. 
This limit to the range of applicability of the area rule 
is undoubtedly different for wings having plan forms 
that are nonrectangular, but lack of data precludes 
making corresponding comparisons for other plan forms 


at the present time. 
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Spectral Theory of Buffeting and Gust 


Response: Unification and Extension 


H. S. RIBNER* 


l’niversity of Toronto 


SUMMARY 


The statistical approach to the dynamic response of aircraft 
turbulent air is generalized and extended The basic element 
the response of the airplane to flight through an inclined 
jnusoidal upwash pattern in (a) three, (b) two, or (¢ without 
clination--one dimension. Case (a) is general; in case (b 
iirplane is idevlized as a lifting surface or a lifting line; in 
c), as a lifting point Earlier results of Liepmann are re 


vered in the cases of a lifting line and a lifting point. An ex 
nsion of the case of a lifting line is made for the sweptback 
fting line, and an extension of the case of a lifting point to 
vide rolling moments due to spanwise upwash gradients 


i) and the lifting surface in case (b) are new 


r 


SYMBOLS 


tangent of sweep angle (lifting line 

VV —-! 
k wave number vector (), Re, & 
implitude of single-wave lift response 
resultant gust lift response 
flight speed 
time 
unplitude of single-wave upwash 
resultant gust upwash 
coordinates relative to stationary axes 


coordinates relative to axes attached to moving 


airplane 
transfer function describing airplane single-wave 
response 


¥ gust spectrum function 


Subscripts (on scalar symbols 


gust idealized as one-dimensional 
2 gust idealized as two-dimensional 


Primes denote spectral quantities associated with Ow /Oy 


INTRODUCTION 


: PURPOSE of this note is first, to provide a gen- 
eral formulation for flight through turbulent air in 
three dimensions—-one that will include and extend the 
basic relations pioneered by Liepmann! for one dimen- 
sion and extended by him? and by Diederich*® to two 
dimensions and second, to introduce for all three cases 
unified spectral 
problems in 


approach, adapted from certain 


turbulence.? formalism 


iS derived in the main from reference 5, but the origins 


The spectral 
cited therein and in reference 4) go back much farther. 


Received February 28, 1956 
Research Associate and Assistant Professor, Institute of 


NVSics 


THREE- DIMENSIONAL GUST RESPONSE 


Lifting Airplane 


The approach is based on the concept that turbulent 
gusts and wakes may be represented as a superposi 
tion or spectrum of plane sinusoidal waves of shearing 
motion of all orientations and wave lengths; this is the 


Fourier 


o 
S* 


physical interpretation of a integral —e 


see reference 5. Considering the upwash component 
alone, an elementary spectrum component may _ be 


written 
awe ne l 


This is a plane wave pattern of arbitrary orientation 
in space: the planes of constant upwash are normal 
to the vector Rk ki, ko, k 
soidally in the direction of Rk with wave number k 
Viki t+ ke +k 


The incremental lift response of 


The upwash varies sinu 


(wave length = 27 k 
an airplane flying 
through this inclined wavy upwash field is of the form 


dLe™™" 2 


when the airplane is assumed to move along the \ 
axis; v is the coordinate of a reference point in the air 
plane. Egqs. (1) and (2) refer to stationary axes; rela 
tive to axes (xv’, y’ 
(2) implies a sinusoidal time variation of lift via x 
ve’ — Ut. 


Evaluation of the lift response |Eq. (2)] serves to 


moving with the airplane, Eq. 


determine a_three-dimensional ‘‘transfer 


I'(k) defined by 


function” 


dL = V(k) dW 3 


The transfer function plays a central role in gust re 
sponse. 

The spatial mean square upwash in the spectrum 
wave (more exactly, group of waves in band dk,, 
dks, dk;) is | dV 
spatially uniform, on the average) this may be written 
Wi Ridkidkodk; where Y(k) is a 


mean square upwash. 


-if the gust turbulence is homogeneous 


“spectral density” of 
(The statistical relations herein 
are exact for homogeneous turbulence and approxi- 
mate for atmospheric turbulence and wakes.) The 
corresponding mean square lift increment is (dL 
Thus upon summing over all of the waves that con 
stitute the turbulent gust pattern, there results 


mean square upwash: 


w? = LdlV 


SS Sv R)\dkidkodk }) 
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Fic. 1. The basic aerodynamic problem for gust response 


flight through a yawed sinusoidal upwash pattern. A closely re 
lated problem is flight through a vawed sharp-edged gust 


mean square lift: 
ee i I(R W(R dkidk dk (0 


in ky, Ro, R3. 


The first of these expressions relates the mean square 


L? = d\dL 


where the integrals go from © to a 


upwash in the turbulent field to its spectral density in 
three dimensions, ¥(Rk). The second evaluates the mean 
square lift increment in terms of the spectral density 


and transfer function. The two together relate the 


mean square incremental lift to the mean square up 


wash velocity in the gust pattern. A prerequisite is 


knowledge of the gust spectral density —e.g., see refer 
ence 2—- and the airplane transfer function. 
Eqs. (1) to (5) describe the lift response of a complete 


airplane flying through a random gust upwash pattern 
This is the most general formulation of the 
(The 
negligible contributions of w and wv fluctuations to lift 
this generality is the 


in space. 
gust-response problem in the engineering sense. 


are excluded.) The price of 
difficulty in calculating the transfer function I'(k;, k», 
k;). The difficulty 
along with the generality, by idealizing the airplane as 
or (3) a lifting 


may be reduced progressively, 
(1) a lifting surface, (2) a lifting line,* * 
point.! 
fwo-DIMENSIONAL GUST RESPONSI 

Lifting Surface 

As a lifting surface the airplane is assumed to be 
sensitive only to the pattern of upwash in the (x, y 
plane; the vertical (s) variation, of wave length 27 ks, 
is taken to have no effect on response. Thus ky may 
be integrated out in the equations with a consequent 
reduction to two dimensions, the (x, v) plane. The re 


sults are 
elementary wave: 


dWee** + ™ (6) 
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lift response 
d Live 7 
transfer relation 
dL. = V(ky, kx)dW ¢ 
spectral density 
Wolk, ke JS WkRidk 
mean square upwash 
= LidW. SS velki, Re\dkdi 
mean square lift: 


S SV (Ri, be) |? Wodkidk 


The upwash pattern in the elementary two-dime 


LL? = Zid. 


sional wave [Eq. (6)| when plotted as a relief m 
has the appearance of a corrugated sheet. The cor 
rugations are yawed relative to the spanwise (y) dire 


(the x axis has beer 
taken as The 
namic problem is thus the calculation of the lift 1 


tion by an angle arctan ky ky 
basic aerod) 


directed downstream 


sponse |Eq. (7)] of the lifting surface flying through 
yawed sinusoidal gust [Eq. (6)]. (See Fig. 1.) This 
suffices to determine the two-dimensional _ transfer 
function I'(k;, k2) by Eq. (8S). The calculation must 
be performed for all values of k; and ky». 
Alternatively, the basic problem may be taken 


the lift response AL(k» k,, &) for flight through a yau 
sharp-edged gust: an extension of the notion of 
dicial response Here tan! (—k» k;) is the angle 


yaw of the gust leading edge, and & is the distance | 
penetration of the reference point into the gust, me 

The lift for th 
obtained 


ured along the —w axis. respoise 


sinusoidal gust of Eq. (6) can be therefro 


by means of Duhamel’s integral —e.g., see reference 6 


Lifting Line 


If the airplane is still further idealized as a lifting 


line the instantaneous chordwise variation of upwas 
is neglected. Under this assumption Liepmann’® an 
Diederich’ have formulated the transfer function for 
wing without sweep in terms of a simple strip theor 
for wing response. Under the lifting line assumpti 
(S) lead directly to the same result [EF 


1S easull 


Eqs. (6) to 
Moreover, the result 


it is found that k» in EX 


(14) of reference 2]. 
generalized to swept wings: 
(14) of reference 2 need only be replaced by k 


+ dk 


where a is the tangent of the sweep angle of the refer 
ence line, and the plus and minus signs are applied 
respectively to the right- and left-hand wing panels 
The derivation is given in the Appendix. It is of 
terest that k» + ak, 
matches the wing-panel sweep. 


vanishes when the wave yaw 


ONE-DIMENSIONAL GUST RESPONSI 


Lifting Point 
The limiting simplification represents the airplane 


a lifting point moving along a line. This point sees 


' 
| 
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nly the upwash pattern along the x axis; the spanwise 


variation, of wave length 27/k», is taken to have no 
effect on lift As was done with k;, ko may be inte- 


rated out in the equations with a consequent reduction 
The results are 


ty) one dimension: the x axis. 


elementary Wave 


dW,e' ) 
lift response 
dLy } 
transier relation 
dl P(k,)d WV, 14 
spectral density 
valk SL” Wolk, ke)dk 15 
or 
i(k; SS vik dkodk 6 
mean square upwash 
>i dW, S Wiki dk, re 
mean square lift 
L >| dL, S Uk vidk S 


Rolling Moments 


Spanwise upwash gradients would seem, at first 
thought, foreign to the one-dimensional idealization of 
the last paragraph. However, they may be injected 
in first approximation by replacing the upwash field 
by the first two terms of its Taylor expansion 1n y: 


: yiOw rehy x0 


w(x, vy) ~ w(x, O) 


Eqs. (12) to (18) refer to the first term. A similar set 


f equations corresponding to the second term are 


where the prime does not signify a derivative 


elementary wave: 


(Ow/OY)(x0) = dW,'e" 19) 
rolling moment response: 
dle" 20 
transfer relation 
dl, = T'(k,)dW,’ 2] 
spectral density : 
Vi'(ki) = S’he*polki, ke)dke 22 
or 
Wi'(ki) = SS keeW(k)dkodk 23 
mean square Ow /Ov: 
‘te 2 j mi — . 
a) >| dw, e Wy (k; dk, 24 
mean square rolling moment: 
>| dl; = S I'’(ky) |?Wy'dky 25 
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The effects of vertical (s) gradients may be evaluated 
similarly, with k; replacing k. as the factor of y in the 


analogs of Eqs. (22) and (23 


CONCLUDING REMARKS 


The present formalism may be enlarged to include all 


three components u, v, w of gust velocity, the corre 
sponding forces, moments, and deflections developed 
by the airplane, and, finally, the dynamics of the re 
sulting airplane motion, with considerations of aero 
The analyses of references 


elasticity. 3 and 7 point 


the way. 
SWEPTBACK 


PRANSFER FUNCTION OIF 


LIFTING LINI 


\PPENDIX 


the increment of wing lift at time 7 
angle of attack at time O over 


Let h(r v be 
following unit impulsive 
t chordwise strip of width dy at y. Then, if the angle 
of attack variation along the reference line in the wing 
(’, the corresponding wing lift is 


eor( f 
IS WU 


| b/2 
lift dy A(t, Vv) wil r. 
i eg ; 


Now consider flight with velocity l’ through the sta 


tionary sinusoidal wave, [Eq. (6)|, such that the air 
plane moves a distance £ in the — x direction in a time 
7; then time and distance are related by 7 E/U. At 
time / the upwash distribution along the swept reference 


line in the right-hand wing panel may be written 
dWe' 


wt, y 


where x = & + ay defines the sweptback line in terms 


of (7, 0), the coordinates of the apex. At the earlier 
time / — 7 the line is a distance & farther downstream, 
whence 


dWeeit(® + & + an) + 5 


Wi — 7, J} = 
This relation may be generalized to apply to both wing 


panels and inserted in Eq. (26), the lift may be written 


dL.e"*, and the ratio dL. dW. = Ik. k.) formed, 
with the result 
I'(ky, ko) = 
| eh/2 . . t oe ' 
dy h eae . dé 27) 
U2. i * a ss 


wherein the plus sign in the exponent refers to the right 
hand wing panel, the minus sign to the left-hand wing 
panel. When the tangent a of the angle of sweepback 


») 


is zero, Eq. (27) reduces to Eq. (14) of reference 2. 
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Some Integrated Volume Properties in 
Linearized Flow and Their Connection with 
Drag Reduction at Supersonic Speeds 


Z. O. BLEVISS* 
Douglas Aircraft Company, Ine. 


SUMMARY 


In some recent linearized studies of the drag of generalized 


(nonplanar) aircraft configurations in supersonic flow the bodies 


(aircraft configurations) were represented by singularities (ele 


ments of lift, side-force, and volume This paper studies one 
aspect of the problem of relating the singularity distribution to 
the body geometry—-namely, the relationship between the total 
singularity strength and the total body volume 

The “body”’ 


mean surface and the stream surface that coincides with it at up 


studied here is the region between a evlindrical 


stream infinity, the mean surface being parallel to the free stream 
and having an arbitrary closed cross section. It is shown that 
(a) in the absence of transverse forces the total body volume is 


equal to the total strength of the volume elements within the 


mean surface, (b) the body volume due to transverse forces may 
be obtained from an integration of the potential in the Trefftz 
plane, and (c) for certain mean surface cross-sectional shapes and 
positions the body volume due to given transverse forces is equal 
in magnitude to the fuselage volume required along a given axis 
to produce the minimum wave drag. 


INTRODUCTION 


I REFERENCES | and 2 the drag of generalized (non 
planar) aircraft configurations in supersonic flow 
was studied using linearized theory. It was assumed 
there that the body (aircraft configuration) could be 
represented by singularities (elements of lift, side-force, 
and volume). It is obviously necessary to relate the 
distribution of the singularities to the geometry of the 
body. One aspect of this problem will be studied here 
namely, the relation of the total volume of the body 


to the total strength of the singularities. 


DISCUSSION 


In principle, a body shape may be constructed from 
the singularities as follows. Take a cylindrical mean 
surface with generators in the free-stream direction; 
in cross section the mean surface is closed and of arbi- 
Consider the stream surface which at up- 
The 


“body” will be the region between the mean surface 


trary shape. 
stream infinity coincides with the mean surface. 


and its associated stream surface between upstream and 
downstream infinity. The volume will be counted as 
positive when the stream surface is outside the mean 
surface and negative in the opposite case. The stream 
Received December 13, 1955 
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theor 


outwas| 


surface should be computed from /inear 


This means that if v, is the normal velocity 
at the point P on the mean surface and a is the slo 
of the stream surface 


(local angle of attack it a poir 


Q such that PQ is normal to the mean surface, the: 


al’, = vn, where l 
body shape will thus depend on the mean surfa 


chosen. It will in general be hollow. 


will not be hollow in the two limiting cases most gener 


ally encountered 
bodies the singularities are located along a line, 


mean surface is a circular cylinder whose radius tends 


to zero; (b) in the theory of thin planar wings t! 


singularities are distributed over a plane, and the mea 


surface also shrinks to a plane in the limit 


venience, the bodies obtained in these limiting cases 


will be called ‘‘local bodies.”’ 
The volume element! 
its axis in the free-stream direction. 


of volume elements is, of course, equivalent to a dis 


tribution of sources and sinks. 


The purpose of this paper is to prove the following 


statements. (1) The volume of a body is equal to tl 


total strength of all volume elements inside the meat 


surface, provided no other singularities are 
(2) The volume of the body is the same for all me 
surfaces containing the same volume elements, provide 
no other singularities are present. [Note that 

follows immediately from statement (1).] (3) TI 
contribution of other singularities—i.e., transverse-fort 
elements 
potential in the Trefftz plane according to Eq. 

The significance of this contribution will be illustrate 


to the volume can be evaluated from tl 


later. 
Before proceeding with the proof of statements 
and (3) it is important to note that since the body shaj 


is constructed from a linearized mean surface as 


scribed above, the volume is additive in the singular 
Hence it is sufficient to consider any simp! 
It is conve! 


ties. 
distribution of singularities in the proof. 


ient to consider the volume elements distributed alon 


a streamwise line, but the distribution of transvers 
Then, in the pro 


the volume elements may be replaced by a local be 


force elements may be general. 


of revolution, according to the well-known slende! 
body relation between source-strength and _ cross-st 


tional area. Higher order singularities, such as qua 





is the free-stream velocity. Th 
However, j 


(a) in the linearized theory of slender 
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rupoles, ete., do not contribute to the volume and will 
The the 


proof will be more general than statements (] 


results obtained in following 


ri ignored 


and (3 
Consider the control surface and local be dy of revolu 


tion shown in Fig. | The control surface has the fol 


wing components: S: is the mean cylindrical sur 


face (with arbitrary cross-sectional shape) extending 


from ¢ to & v; 5S; 1s normal to the free 
stream direction at & -o and is equal to the 


ross-sectional area of the mean surface; .S3 1s normal to 


the free-stream direction at & v and is equal to 


S h where S;x, 1s the cross-sectional area of 


the local body at &€ = x; 5S; closes the control surface 
following the surface of the body. 

For steady flow, the conservation of mass requires 
that the mass flux out of S. and S; be equal to the mass 


flux into .S;—1.e., 


where p’ 1s the density perturbation, uw is the velocity 
perturbation in the the 
velocity (perturbation) at and normal to the mean sur 


free-stream direction, 7, 1s 


lace, and the subscript denotes free-stream condi 


tions. Eq. (1) is general. In subsonic and supersonic 
linearized flow the perturbations are small and the 
products of perturbation quantities may be neglected 


Eq. (1) then reduces to 


where the dependence on x of certain quantities is ex 
plicitly exhibited for clarity. the fact that 
p p...M.*u, the free-stream 


Mach Number, Eq. (2) may be written 


[{ iia heey’ = — 1 | ait. i 


lhe inte gral over S» can be interpreted as the cross 


Using 
where .\/ is 


Ss 


sectional area S, Fig. | 


Briefly, 


of the body at & x (see 


this can be seen as follows. At any point on 


So, v,/l 


is the slope of the stream surface in the 


Jane containing v, and Integrating this slope 
~ 


along a generator of Sy from & to & x yields 
the body thickness at & Integrating this thick 


ness around the mean surface cross-sectional contour 


become Ss 


M * ig 
= fms a 
| 


gives Sz(xv) and Eq. (3 


Spa Srp 


where the sign convention for Sz, is the same as for 
see Fig. |] 
The volume of the body from ¢& to & vis 


| \ S fidé 


Interchanging the S and € integrations and noting that 


“uaé Q 


where ¢ is the perturbation velocity potential, there re 


. Vu rc 


I,» is the volume of the local body up to \ rhe 


sults* 


sign convention for |’, is the same as for ¢ 

The case of particular interest in this paper is 4 
to, Then, in Eq. (5), V, is the body volume 
between upstream and downstream infinity, [77% 
sum of all vol 


is the volume of all local bodies (or the 


ume elements) contained within the mean surface, and 
the last term is the integral of @ in the Trefftz plane 
over the area enclosed by the mean surface contour. 


If the flow field is created by volume elements alone 


i.e., there are no transverse forces—@ 0 in the 
Trefftz plane and I’z = I Chis proves statements 
1) and (2). 

Assume now that, in addition, transverse-force 


In general, @ + 0 in the Trefftz 
Then 
represents the contnbution of 


elements are present. 


plane and the integral in Eq. (5) 1s not zero. 
the last term in Eq. (5 
This proves 
the 


In particu 


the transverse forces to the body volume 
This 
mean surface chosen and may be negative 


lar, it should be noted that the contribution approaches 


statement (5 contribution depends on 


zero as the area S; approaches zero. Thus, transverse 
forces do not contribute to local body volumes 

To illustrate the above remarks consider the axially 
ZeToO 


symmetrical supersonic flow about a fuselage at 


angle of attack surrounded by a ring wing (see refer 
ences | and 2). Let the ring wing have zero thickness, 
infinitesimal Ac, 


loading per unit length around the wing 


chord and uniform outward radial 


re must be exercised 


In this interchange of integrations ca 
whenever the £ integration intersects the body upstream of 
rhis leads to an additional term in Eq. (5) which may be neglected 
in the linear theory However, for nonslender bodies in inc 
pressible flow (.1/ (}) Eq. (4) applies when Sy is far enough 
from the body, and, in general, the additional term in Eq. (5 
must be retained. For is case the additional term is equ ) 
the ‘‘volume of additional apparent ss’’ when .S; is downstre 
of the bod, 
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—e nN 


MEAN SURFACE — 


TREFFTZ PLANE 
(x = + @) 





F = Cp(1/2)p.U 7h 
The potential distribution in the Trefftz plane is 


j 0 outside the ring 


o= . ' — meee: 

(\—F/p.U. = —(1/2)CrlU..Ac_ inside the ring 
Consider any mean surface enclosing the ring wing. 
The contribution of the fuselage to I’, is the fuselage 


volume Vr. The ring wing contributes a volume 


Vew = —(1/2)Cr(.VW.? — 1)rR?Ac 


where F& is the radius of the wing. The total body 


volume is 
lp ong Vy + Ven 


A necessary condition for zero wave drag is I’, = 0. 
A necessary and sufficient condition is that the body 
volume be zero identically at all points of the mean sur- 
face—i.e., V, is zero at every point of the mean surface. 
Since the ring wing contributes negative volume it is 
clear that |’, can be made zero for a given wing loading 
F by adjusting I’y. Furthermore, the fuselage volume 
can be distributed along the axis of symmetry in such a 
‘way that I’g(x) = O for all x. Using the axial sym- 
metry, this means that the body volume is zero identi 
cally at all points of the mean surface and, hence, the 
Since there is no trailing vortex 
This is the ‘Ferrari 


Wave drag is zero. 
system, the total drag is also zero. 
ring wing plus central fuselage.” It is the three- 
dimensional counterpart of the Busemann biplane. 

In the general case, the physical significance of the 
volume contributed by transverse forces is not known. 
However, for certain cases in supersonic flow a physical 
significance can be attached to this volume. The 
above example is a particularly simple case that illus 
trates most of the essential points. It is seen that, if 
this volume is negative, a fuselage may be introduced 
to produce favorable interference with the transverse 
forces and, hence, to reduce the wave drag. The fuse 
lage volume required for minimum (not necessarily zero 
wave drag is equal in magnitude to this volume. This 
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fuselage volume must be properly distributed to att 
the minimum. 

In the more general case of a given transverse-for 
element interfering with a fuselage that lies along 
given axis, the correct result for the fuselage volume ; 
quired for minimum wave drag is obtained from f 
»)) if a particular mean surface is used. Very brief 
the result is obtained as follows. Consider a liftj; 
element (force in +2 direction) with infinitesimal cho; 


Ac and span As carrying a load 


Choose a mean surface with circular cross section tha 
encloses the lifting element and whose axis of revoly 
tion coincides with the given fuselage axis. The g 
ometry in the Trefftz plane is shown in Fig. 2. Th 
potential distribution in the Trefftz plane corresponds 
to the two-dimensional incompressible flow due to 


dipole (see, for example, reference 2) and is given by 


@ = (1 4n)CrU.AcAs[z/(y? + 2? 


To obtain the body volume due to the lifting element 
@ must be integrated over the circular region enclosed 
by the mean surface. It can be shown® that the in 
tegral is zero for the annular region between the mea: 
surface and the concentric circle (dashed in Fig. 2) that 
passes through the dipole. The integral of ¢ over th 
area enclosed by the dashed circle can be readily ol 


tained by using the mean value theorem of potential 





theory,* which for this case can be stated as follows 
if @ has no singularities within the dashed circle th 
average value of @ over the area enclosed by the circk 
is equal to @ at the center of the circle. Then the in- | 
tegral of @ over the area enclosed by the circular meat 
surface is the value of ¢ at the fuselage axis multiplied 
by 7?, and the body volume due to the lifting element 


is {with s (vy? + 22) = —[(cos y) R]! 
| 
I ae [(./ — | t]/CrAcAsR cos 7 


Setting the total body volume equal to zero— 
ly = Ve + Vie = O—the required fuselage volun 
This is the result given in references 
5, 6, and 7 for the fuselage volume required for mim 
mum wave drag. It should be noted that, since th 
flow field is not axially symmetrical about the fuselag 
axis, the minimum wave drag cannot be zero; henc 
the body volume is not identically zero at all points 0! | 
The significance of the body vol- | 





the mean surface. 
umes due to the lift element for mean surfaces with other 
cross-sectional shapes is not yet clear. A more C0! 


plete discussion of these results is given in reference 
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Influence of Aerodynamic Heating on the 
Effective Torsional Stiffness of Thin Wings 
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SUMMARY 


of the loss of torsional stiffness that may be 


subjected to axial stresses induced by 


Theoretical results are presented for the 


effective torsional stiffness of solid wings of symmetrical double 


cro section 


iccelerated to supersonic speeds. It its 


ted that serious losses of torsional stiffness may occur 


onstr ed 
wings 
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uti 


the influence of various parameters such as 


celeration, and thickness ratio are discussed 


SYMBOLS 


cross-sectional area 

quantity defined by Eq. (17 

Young’s modulus 

shear modulus 

function defined by Eq. (23 

mass polar moment of inertia per unit length 

torsional stiffness constant 

Mach Number 

Prandtl Number 

function defined by Eqs. (22 

Reynolds Number 

temperature 

total torque on cross section 

torque on cross section due to axial stresses 

torque on cross section due to St. Venant shear 
stresses 

complementary energy density 

velocity of sound 

span 

chord 

specific heat of wing material 

maximum thickness of wing 

acceleration due to gravity 

heat-transfer coefficient 

heat-transfer coefficient based on final fligh 

M 

parameter defined by Eq. (B-17 


velocity 


thermal conductivity of air 

function defined by Eq. (27 

function defined by Eq. (B-22 

radial distance from center of twist to point in 
cross section 

Av A 

wing thickness 

skin thickness of hollow wing 

rectangular coordinates 

dummy variable 

function of time parameter A defined by Eq. (9 

coefficient of thermal expansion 

coordinate, 2x /t 


nondimensional chordwise 


h constant along the chord 
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Superscripts 


ratio of specific heats 


shear strain in xy-planc¢ 


spanwise strain 
nondimensional spanwise coordinatc, 24 
temperature-recovery factor 
temperature parameter, | 7 ] 1 
1 or 1 7 1 1 if 
1 = 7 
time parameter 
value of time parameter when 1/7 = 1/ 
kinematic viscosity 
nondimensional chordwis¢ coordinate, 4 


h variable along the chord 
mass density of wing material 
normal stresses 
time 
time at which JJ, is attained 
shear stress in xy-plane 
angle of twist 
complementary energy 


frequency 


average 

adiabatic wall 
effective 

elementary 

final flight condition 
maximum 

minimum 

initial flight conditior 


free stream 


initial flight condition 
final flight condition 


INTRODUCTION 


’. EFFECT of aerodynamic heating that has received 
little attention is the apparent change in the 


structural stiffness of aircraft components that can be 


caused by thermal stress. 


In the present paper at 


tention is centered on the loss of torsional stiffness that 


may be incurred by thin wings subjected to thermally 


induced axial stresses 


The physical mechanism under 


lying this tvpe of stiffness loss is reviewed, and theoreti 


cal results are presented for the effective torsional stiff 


ness of solid wings of double-wedge cross section ac 


celerated to supersonic speeds 


Details of the analyses 


are given in appendixes. 


EFFECTIVE TORSIONAI 


STIFFNESS 


It is well known that the torsional behavior of beams 


having thin cross sections can be strongly affected by 


axial 


stress. 


For example, the torsional buckling of 


bars is attributed to the action of axial compressive 
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stress in effectively reducing torsional stiffness.' In 
an entirely analogous fashion, the presence of axial 
thermal stress in a thin beam can produce an apparent 
change in its torsional stiffness. 

The torsional effect of axial stress can be calculated 
simply. Fig. I(a) shows a longitudinal element of a 
beam, with axial stress of magnitude a, at some point 
P located a distance 7 from the axis of twist AA. If 
the beam is twisted, the fiber PQ, originally parallel to 
the axis AA, assumes the new position P’Q’ shown in 
Fig. l(b). It may be reasoned intuitively that the stress 
a, follows the fiber’? and hence assumes the direction 
P’Q’ when the beam is twisted. There is then a com 
ponent o,r(dg dy) of the axial stress that is perpendicu 
lar to the axis of twist. The total torque on the cross 
section due to the initial axial stress over the entire 


cross section 1s then 


fi. | o,r(de dy)dA 
JA 


This torque due to the axial stresses must be added to 
the torque due to the St. Venant shear stresses 


Ist. v GJ(dg/dy) 2 


to vield the total external torque on the cross section 


T, = g 1 | a rdA | (de dy 5) 
JA 


The quantity in the brackets is the effective torsional 


stiffness of the beam; thus, 
GT ess GJ = l on | o,r°dA GJ (4 
JA 


The same result has been derived by Goodier,’ without 
recourse to the intuitive approach, on the basis of gen 
eral equations of nonlinear elasticity. 

When the axial stresses ¢, on a free-end beam (such as 
a wing or stabilizer) are due to temperature, they are 
necessarily self-equilibrating; that is, the net thrust and 
bending moment due to ¢, must vanish. Consequently, 
as can be easily verified, the magnitude of the torque 
[see Eq. (1) | due to such thermal stresses is independent 
of the axis of twist, as is the St. Venant torque [see Eq. 
(2)].. It may be worthy of mention that, within the 
framework of elementary beam theory,  self-equili- 
brated axial stresses have no influence on bending stiff- 
ness. 

It can be seen from Eq. (4) that a distribution of 
thermal stress o, that is compressive (negative) in the 
vicinity of the leading and trailing edges of a wing would 
cause a decrease in the torsional stiffness. As will be 
shown, such distributions of thermal stress are produced 
by the aerodynamic heating encountered in accelerated 


flight. 


SoLID DOUBLE-WEDGE WINGS 


On the basis of various simplifying assumptions, the 
temperatures, stresses, and torsional stiffness changes 
in solid wing sections of symmetrical double-wedge 
profile [see Fig. 2(a)] undergoing various flight histories 
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are calculated in Appendix A. Before the detajk 
presentation of final results of some generality, numer 
cal results for a particular example will be shown in orde, 
to exhibit the essential features of the problem 


Illustrative Example 


Consider a double-wedge airfoil of solid stainless st, 
instantaneously accelerated from Mach Number zero t 
Mach Number 3 at 50,000 ft. altitude. The chord is 
in. and the midchord thickness——3 per cent of 
chord—is 1.08 in. On the assumption that the init; 
temperature of the wing is the stratospheric free-stregy 
temperature of 67°F., the temperature variatio; 
along the chord is shown in Fig. 3(a) for several values 
of time 7 following the acquisition of Mach Number: 
'n the calculations, the temperatures were assum 
constant through the thickness, and heat conduetio 
along the chord was neglected. The symmetry of th 
temperature distribution about the midchord is a consi 
quence of the additional assumption of a constant heat 
transfer coefficient in the boundary layer along th 
chord. The variation of temperature along the chord is 


simply a consequence of the fact that the massive mid 


chord region takes longer to heat up than the thinner 
leading and trailing edge regions. An elementary calcu 
lation of the axial stresses resulting from the tempera 
tures in Fig. 3(a) yields the distributions shown in Fig 
3(b) for the corresponding instants of time. Eventu 
ally, for very long times, the temperature in the wing 
would become uniform at the adiabatic wall tempera 
ture of 569°F., and, correspondingly, all thermal 
stresses would vanish. But at the intermediate times 
shown it is seen that compressive stresses occur near th 
leading and trailing edges, with equilibrating tensile 
stresses around the midchord. At each instant, then 
the torsional stiffness is reduced according to Eq. 
with the resultant time variation of effective stiffness 
shown by the intermediate curve in Fig. 4. Thus 
about 1.25 min. after the sudden acquisition of Mach 
Number 3, the wing section sustains its maximum loss 
of torsional stiffness having effectively only 22.5 per 
cent of its original stiffness left. As flight at Mach : 
continues, the stiffness gradually increases, tending 
toward its original value as temperatures become uni 
form and thermal stresses vanish. It should be men 
tioned that the shear modulus G has been assumed 11: 
variant with temperature; any reduction of G with 
temperature would, of course, lower the torsional stiff 
ness further. 

Also shown in Fig. 4 are the variations with time 
GJ,,;, GJ for instantaneous accelerations to Mach Nun 
bers 2. and +. The indicated reduction to zero stiffness 
after about 35 minutes of flight at Mach Number 
may be interpreted as torsional buckling. According t 
Fig. 4, the wing at Mach Number 4+ would unbuckle at 
time 7 = 2.3S min. and then gradually regain its tot 
sional stiffness. 

Since instantaneous acquisition of velocity is phys! 
cally imp ssible, it is of interest to determine how ther 


sults would be modified by the imposition of mort 
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INFLUENCE OF AERODYNAMIC HEATING LOSS 
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realistic flight conditions. Fig. 5(a) shows a family of Tl — T,,° — 
five flight histories: the instantaneous acquisition of FD ne Te ia 
Mach 3 already considered, and accelerated flights that O<x<c2 
reach Mach 3 in 1, 3, 6, and 10 min., respectively. The — 4 : 
calculated variation with time of GJ,;,G/J for each — 9 P 
flight history is shown in Fig. 5(b). It is of significance 2 SxS (6b 
that, at least for the particular wing considered, the i 
maximum loss of stiffness for moderately rapid ac EalT, = 1 - . : — (lL = Ade y 
celerations is only slightly less than that calculated for NEI r 0 c/2 (7 
instantaneous acceleration; even for the very sluggish F 
acceleration that takes 10 min. to get up to speed a ¢ i | A )e + 
stiffness loss of over 50 per cent would be incurred. NEU —X ce2SxSc (tb 
The next section presents the general results on which GJ.../GI | (Ea/G) X 
the preceding examples were based, and also discusses 117 = d/o) AA) (8 
the assumptions with which the results were derived. 
where 
GENERAL RESULTS AND DISCUSSION Ar) 1/8)fe7 (4X + 7A2 + D 
Eu(—2r)(10A* + SA*® + A 


Instantaneous Acceleration 


The following initial and flight conditions are speci 
fied: 

(1) For time 7 < O, the wing flies at Mlach Number 
Vo. 

(2) The initial wing temperature at 7 = OU is 7, 
the adiabatic wall temperature corresponding to .l/,; 
in general, adiabatic wall temperature is given bv 

Tou Vj 1 + [nly — IMP /2); 5 
where 7°. is the free-stream temperature. 

(3) At time 7 = O, the wing suddenly acquires a speed 
of Mach Number .1/;. 

(4) The free-stream temperature 7 
stant throughout the flight; this condition may be con- 


is assumed con- 


sidered to correspond to flight with any altitude varia- 
tion in the isothermal stratosphere between 35,000 and 
105,000 ft., or to constant altitude flight. 

The temperatures in the wing cross section are calcu- 
lated for time 7 > O on the basis of the following as- 
sumptions: 

(1) The heat-transfer coefficient which governs the 
flow of heat from the boundary layer to the wing is con- 
stant along the chord, but may vary with time (as the 
altitude changes, for example). 

(2) Temperatures are assumed constant through the 
wing thickness at each chordwise station. 

(3) Heat conduction along the chord, as well as span- 
wise heat conduction, is neglected. 

(4) Radiation effects are neglected. 

Finally, the axial stresses on the cross section are 
calculated c» the basis of the elementary theory that 
assumes plane sections remain plane, with no lateral re- 
straint, and torsional stiffness changes are computed 
according to Eq. (4). All material properties, such as 
elastic moduli and coefficient of thermal expansion, are 
assumed invariant with temperature, and all deforma 
tions are assumed to be elastic. 

It is shown in Appendix A that, for 7 > 0, the tem- 
peratures, stresses, and effective torsional stiffness are 
then governed by the following equations: 


for constant fh = h 10a 


-\) is the exponential integral 


(¢ s)dz 


An alternative form of the result [Eq. 


In Eqs. (7) and (9) Ea 
function defined by 


Eu(—A 
and tabulated in reference 3. 
(S)] for the 
effective torsional stiffness is obtained by the use of Eq 
(5) for the adiabatic wall temperature; thus, 


(GJ,5,/GJ) = 1 — (Ea/G)T 2.) (MP — Mo?]/(d/c)*} X 


A plot of A(A) versus J is shown in Fig. 6; the maximum 
value of A(A) occurs for \ = 0.386 and is equal to 0.0456 
Since \ is a monotonically increasing function of time 
this maximum value of A must be attained at some time 
during the flight. Consequently, the minimum value 
of GJ, 7,/GJ, from Eq. (8), is 


(GJ e57/GJ) min = 1 — (0.0456)(Ea/G) X 

[Taw — Taw |/(d/c)*} (12 
or, from Eq. (11) with 0.9 and y = 1.4, 
GJ.57/GJ)min = 1 — (0.00821)(Ea/G)T.. X 


;([M, - M bf (d ¢ 15 


From Eq. (13) it is seen that the maximum per cent 


reduction in torsional stiffness is dependent of 


parameters that may be characterized as pertaining t 


(a) material, Ea/G; (b) geometry, d/c; (ec) altitude 
[..; and (d) velocity, MW, — M,’. 
It is important to note that the result for 


GJ .¢;/GJ) min is not dependent on the magnitude of the 


heat-transfer coefficient. 


min 


This fact is actually a comnsé 





C2 (6a 
( 6b 

) " 

Cc = ‘é 
Nje * 4 
ib 
L(A X 
A J 
10a 


1 integral 


| for the 


se of Eq 


laximum 
0 0.0456 
of time 
yme time 
m value 


( 
\ 


per cent 
ent on 
ining t 
iltitude 


ult for 
e of the 


1 conse- 





INFLUENCE OF 


yence of the assumed infinite acceleration, but it is 


significant nonetheless, since, as will be shown, the re 
suits for infinite acceleration remain nearly valid in 
many cases of moderately rapid acceleration. 

It may be remarked that if interest centers on ac 
celeration from a subsonic speed to a moderately high 


supersonic speed, little error is introduced by setting 


v= 0 For example, if .\/) is actually 1 2, then l/, 
iffers from 7 — A/,? by less than 7, 3, and 2 per 
centfor J, = 2, 3, and 4, respectively.) On the basis of 
Eq 13), with 14, = 0, Fig. 7 shows how (GJ,;; GJ)» 
varies with the thickness ratio d c, at several values of 
for the case of stainless steel wings (a G = 24.2 X 
0-*/°F and stratospheric flight (7 393 °F. 
ibsolute). The Figure shows how susceptibility to loss 


f torsional stiffness increases sharply with decreasing 
thickness ratio. The effect of varying material is indi 
cated in Fig. 8, wherein (GJ/,;, GJ), 
V, = 3, for titanium (Ha/G = 11.4 & 10-*/°F.) and 
juminum (Ka G = 33.4 & 10~-°/°F.) as well as for 
steel. (These assumed values of a/G, incidentally, are 


is shown, for 


only to be regarded as reasonable estimates for illustra 

tive purposes; in any particular design problem, a care 

ful choice of the value of Ea G appropriate to the tem 
perature range involved should, of course, be made. 

A useful representation of Eq. (13) is afforded by Fig. 

, which shows (for steel wings in the stratosphere and 

0) plots of d c against .\/, for several values of 

GJ 4, GJ 

thickness ratios at any Mach Number; for example, a 


This Figure indicates the limitations on 


steel double-wedge wing cannot have a thickness ratio 
less than about 3! » per cent if it is to retain more than 
half its original torsional stiffness after instantaneous 
eceleration to Mach 8 in the stratosphere. 

fhe results previously shown in Figs. 3 and 4 for the 
illustrative example of a 3-ft. chord, 3 per cent thick 
steel airfoil at 50,000 ft. altitude were obtained from Eqs. 
‘in order to show the results in terms 


6), (7), and (11 
of time in minutes, it was necessary to assign specific 
values to the quantities involved in the time parameter 
\given by Eq. (10a). The major uncertainty in evalu 
iting \ lies in the specification of the heat-transfer co 
eicient /,. For the illustrative purposes of the present 
report, the simple empirical formula for the case of a 
turbulent boundary layer (see, for example, reference 4 

h = [(0.030)k,/x]R,*/?Np,'? (14 
was used as a basis for specifying #. Since the present 
calculation requires a value of h constant along the 
chord, rather than the variable 4 given by Eq. (14), the 
value given by Eq. (14) at the midchord, « = c 2, was 
hosen. Thus, 


(lo 


(0.030)R,(M,a)*?Np,' 3 
1, = a 


1/5;_.\4/5 
(¢/Z (y) 


The time + in terms of the time parameter \ then be 


mes, from Eqs. (10a) and (15), 


T= l(d Cc’ ”( Pinly )/BM + aN (16 
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where B 0.0689)k, (a/v) °Np,* 3 17 


For the assumed altitude of 50,000 ft., the following 


values were assigned to the parameters in Eq. (17 
a 971 ft./sec. 
vp S.16 0 sq.ft sec 
Ry 3.16 KX 10°° B.t.u., ft.-sec.-°F 


| Pe 0.9 
The value of Pin l for stainless steel was taken as 59.1 
B.t.u./it.*-"F. 

It should be remarked that recent boundary-layer 
studies indicate that Eq. (14) for h may be seriously in 
error at supersonic speeds (see, for example, reference 5 
However, as has already been mentioned, the magnitude 
of h does not influence the minimum value of G/J/,,,/G/ 
where instantaneous acceleration is assumed; only the 
time at which the minimum occurs is effected by i. The 
magnitude of /: does influence (GJ,,;, GJ for the case 
of finite acceleration, but, as will be shown in the next 
section, the magnitude of this influence may be quite 


small for moderately rapid accelerations 


Finite Acceleration 


Appendix A also contains a derivation of the tem 
peratures, stresses, and torsional stiffness changes of a 
solid double-wedge wing subjected to the following 


initial and flight conditions: 


(1) For time 7 0, the wing has zero velocity. 


2) The initial wing temperature is 7, the free 
stream temperature. 
(3) For time 0 < 7 < 7, the wing flies at censtant ali 


; 


tude with the Mach Number variation 
Ul MAr/r)" IS 


Here .\/, is the final and largest flight Mach Number 
and 7,;1is the time at which it is attained 


? 


(4) For time 7 > 7r,, the wing continues to fly at 


constant altitude at \lach Number J\/ 


The simplifying assumptions outlined in the last sec 
tion continue to apply; but one new assumption is 
needed, namely the variation with time of the heat 
transfer coefficient during the accelerated portion of the 
flight. 
variation; thus, while the Mach Number is changing 


Eq. (15) was used as a basis for this assumed 


h = hf M/M,)*”* 19) 
where /i, is the value of /: corresponding to \/,. From 
Eqs. (18) and (19), 

h hr/rs)° 20) 


during the accelerated portion of the flight. 


The final results for the effective torsional stiffness are 


, (To? — T.W.)/(d/c)*} P(A, dy 2] 


where 
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Here, as in the previous section, A is the time parameter 


lefined by Eq. (10a 1.€., 

r 2h pT) PmCmA 10a 
ind \,, the value of the time parameter at which Mach 
Number .|/, 1s attained, is 

Ny Dh i ¢/ Dotull 24 


It should be true that as \,— O, the results for effective 
stiffness reduce to those for instantaneous acceleration, 
with .\/ O (or 7 / 


follows from the fact that, as can be readily verified, the 


That this is indeed so 


first derivative of /T is 


H"(2 IPOA 
Consequently, for A nN 
n P(X, A 2/5)H"(r F1"(\) | 72 A(A 
Thus, for Ay > O 
( Kka’G ; / l ad tA \ 
0) 
which is precisely the result given by Eq. (8) with 


rhe effect of various rates of finite acceleration to a 
ven Mach Number J/, on the minimum value of 

G/ that occurs can now be determined. The 
larger \,, the longer is the time to reach .1/,, and hence 
the slower is the acceleration. It is found that, for any 
the largest value of P(A, A 


alue of \ greater than \,,; in other words, the minimum 


always occurs for a 
fective torsional stiffness always occurs after the 
steady portion of the flight begins. Fig. 10(a) shows the 
ation with A; of the values of \ which yield the 
maximum value of P(A, A;), as given by Eq. (22b). Itis 
seen that as A, 


ceur closer and closer to A,. Fig. 10(b) shows how the 


increases, the maximizing values of A 


rot 
10 


m(X [P(A, Xp) mar [POA, 0) Inn 27 
aries With A ,; here, of course, [| P(A, 0) |,,.2 1s the same as 
A(X 0.0456. For a given wing, altitude, and 


maximum Mach Number .\/,, the ordinate in Fig. 10(b 
sequal to the ratio of the maximum torsional stiffness 
hange due to a finite acceleration to the change due to 
acceleration. Consequently, if A, is 
10(b 
nite acceleration on the results previously found for 


istantaneous 


known, Fig can be used to assess the effect of 


me 


ODT} MIC HEATIN( s 
the minimum effective stiffness for the case of instan 
taneous acceleration. Thus, from Eqs. (12) and (13) for 
instantaneous acceleration 
( GJ n1lA 0.0456 Ma/u x 

i {Z Zo 
Gy GJ 
m(XA7) {0.00821 |[Aa GI[7. MV d 2Sb 


The curves of Fig. 5 for the illustrative example were 
obtained from Eq. (21), and showed that, for this ex 
ample, moderately rapid accelerations produced maxi 
mum stiffness losses only slightly less than those in 
curred by instantaneous acceleration. Phe general ex 


pression Eq. (2Sa) and Eq. (2Sb), and Fig. 10(b), indi 
cate that, in general, the relieving effects of finite ac 
is about one or less 


celeration are quite small if A 


With the use of Eq. (15) for /y, it is found that in terms 
of the average acceleration .\/, 7,;in Mach Numbers per 


unit time 


and 1s given by Eq. (17 


where B depends on altitude, | 
7,in Mach Num- 


Fig. 11 shows curves of A, against .\/ 

bers per minute, for an altitude of 50,000 ft., d « 0.03 
and \/ 3. The material is stainless steel, with 
Pmt o9.1 B.t.u. ft F.; separate curves are shown 


for chord lengths of 3, 6, and 12 ft rhe corresponding 
values of A, for other thickness ratios, Mlach Numbers 
chord length, and materials may be readily found by 
scaling according to Eq. (29 rhe value of (p,,< for 
both aluminum and titanium 1s about 36.6 B.t.u. ft 


F. An auxiliary abscissa in Fig. 11 shows average 
icceleration in g’s. 
Inspection of Fig reveals that for a very wid 


are quite small, s« 


range of parameters the values of 
that the 


corresponding values of m(A are close t 


unity. Consequently, the results for infinite accelera 


tion have practical validity in many cases of finite 
acceleration. Of course, when accelerations are so slow 
as to modify the instantaneous acceleration results ap 
preciably, the stiffness effects can be estimated quant 
tatively by the use of Fig. 10(b) for m(A 

Another important conclusion that may be drawn 
from Figs. 10(b) and 11 is that inaccuracies in the est1 
mate of the heat-transfer coefficient are not likely to be 
serious insofar as the effect on torsional stiffness is con 


cerned. Since, as Eq. (24) shows, A, is proportional to 
h,, a change in h, by even a factor of two or three is un 
important if A, is well within the range where the in 
stantaneous acceleration results are nearly valid 

These considerations concerning the magnitude of the 
assumed value of /i, are without reference to the errors 
introduced by the assumed constancy of / along the 
chord. The next section attempts to shed light on the 
nature of the modifications of the results that are intro 
duced by permitting / to vary in the chordwise direc 


tion 
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Effect of Heat-Transfer Coefficient Variable Along the 

Chord 

The constant-thickness wing of rectangular cross 
section shown in Fig. 2(b) would not suffer any loss of 
torsional stiffness if, as was assumed for the double- 
wedge, the heat-transfer coefficient were constant along 
the chord; the wing temperature would simply rise 
uniformly and no thermal stresses would develop. 
Consequently, it may be instructive to consider the 
effect on the constant-thickness wing of a variable 
heat-transfer coefficient; in this way, an idea can be 
gained of the errors involved in the preceding calcula 
tions for the double-wedge wing which ignored the 
chordwise variation of /. 

From the approximate Eq. (14) for the heat-transfer 
coefficient of a turbulent boundary layer, it follows that 
h varies inversely as the fifth root of the distance from 
the leading edge; thus, 


his) = héc)/ (a/c) 30 


Consequently, in accelerated flight, the wing tempera 
ture in a constant-thickness wing would be largest at the 
leading edge and decrease continuously toward the 
trailing edge. In Appendix A, the temperatures, 
stresses, and stiffness changes are calculated for an in- 
stantaneously accelerated constant-thickness wing; 
except for the use of Eq. (30) for /, all of the conditions 
and assumptions used for the double-wedge wing were 
applied. The final result for the minimum effective 


stiffness 1s 


(GJ .7;/GJ); 

1 — (0.00606)(Ea G)) [Tin — Ty? | (dc) Bla 
or 

GI. 4 GJ), 


1 — (0.00109)(La G)[7T (MW / Mo*)/(d/c)?] (31b) 


In Fig. 12, the results of Eq. (31) are shown for con- 
stant-thickness stainless steel wings, with .l/) 0 and 
My, = 3, and are compared with those found for the 
double-wedge wing on the basis of a constant /. It is 
apparent that the stiffness loss of the constant-thickness 
wing is quite substantially less than that of the double- 
wedge wing. It is reasonable, therefore, to expect that 
the results found for the double-wedge wing with con- 
stant / would be modified only slightly if the variation of 
h along the chord were taken into account. This ex- 
pectation is all the more plausible when it is realized 
that the peak influence of a variable /; need not coincide 
in time with the occurrence of the maximum stiffness 
loss found on the basis of a constant /. 

As a matter of auxiliary interest, the results found for 
the solid constant-thickness wing are applicable with 
only a slight modification to a hollow wing of constant 
skin thickness, with no internal structure. For the 
hollow double-wedge wing shown in Fig. 2(c 
(GF te GF aun = 1 (O.00S08)(/2a G) X 

ty ae Tae |/(d/c)?}{ = (32a) 


or 
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and, again for stainless steel wings in the stratospher 
with MM, = 0, and M, 3, the results for the holloy 
wing are shown in Fig. 12. For the hollow wing of co; 
stant-skin thickness, as for the solid constant-thickness 
wing, no stiffness loss would be found if the heat-trans 
fer coefficient were made constant along the chord. | 
should be remarked, however, that these results for the 
pure hollow wing may be somewhat academic, since thy 
addition of internal structures such as shear web 
would change the results appreciably. Depending or 
the precise distribution of web material, the multiwel 
wing would probably suffer stiffness losses substantially 
greater than those shown for the pure hollow wing 


Free-End Effect 


The assumption that plane sections remain plan 
used in the calculation of thermal stresses, is stricth 
valid only far from a free end. Ata free end, of course 
no axial stresses can exist, and the effect of the free end 
is felt for some distance inboard from the tip. Thus, th 
effective torsional stiffness remains at its original valu 
at the tip, and is actually variable along the span, ap 
proaching the values calculated in the preceding sections 
with increasing distance from the free end 

In Appendix B an approximate calculation is made of 
the spanwise variation of the effective torsional stiffness 
of instantaneously accelerated, double-wedge wings of 
constant chord and finite span (see Fig. 15). In addi 
tion, the change in the first natural torsional frequency) 
of the wings considered as free-free beams is calculated 
approximately; the change in frequency or, more pr 
cisely, the change in the frequency squared —consti 
tutes a measure of the integrated effect of the torsional 
stiffness changes along the span. The results of th 
calculation are presented in Fig. 14, which shows, as a 
function of the aspect ratio b,c, the ratio of the maxt 
mum change in frequency squared to that calculated on 
the basis of the constant effective stiffness given by the 


elementary theory [see Eq. (S)]. It is evident from this 


stiffness changes are concerned little error is incurred 1n 
the use of the elementary thermal stress theory for 
reasonable aspect ratios; even for a value of wing aspect 
ratio as low as 1.5, the error is only slightly greater that 


10 per cent. 


Effect of Initial Twist 

All of the foregoing results for effective torsional stiff 
ness are based on the assumption that the wing 1s 
initially ‘“‘perfect’’—that is, untwisted-—and hence re 
mains untwisted while it experiences thermal stress. A 
wing with some initial twist would undergo additional 
twist as a result of the thermal stress distributions con 
sidered; indeed, as the condition for thermal buckling 
of the perfect wing is approached, this additional twist 
would become substantial (though not infinite, becaust 


of nonlinear large deflection effects). It follows from 
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13. Wing geometry and coordinate system for finite-span 
wing considered in Appendix B 


the theory of twisted beams'® that for a sufficiently large 
inal twist, the effective torsional stiffness would be con 
siderably higher than that of a perfect wing subjected 
to the same thermal stress.* On the other hand, the 
magnitude of such a final twist would probably be un- 
desirable from both structural and aerodynamic con 
siderations. In summary, then, the behavior of a pre 
twisted wing would be qualitatively as follows: for 
moderate initial twist, the effective torsional stiffness of 
ithermally stressed wing should be approximately the 
same as that calculated for a perfect wing as long as the 
condition for thermal buckling of the perfect wing is not 
ipproached; when the thermal buckling condition 7s 
pproached, the pretwisted wing would probably suffer 
excessive deformation, even though its final torsional 
stiffness would not be as low as that predicted for the 
perfect wing. 

These considerations are consistent with the tests 
ind calculations of Vosteen and Fuller’ who dealt with 
i flat cantilever plate subjected to thermal stress dis 
tnbutions similar to those considered herein. The 
thickness ratio was quite low (1/SO), aggravating the 
elects of initial displacements, as evidenced by the 
lact that appreciable distortion occurred during heat 
ing, with a consequent large deviation of measured 


torsional frequencies from those predicted by an 


"This observation was originally made to the authors by 


Joseph N. Kotanchik of the NACA 
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analysis in terms of flat plate theory analogous to the 
beam theory of the present paper. 


CONCLUDING REMARKS 


The results exhibited show that serious losses in effec 
tive torsional stiffness can be produced in thin solid 
wings by the thermal stresses that arise during and afte 
accelerated flight to supersonic speeds The thermal 
stresses in solid wings are mainly attributable to chord 
wise mass variations; consequently, solid wings having 
thickness distributions more uniform than the double 
wedge profiles considered in this paper may be expected 
to suffer smaller losses of stiffness. Similar losses of 
stiffness can occur in built-up wings, wherein the 
origin of thermal stress resides primarily in temperature 
differences between internal and external elements 
Stiffness losses for built-up wings would probably be less 
serious than for solid wings; but, in any case, the mag 
iitudes of the torsional stiffness changes of any shape ol 
wing cross section can be estimated by means of the 
simple formula 


J ord 


once the thermal stresses o, on the cross section are 
known. 


APPENDIX A \NALYTICAL CALCULATIONS 


Let it be assumed that, at station x along the chord of 
a thin solid wing of arbitrary cross section, an element 
fdx) retains the heat transmitted to it from the bound 
ary layer at the upper and lower surfaces. Then, with 
the additional assumption that the temperature of the 
element is uniform through its thickness, the heat bal 


ance equation governing the temperature / 1s 
2h| 7 — 7] | 


j 


PmnCnt(Ol Or \ l 


where, in general, the heat-transfer coetlicient / is 
variable in both space and time. The adiabatic wall 


temperature 7), 1s taken as 





| Tu) 1 + [n-(9 1) /2)M*} \-2 
O; ™ 
} - 
8} 
e 6} 
‘Max 
= = | ¢ 
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2) S 0 5 2 eo 30 
D/C 
Fic. 14. Effect of wing aspect ratio on the ratio of the actual 
change in first torsional frequency squared to that based on 
elementary theory (effect of free end neglected) for solid double- 
wedge wings considered as free-free beams E/G 2.54 
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where .J/ and 7°., are the free-stream Mach Number and 
temperature, respectively. In the ensuing calculations, 
/’,, is assumed constant throughout the various flight 


histories to be considered. 
Instantaneous Acceleration 
Let the Mach Number history of the wing be 
M/ Mo tor rT <Q 
M, for 72°90 


0 be 7,,, 
the adiabatic wall temperature corresponding to .\/ 
M>. Then, for 7 > 0, Eq. (A-1) 


and let the wing temperature at time 7 


nay be written as 


(Ce 1 (1 Fae”) 


(Pm€nl/2h)(t/d)[o(7 Taw”) /Or] (A-3 
where d is the maximum airfoil thickness and 77, 
corresponds to JJ = J/,;. Write 
6= [T — Tan ]/ [Tox Ta | 
rn (2 Py d) hdr 
Then Eq. (A-5) becomes 
l 4 (t/d)(O0/OX) (A-4) 


The solution of Eq. (A-4), with the initial condition 


a OatA 0, is 


a l e : F A-5 


This result for the temperature history is for thin solid 
wings of arbitrary cross section and with arbitrary 
variation of / along the chord. Thermal stresses and 


torsional stiffness changes will be calculated for two 


special cases. 


Solid Double-Wedge Section; h Constant Along the Chord 


For the symmetrical double-wedge airfoil shown in 


Fig. 2(a) 
tid 2x /¢ ie 2/2 
2 {1 (X/ ¢ c= 6/2 
Hence, 
6 e HF= 8 = 1 \-6 
| ./211—(B/2 1<g <2) (AA-0 
where @ 2x/c. If hk is assumed uniform along the 


chord, then A is independent of x (or 6). For the sym 
metrical temperature distribution given by Eq. (A-6), 
the assumption that plane sections remain plane re 
quires that the longitudinal strain distribution €, be con 


stant over the cross section. Then, by Hooke’s law, 


€, = constant (o,/E) + all — Ty,© (A-7 
where any effect on ¢, of stresses o, and a, in the plane 
of the cross section is neglected. Application of the 


equilibrium condition of zero thrust 


to,dx 0 
0 


LERONAUTICAL SCIENCES DECE 


or Oya 0) 


permits the evaluation of the constant in EK ‘7 
with the following result for the stress 


0 loa | ( / / / / 
or 
a] Ka(l / 6 a ew 
where 
>] 2] >] 
Gow | 30(3)d3 bd 3 2 B8( 3)d3 
JZ WV /J WV / 0 


Evaluation of Eq. (A-S) with the use of Eq. (A-6) for 


provides the results for the stresses viven by Eqs ia 


and (7b) —1.e., 


0; 
: poe Nye 
EalT, ha 
A fa \ Q) CiZ 
e 1 — Aj 
The effective torsional stiffness given by Eq. (4) is 


evaluated by approximating the radial distance r fron 


the midchord by | (c 2) v ; then, for the symmetrical 


stress distribution under consideration, 


“ES G/ T 2 \¢ » vi-fo.d\ 
. if] 
or, with 
f 4 Ov da 
GJ, G./ cid /4 3(1 3 3 


or 


GJ «+1 GJ +L (3 Ge d 3( | O)-0,a0 A-9 


Evaluation of the integral in Eq. (A-9), with o, given by 


Eq. (7), provides the result for G/ GJ given by E 


io) te 
GJ .77/GJ (Ea/G) X 
(2a / | (dc)*h AQ) (8 
where 
A{\) = (1/8)le “(4A + 7A dh?) + 


Hi(—X)(10A* + SA® + A J 


Constant-Thickness Section; h Variable Along the Chord 


With the thickness ¢ constant (= d) and the variation 
along the chord of the heat-transfer coefficient governed 


by the equation 


whet 


lhet 


wher 
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h(x) = hi(c)/(x/« 30 
he te! iture history at each chordwise station x is 
described by 

id | ( 


where the time parameter A is now defined as 


> 


[he temperature is not symmetrical about the mid 
chord and the assumption that plane sections remain 
plane permits the possibility of rotation as well as 


translation ol the cross section. Hence, 


where dp and ad, are constants Setting 
é dy + ayn (o,/f£) + all l 


nd imposing the condition of zero thrust 


. 
a,x 8] 


J 0 


iid of zero moment 


va,dn 0 
a 


f the constants a 


permits the evaluation « and ay 


With the introduction of the notation &€ vc, the 


stresses then turn out to be given by 


i 
Eal/ / | | Adé 6 cAdé 
fo) 


>] . ] 
O£& Adé » cadé ( 
Le/’O0 JW a 
here 0 | ( 
Phen 
ka(l / LA (Xr + 6B \ T 
GEA (A 2B(A)|] + e (A-10) 
were AQ) = fe OM 
) 


>i 
Bin ger ae 
i] 
rhe quantities A(A) and B(A) can be evaluated by 
means of successive integration by parts to give 


+ 


1) /k1(—A 


A(\ e “/24)(24 6A + 2d" A? + AS] 
(A°/2 
B(\ N '30240)A (A) + (e>*s/6048) X 
3024 336A + 42d 
With GJ =(G/3)ct*, and r°dA taken as t|x — (¢ 2) }? X 
the expression for the effective stiffness [Eq. (4)] be 


mes 


2 a dé G(d/c)? A-1] 
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Substitution of Eq. (A-10) for ¢,, and evaluation of the 


integral vields 


GJ f GJ l lea G) X 
/ 1 d/c)*|Q(A 4-12 
where Q(A 1/2)A(A BX IC(A 
C(A | E dé 
\ die 
B(A) + =~ [24024 I716\ 4 
SO0860 (2002 
142A LIA \4] 


Phe maximum value of Q(A) occurs at A 0.60 and is 
equal to 0.00606, from which Eqs. (31) for (GJ, ¢,/GJ 
follow. 

The result quoted in the body for (G/,,, GJ of a 
hollow double-wedge wing of depth d and constant skin 
thickness follows from the observation that Eq \-10 


remains valid for the hollow wing if \ is redefined as 


where /, is the skin thickness. The torsional stiffness of 
the hollow wing is approximately 


y 


GJ G2 )t.cd 


so that, with r°dA Pt,|a c 2)\/?dx in the formula 
for effective stiffness 


GJ .¢7/GJ | 


{ . 
Comparison with Eq. (A-11) shows that the maximum 
relative stiffness change for the hollow wing is 4 5 times 
that of the solid wing, whence Eq. (82) of the text 


fc 1c WS. 


Finite Acceleration 

Phe solid symmetrical double-wedge wing 1s now con 
sidered to accelerate at constant altitude from Mach 
Number zero to Mach Number .\/, according to the 


relation 
VW Uri r 


For r > r,;, the Mach Number remains constant at the 


value of .J/,. 
transfer coeflicient, 


On the basis of Eq. (15) for the heat 


h hA(M/M,)*"" 19 


where /1, is the value of /) corresponding to .\/,, and, as 
before, is constant along the chord. Hence, from Eqs 


IS) and (A-14 


Che initial wing temperature is assumed equal to the 
free-stream temperature / 
For the double-wedge wing, with 0 < x < c 2, the 


heat balance Eq. (A-1) may be written 
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O(7/ / Or QN/ PCr) /(2x/c)| X 

((7 / (7 i ieee \-14) 
But, from Eq. (A-2), 

1 / CT . / U/M 
Hence, by Eq. (18), 
/ / ‘ae Mala) Os7r87r 
/ | Ter 

Now, let 

r Qh t/ Pm Cal 

r 2h pt 7/ PmCma 

s AA 
4 (7 Pot! tl - / 

Then, for O S x S c¢/2, the formulation in nondimen 


sional terms of the differential equation governing the 
temperature history becomes 

00/Os = [d,/(2x/c)][s°? — 8°70] QO<s<1) 

O# Os : 


The solution of Eq. (A-15) satisfying the initial con 


dition 6(0 0 is 
53 . 
A , (1 ( z USss il 
>Xr 
D6 . , (A-16 
[, 2/5)}] 
>A 
2 1 
where 3 2x/C. he values of 6 forx = c 2, by svi 
metry, are given by Eq. (A-16) with 3 replaced by 
2[1 (8/2) ]. 
From Eq. (A-S), with 7 replaced by 7°., the 


thermal stresses are given by 
a, Kall, / 6 6 (A-17 
where, as before, 


4, 2 


80(3)d3 


Then, from Eq. (A-9), 
Sha \ / l P / 
GJ G (d/c) 


*! 
} B(1 — B)*[Oare — 0(8)]d8 (A-18 
ae ( 


) 
The use of Eq. (A-16) for @ in the evaluation of Eq. 
(A-1IS) provides the results for G/J,,;, GJ given by Eqs. 


(21), (22a), and (22b). 


FREE-END EFFECT ON THERMAL STRESS 


AND TORSIONAL STIFFNESS 


APPENDIX B 


The constant-chord double-wedge wing of finite span 
b shown in Fig. 13 is considered to be instantaneously 
accelerated from J = JJ, to AZ A\7,, and to undergo, 


at each spanwise station, the temperature history given 
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by Eq. (6). An approximate calculation will be mad 
the accompanying thermal stress distribution in th 


wing, with cognizance of the fact that at the wing tip 
+ (b/2 | 


, the axial stress o, must vanish 
For the sake of simplicity, the analysis is based on t] 
type ol approximate elasticity theor 


Phat is, t} 


‘sheet-stringer’’ 
common in aircraft structural analysis. 
effects of normal chordwise stresses on distortion 

entirely ignored; only the axial stress o, and the shea 
stress rT both assumed constant through the thick 
ness, are considered pertinent, and obey the equati 


of equilibrium 


t(Oa, OY (O Ox) (tr 0 B 
The strain-temperature-stress relations are assumed t 
be 

€ l I: 0 ( [ [ B-9 

1 l'G)r K 

The associated complementary energy density is 
LV, G 2/6) + all / lo, + T 2G B-4 
since é Ol’. O« 


A solution for the stresses o, and 7,, may be found by 


mininiuzing the total complementary energy of the wing 
with the only stress states admitted to competition it 
the minimization process those that satisfy the equ 
librium Eq. (B-1) and the boundary conditions 


a,|x, #+(b/2)] 0 B-5 


By symmetry, attention may be restricted to o 


quadrant of the wing, say 0 S v Sc 2 and 0 


b/2. 

variables 8 = 2v c, n = 2y_ b, the equilibrium Eq. (B 

mav be written, forO < x Sc 2, with? Gd, as 
(c'b)(Oa, On) X (18) [0(Br Oo | 0 B- 


An approximate solution will be found by 


that 
B-7 


o Oo Y UI 
where (¢,),, is the elementary stress distribution give 
by Eq. (7), and g() is a function that describes how th 
magnitude of this distribution varies along the spa 


From Eq. (B-5), 


g () B-S 
From Eq. (B-7), then, 
0, . 
: p(B) g(y B- 
Ea|l 1 a 
where 
pla ee nw ff Nye * + ME —A B-10 


From the equilibrium Eq. (B-6), 


With the introduction of the nondimensional 


assulimly 
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{7 la 
( . * . 
~ |} g (n Bp(B)dp + e(n B-1] 
But 7 Q at the midchord, 8 1; also, it must be 
true that 


B8p(8)dB 0 


by virtue of the fact that (¢,),; produces no net thrust. 


Hence, ¢(n 0, and, carrying out the integration in 
Eq. (B-11) yields 
T ( e | | 
» ) ) 
( )f Uf] > “ 
Ka / / | b - 
where 
aie ER d De Bo /B)1/B} 
3} [1 AJe~* — AME —A); B-13 


fhe complementary energy in one quadrant of the wing 


1 | 


Is 


] . 


bed (°° oe 
BIVdgsdn 
} 0 0 


(WW .dxdy 
With the use of the temperature relationship 


/ l , 


and Eq 6 
for the stresses, the energy may be written in nondimen 


, as well as the expressions (B-9) and (B-12 


sional form as 


G ,) U)E (MT ¢ bad odn 
1 \O 


1 4)(2/G)(e b)200(\) [2'() [2h dn 


whert 
\ ? \ le 2D od De 2r 
1/2 Ke" \2Fi(—dX B-14 
*1 
\ ; ) Ne ; N*E 1 Np } 
3( Ne + OX 1 \ de B-15 


Numerical integration appears necessary in the evalua 
tion of co( 
The unknown function g(7z) is now determined from 
the condition that it minimize F. Application of the 
calculus of variations yields the Euler equation 

k?g k B-16 


S NS 


where 
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Eq. (B-18) merely states that 7,, = 0 at n 0 as it 
should by symmetry. The solution of (B-16) that satis 


fies the boundary conditions (B-S) and (B-1S) is 


(7 l cosh kn cosh k B-19) 


rhe factor g(7), which modifies the stresses o, according 
to Eq. (B-9), 
factor on the reduction of torsional stiffness 


may also be used directly as a modifying 
Conse 
quently, applying the factor g(m) to the second term in 


Eq. (8) gives 


> 17 / d-c)*| A(A\)g(n, A B-20 


where g has been shown as a function of A as well as 7 
to emphasize its dependence on time through the 
quantity k defined by Eq. (B-17). From Eq. (B-20) it 
follows that the reduction in stiffness is zero at the wing 
tip, where g = 0; at the midspan, the end effect changes 
the stiffness reduction from that previously calculated 
without end effect by the factor | | cosh k), and 
for k large, this factor is very nearly unity. For the 
limiting case of infinite aspect ratio, Eq. (B-20) reduces 
to the result found by the elementary theory. A meas 
ure of the overall importance of the changes in stiffness 
will be found by calculating approximately the change 
in the square of the first torsional frequency of the wing 
considered as a free-free beam. On the basis of the as 


sumed vibration mode shap« 
¢ gi sin (ary b 


the frequency w may be estimated by a simple Rayleigh 


solution that equates the maximum potential energy 


The result for w” is 


r*|GJ 1,b7| X 

l )Eal7 / d/‘c)*| A(X) P(A B-21 
wher¢ 

| tanh k B-22 


For the unheated wing, the frequeney squared is 


7 


so that Aw", the change il 


Aa Ww lea l i d ( x 


When the elementary stress distribution is considered 
to act all along the span, the result for Aw? is given by 


Continued on page 1108 











The Tip C ion for Wind-T | 
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1e Lip Correction to1 ind- | unne 
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ests of Propellers 
\\ 
THEODORE R. GOODMAN* 
Cornell deronautical Laboratory, Ine. | 
INTRODUCTION blade element theory and assumes a lightly loaded 
; ; ; propeller with a large number of blades. The analysis 
FPVNE DIAMETER of propellers of current interest ss é ee oe 
. aia appears in Chapter VII of Glauert’s article in Durand [ 
nay be as large as 19 ft., and for VTOL applica : : ; 
: : ? - ’ and is reproduced here almost verbatim from. that ce 
tions even larger. It 1s clear, therefore, that with 
ig ee , : P , source. k 
present-day facilities it 1s impossible to test such 
propellers in a wind tunnel whose cross-sectional area Define ‘a 
is large compared with the cross-sectional area of the B number of blades 
. rin . . . ] 
propeller disc. The tunnel walls will be near the A local circulation : 
propeller tips, and, as a consequence, the local flow ° ne rads 
° ° y ° 6 ) air density 
near the tips will be altered. Whereas this effect 
i : 2 Q angular velocity of propeller Ww 
was present for smaller propellers as well, it was of “5 angular velocity of the air behind propelle: 
7a . » ; ni hh 111 
negligible importance, and the bulk of the wind r radial distance along black 
tunnel wall interference could be computed as a simple / = thrust 
blockage correction. With a “large’’ model in a —v torque 
a af ° | forward speed 
small” tunnel, the altered flow at the tips can cause ' 
=P Y : : ; ul axial velocity at propeller dis¢ re 
an additional correction of the same order of magnitude Or/V ; 
: ge C 
as the blockage correction. ' 
iia : . ‘ : ? ‘he ex SSI for » eleme Ss d al 
rhe primary reason for testing full-scale models is Phe agin lor the elements of thrust and 
the difficulty of designing models with dynamic simi torque on a blade element are obtained from cascade 
. . , ‘ . -orv (see reference “hi rr V Si 
larity at high speeds. Nevertheless, as a first step theory (see reference 1, Chapter \ 
‘ : : : ‘ le 
in computing the tunnel-tip correction, a low-speed dT /dr = BKp[2 (w (2) ]Q tl 
é ; epee ; " 
(incompressible) analysis will be performed. Further i0/d BI “ 
‘ 7 . ‘ . ( r 3 ) y | 
more, it will be assumed that the propeller is lightly an —_ . 
° ° - 7p Ci 
loaded and consists of a large number of blades. The also BK 2trw! 3 f 
. . . ° ec 
per cent error introduced by these approximations is sa ; 
a aS a: . . . . and an alternative expression for the thrust is obtained tu 
certainly less important in computing interference ; ; ; 
Riga re bsietta from a consideration of the axial momentum as CO 
corrections (which themselves are presumed small) 
than they would be in computing the aerodynamic dT /dr barr pu(u |” { 
characteristics of the propeller itself. Suppose, for ' ' ; 
. . : . . sy equating the two expressio1 ; > element ol 
instance, that the interference correction is 10 per cent, h ~. . sa ' | ‘ : caren ie = _ . n ; i 
. . . thrust, it appears that tne axial and roti nal velocities 
and the error in the correction is 20 per cent, then the a a : : er id 
. : are related by the equation 
error in the corrected data will be only 2 per cent. lass atl metas sh 
The analysis which follows is accompanied by graphs uu |’ (1 2)fQ w/2) lw? re 
for the computation of corrections in thrust and " 
— ae : : and, finally, the loss of energy due to the reaction oi a 
efficiency for two-, three-, and four-bladed propellers ; . ae 
: . . we . . the blade elements under consideration is 
as a function of advance ratio and ratio of propeller 
radius to tunnel radius. The theory might possibly dE/dr Q(dO d VidT dr 6 
be extended to the analysis of shrouded propellers 


and ducted fans. 


THE PROPELLER OF HIGHEST EFFICIENCY 


In order to make the analysis which follows perfectly 
clear, a preliminary glimpse will be taken at the pro- 


Passing now to the approximation of a lightly loaded 


propeller, Eq. (5) may be taken to be 
V(a V (1, 2)Qa» 


or, by virtue of Eq. (3), 


peller of highest efficiency. This analysis is based on u 1” (BKQ/4arV 
ceive Yecember 28, 1955 . > - 
Receives Decewm a and Eq. (1) for the element of thrust becomes 
* Formerly, with the Aerodynamic Research Department : tel 
Now, with Allied Research Associates, Inc se sees 
The author is indebted to the students and staff of the Aero dT /dr BK p®r 8 
nautical Engineering Department of Cornell University who 
offered a suggestion which led to a simplification of the analysis Also, by virtue of Eqs. (3) and (7), Eq. (6) for the loss 
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of en 
1} BKp|(BAK@& ri" BRV 42?) |r] 
9 
BK pQ 4m) [(Qr/ I t V/ Qr \ 
Writ 
7 BKQ/2rV 10) 
[hese last two expressions become 
dl‘/dr (2rpV4/Q)yx 
dk dr 2rp Vy? /Q) (C1 + x?) 2a (11 
[The optimum condition is now derived from the 


condition that the ratio of the increments of |’7° and 
Ek, due to at any 
radius, must be independent of the radial coordinate 
Chapter VII). 


an arbitrary small increase of y¥ 


, (see reference 1, The condition is, 


therelore, 


ly(1 + x = AS 
where A is a constant, and this condition reduces 
immediately to 
BKQ/2rV2 = y = [Ax?/(1 + x?)] (12) 


Chis completes the quotation from Durand with 
regard to the efficiency. In 
Chapter VII, Section 4 of his article, Glauert presents 
-randtl 
for accounting for a finite number of blades by con- 
The author 


propeller of highest 


an approximate method originally given by 


sidering the radial flow near the blade tips. 
has reworked the Prandtl solution taking into account 
the presence of a nearby wall. The distribution of 
thrust near the tip is markedly affected if the tip is 
the correction to thrust and 


close to the wall, and 


efficiency, which it is necessary to apply to raw wind- 


tunnel data, may, under certain circumstances, be 
considerable. 
[THE PROPELLER Tip NEAR A WALL 
The model of the propeller will be taken to be 


identical to the Prandtl model. The system of vortex 


sheets is replaced by a series of parallel lines at a 


regular gap and extending indefinitely to the left, as 
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in Fig. 1. The gap represents the normal distance at 


the boundary of the slipstream between the successive 


vortex sheets, and hence 


s = (27R/B) sin @ I2erR B)(A/V1+A 13 


where ¢; 1s the angle of the screw surface at the bound 


ary of the slipstream, defined by the equation 


tan ¢@ = V/QR =r 
At a distance to the nmght of the boundary of the 
slipstream is a vertical wall infinite in extent. The 


distance d is thereby seen to represent the difference 
between the tunnel radius and the radius of the pro 
peller disc. 

The flow around the edges of the vortex sheets of 
the slipstream in the presence of a wall may then be 
estimated approximately from the flow around the 
edges of the system of semi-infinite parallel lines as 
shown in Fig. 1, when the system of blades moves 
downward at right angles to the lines with uniform 
velocity. From considerations of symmetry, it is 
seen that the lines BC and ED are equipotential lines, 
and hence the entire solution, which is periodic in 
nature, can be synthesized by considering merely the 
polygon ABCDEF in the 2(= x + By 
using the Schwarz-Christophel theorem, the inside of 


iy) plane. 


this polygon can be mapped on to the upper half of 
the ¢ plane. This plane is shown in Fig. 2, and the 


mapping function is given below. 
¢ = ()cosh[x(d s) ‘s]{ /[eosh( md /s) | 14) 


In the ¢ plane, the complex velocity admits square 
root singularities at the points B C D E, and, at the 


same time, the velocity at infinity must be finite. 
If w represents the complex potential, then the solution 


must be of the form 


dw/dé = (bF + o)/ V (&? 1) [sech*( ad /s) et 15) 


where } and c are constants. 
Since dw ds = (dw dt)(dt/dz 
be obtained 


the complex velocity 
with the aid of Eq 


the 


in the z plane 
(14 


expressed in terms of 


may 


and also by using Eq. (14 result may be 


The resulting equation is 


dw b cosh | a(d s)/s| + ¢ 
16 
dz V cosh?|r(d cosh?( ard 
As Rez —~ ~, dw/dz —~ 0. If this condition is 


applied, then the result is that 5 0, and, Eq. (16 


reduces to 


dw/dz = c/V cosh?|r(d Z)/S cosh?( ad/s l7 


In order that the solid lines as indicated in Fig. | 
be streamlines, and the dotted lines be equipotential 


The 


lines, the constant c must be a pure real number 
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circulation fall-off near the tip is obtained by integrat 
ing the tangential component of velocity, u, from 
P to E, and then from E to P’. 
set 3 x and integrate Eq. (17) from —a to 0. The 
factor 2 takes into account both the upper and lower 


surfaces. 


" dx 
I 2 IS 


“-4 V cosh? [r(d x)/s cosh?(ad ss 


This integral can be evaluated by the substitution 


cosh | a(d x)/s| = 0. The result is 


F (2sc/m) cosh(md ‘s)cd— X 

({cosh(ad ‘s)|/ jcosh [m(d + a) s]}) (19 
where Glaisher’s notation for the Jacobian elliptic 
function has been used, and the modulus is sech (ad/s). 


When a becomes infinite there is no fall-off, and F 


must be unity. This condition determines c. Using 
the identity cd~'g K sn 
plete elliptic integral of the first kind, and also using 


the fact that sv ~!(0) 


1g, where A is the com 


0, there is obtained for ¢ 
( wt 2 sK cosh (ads 20 
whence 


F (1 ,cosh|[ar(d + a)‘s]y) 


K\cd~'({eosh(ad s 


In this formula a denotes the distance of the section 
under consideration from the edge of the lines, and must 


be replaced by (R — r) in applying the result to the 


system of vortex sheets of the slipstream. Write, 
therefore, in accordance with Eq. (13), 

7 Tas [B(R r)V 1 + A7I/(2RA (22 

g rd (BdV 1 + XA?) /(2RXr (23 

F 1 AK cd~' j (cosh g)/[eosh(f + g) |; (24 


The quantity F represents a reduction factor which 
must be applied to the circulation about the blade 
element at any station due to the proximity of that 
station to the tip. If this correction is applied, then 
by virtue of Eq. (12) the optimum distribution of 
circulation along the blades becomes 


y¥ = AXF/(1 + xX (29 
where F is given by Eq. (24) and f may be expressed as 


f = [BU — AwX)V 1 4+ A?d/20 (26) 


» 


The relationship between F and f is given in Fig. 3 


for various values of g. Since the quantity g is directly 


LON AYU TI< 


For this purpose 
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tions shown in Fig. 3 indicate that the presence of th 
wall tends to make the tips less effective as tips, i, 
the circulation at any station is larger when there ; 
a wall present than when there is not. This behavyio; 
might have been anticipated; because of the constrain; 
of the wall, less air is able to flow around the propeller 
disc and more air must flow through it. The gain } 
thrust caused by the presence of the wall is purchas, 
at a price, however; since the distribution of thrust 
is such that the propeller operates at highest eflicienc 


Ils 


wind-tunnel w; 


in free air, the presence of the 
necessitates a loss of propulsive efficiency. 

It will be assumed that the thrust distribution o; 
all propellers can be represented by the ideal optimum 
distribution as given in the preceding section modifie 
in accordance with the tip correction as computed i1 
this section. This assumption can be justified on tw 


grounds. First, any particular propeller will b 


designed so that the thrust distribution will tend to b 


optimum as a principle of good design. Hence, th 


optimum distribution cannot be extraordinarily far 


from the actual distribution for any particular pri 


peller. Second, it is presumed that the compute 


correction to measured data as obtained in the wind 
tunnel is small in order that conditions in the wind 


tunnel approximate free air conditions. If this is so 


a rough approximation to the wind-tunnel wall corre 


tions is sufficient for practical purposes; any small 


errors in the corrections are of higher order in the data 


THE CORRECTION TO THE THRUST 


The ultimate purpose of the present study is t 
determine quantitatively the difference in behavior 


propeller 


of a propeller in a wind tunnel from that of : 
in free air. The quantity F for the propeller in fre 
air can be obtained by allowing d(or g) to approacl 
infinity. 
with the result given by Glauert 
; then, by virtue of Eq. (24 


Denote the value « 
F for g—> © as F 
F (2/7) cos le 2 

If reference is made to the first of Eqs. 
with the aid of Eq. (25), the thrust of the propelle 


in the wind tunnel is given by 


] » 
an DS 


[ a 2rpl ‘Ay Ff 


Q7(1 + x 


The equation which corresponds to the propeller 
in free air is identical to this, except that F is replace 
by F The difference between the thrust in the 
tunnel ‘and the thrust in free air divided by the thrust 


The result of this limiting process coincides 





proportional to d, it is seen that g = © corresponds in free air gives the per cent increase in thrust due t 
to the case of no wall at all, and the smaller the value the presence of the walls. Let / Ax, and then, 1 
of g, the nearer the wall is to the tips. The distribu use is made of the identity cd~'g K — sn~'q the 
result is 
f tedt ). (ll T 7 cosh g | 
ta sin”! e sn 
67 Jo y° + 2K cosh(u ui +g) § 29 
1 ti dt — 
( 1 e “ 
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where 7 BV 1+ X27/2dr 30 


These integrals have been evaluated by Simpson’s 
rule using the Smithsonian Elliptic Functions Tables. 
In Figs. 4 and 5 is shown a plot of 67° 7 vs. X for a 
2-, 3-, and 4-bladed propeller, in one case when the 
ratio of disc radius to tunnel radius is 0.87, and in the 
other when this ratio is 0.93. It is readily ascertained 
that the per cent increase in thrust due to the proximity 
of the tips to the wall may be larger than tolerable 


experimental error. 
THE CORRECTION TO THE PROPULSIVE EFFICIENCY 


The propulsive efficiency is defined to be 


n = VI/Q2Q 31) 


E = 090 —-— VT (32) 


In view of Eqs. (11), it is convenient to express the 


efficiency in terms of the thrust and energy 


7 = VI/(E+ VT) 3: 

A small change in 7 or F will cause a corresponding 
change inn. Denote this change by 6. The operator 6 
will be interpreted to mean the small difference between 
a quantity in the tunnel and in free air. By logarithmic 
differentiation it is found that 
bn/n = (67/T) — [(6E + V6T)/(E + VT)) 

= (67/T) — (6E/E)(E/(E + VT)] —- 

(67/7T)(VT/(E + VT)} 
= (67/T) — (1 — n)(6E/E) — (67 /T) 34 


or dn/([n(1 — n)] = (67/T) — (6E/E) (30 
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The quantity 67° 7 has already been determined j 
Eq. (29). 


with Eq. (25), 


By using the second of Eqs. (11) « upled 
the energy may be found. Once again 


the difference is taken between the value in the tunnel 


and that in free air, and divided by the value in fre, 


air. The result is 

>! fd -_ fe >] dt F 

6k: E = f fk f \ 36 
J0 A- + [- J 0 \* +— } 

where F and F.. are given by Eqs. (24) and (27 


respectively, and if f is expressed in terms of ¢ the 
result is f = Bil — t) V1 + A?/2A. Once again the 
integrals have been performed by Simpson's rule, 
The results have been applied to Eq. (35) and graphs 
for the efficiency correction have been plotted in 


Figs. 6 and 7. 
CONCLUSIONS 


When a full-scale propeller is tested in a wind tunnel 
whose cross section is about the same size as the 
propeller disc, it is necessary, in addition to the blockage 
The effect 
of the walls is (1) to increase the thrust, and (2) to 


correction, to apply a wall-tip correction. 


decrease the propulsive efficiency. 
The two corrections have been obtained quantita 
tively in an approximate manner for incompressible 


flow. 
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Normal Forces on Flat 


Plates in Low-Density 


Supersonic Flow! 


D. M. TELLEP* ann L. TALBOT** 


University of 


SUMMAR\ 


Normal forces on rectangular flat plates in supersonic, rarefied 
measured over the range 200 Re < 2,400 and 3.7 


$1 at angles of attack up to ¢ ight degrees In addition, the 


fect of aspect ratio on normal forces was determined with a 
f models having an aspect ratio range from 1.0 to 2.0 
The results were compared to linearized inviscid theory and indi 


ted normal forces which were 33 per cent to 77 per cent greater 
in values predicted by theory, the experimental normal forces 
] I 


lecreasing With increasing Reynolds Number The decrease in 


rmal force with decreasing aspect ratios was found to be from 
}to 7 per cent greater in magnitude than predicted by inviscid 
of Reynolds Number over 


e range 


tested 


Weak interaction theory was used to derive a simple expression 


for the normal force and was found to predict the order of magni 


ude of the incre ised normal force as well as the trend of the nor 
| force with Revnolds Number 
SYMBOLS 
platform area of plate, sq.in 
span of pl ite, in 
( constant in linear viscosity temperature relation dé 
fined by Eq (S 
( normal force coefficient 
Ng = normal force curve slope, per radian 
( = pressure coefficient 
D drag force 
d variable in weak interaction theory, defined by Eq 
Q 
slip parameter, (2/3 )(/ Re VW 
chord of plate, in 
Vf Mach Number 
\ net normal force, grams 
\ gross normal force, grams 
\ tare force, grams 
static pressure 
stagnation pressure 
dynamic pressure 
Re Reynolds Number based on chord of plate 
Rein free-stream Reynolds Number based on 1-in. charac 
teristic length 
constant in Sutherland formula, see Eq. (8 
thickness of plate, in 
1 free-stream temperature, “Rankine, or tangential 
force 
J stagnation temperature, “Rankine 


temperature of plate surface, “Rankine 
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= velocity, ft. /sec 


it = pan force, grams 
x = distance along plate, measured from leading edge, 
in 
a ingle of attack, radians 
r = ingle of attac k, degrees 
= ratio of specific heats 
Lu = absolute viscosity, Ib.sec. /ft.* 
i = kinematic viscosity, ft.?/sec 
? = density, lb.sec. /ft.* 
= Prandtl Number 
x = interaction parameter, JJ°yV/ C/y Vx 


Subscripts 


conditions outside the boundary layer far downstre 
from the leading edge 
= refers to upper surface of the plate 
refers to lower surface of the plate 
) free-stream conditions 


INTRODUCTION 


R' CENT INTEREST in high-altitude high-speed flight 
has given rise to considerable interest in the 
aerodynamic characteristics of bodies of simple geom 
interaction el 


etries in flows where rarefaction and 


fects exist. Aside from the obvious practical impor 
tance of obtaining aerodynamic data in flow regimes 
where these effects are present, there also exists a need 
for information and details to aid in understanding the 
fundamental nature of such flows. Theories account 
ing for the effects of rarefaction and interaction are 
largely in the early stages, so that empirical informa 
tion is required to assess the accuracy of these theories 

The purpose of this report is to describe the experi 
mental procedures and results of a wind-tunnel study 
of normal forces on rectangular flat plates in supersonic 
rarefied flows. The scope of this study was limited to 
a preliminary investigation of the effect of low Reynolds 
Number on the normal force and normal force aspect 
ratio correction for plates at small angles of attack. 

To provide a physical reference for the results con 
tained herein it may be remarked that the experimental 
but 


flow conditions duplicate the fluid dynamic not 


thermal) flight condition of a plate with a 10 ft. chord 
flying at a Mach Number of four at altitudes of 57-68 


miles. 


Theoretical Considerations 


rhe force normal to an airfoil in an inviscid flow may 
be readily obtained providing the pressure distribution 
the With a flat plate the 


problem is somewhat simplified because of the absence 


over surface is known. 


of surface curvature, and in this case an integration ot 








1100 JOURNAL OF THE AERONAUTIC 


the pressure acting on the surface of the plate yields 


the normal force directly. 


Busemann! obtained, for the isentropic approximation 
to the solution for the pressure coefficient C, on a 


plane surface inclined at an angle a to a uniform super 


sonic stream, 


VM — 1 
) | — 1)/2]M4* — 2M?+4+ 2; 
(iM? — 1)? 


LY ' 
a” 


where C, = (pf — fo)/(1/2) p0 V0". 

The normal force, V, acting on a flat plate of chord L 
and width b, may be written as the integral of the 
pressure difference between the upper and lower sur- 


faces 
*] 
N —b | (bu — pi) dx (2) 


where the subscripts uw and / refer to the upper and 
lower surfaces, respectively, and x is the coordinate 
distance measured from the leading edge. Using the 
definitions of the normal force coefficient and the pres- 
sure coefficient, Eq. (2) may be written as 


v/ 
Cx = KC oats Cc. \d(x/L) (3 
normal torce coefficient = 
N/(1/2) po Vo? bL 


C, = pressure coefficient on upper surface 


} 


where Cy 


C,, = pressure coefficient on lower surface 


By retaining only terms of order a in Eq. (1) and 
substituting in Eq. (3), the normal force coefficient for a 
two-dimensional flat plate in an inviscid supersonic flow 


may be written as 


This is the familiar result of linearized theory. 

The above expression for the normal force coefficient 
may be modified to account for finite plate size. In 
reference 2 it is shown that the normal force developed 
in the areas that lie within the tip Mach cones is half 
that developed over the rest of the plate. Accounting 
for this fact, the normal force coefficient for a finite 
aspect ratio plate in an inviscid supersonic flow may be 


written as 


Cn = la/V AM? — 1 [1 — (0.5 RV Mf? — 1) 


where MR = aspect ratio = b/L. 

Experiments have demonstrated that the inviscid 
theory is not adequate for hypersonic flow. In high 
Mach Number low Reynolds Number flows the bound- 
ary layer becomes quite thick due to the high stagna- 
tion temperature - static temperature ratio and the vis 
cous nature of the flow. The thicker boundary layer 
results in a zone of interaction between the leading 
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edge shock wave and the boundary layer. The pres 
sures in this zone are considerably different from thoy 
predicted by the inviscid theory. In addition to th 
zone of shock-wave - boundary-layer interaction, th 
boundary-layer growth causes significant changes 
the pressures in the inviscid external flow field down 
stream of the shock wave. The influences from th, 
inviscid external flow in turn feed back into the bound 
ary layer to affect its rate of growth, so that the results 
of these mutual influences are very complex. Thy 
problem requires a simultaneous solution of the bound 
ary-layer and external inviscid flow equations. 

external flow interaction has 
Maslen 


Lees and Probstein wer 


The boundary-layer 
been studied by Lees and Probstein,®" 
Mack,* Kuo,'’ and others. 
interested primarily in the interaction effect, whereas 
Maslen and Mack considered the effect of slip on th 
boundary layer. Maslen showed that with use of the 
conventional boundary-layer equations and slip bound 
ary conditions there was no first-order s//p correction to 
the pressure distribution along the plate, and he ob 
tained a pressure interaction correction (due to the 
displacement effect of the boundary layer) which 
agrees with the first-order correction of Lees and Prob 
stein, who used the no-slip boundary conditions 
On the other hand, Mack used an Oseen-type sli 
boundary layer in his analysis, the same as has been 
used with some success by Schaaf and Sherman? to 
predict the effect of rarefaction on the skin friction 
Mack found a correction to the normal force which 
falls rapidly to zero as the Mach Number becomes 
large. 

The theory which will be used here for comparison 
with the experimental results is the first order ‘‘weak 
interaction”’ result of Lees and Probstein, which is 
meant to apply far from the leading edge of a two 
dimensional flat plate in supersonic flow. Their re- 
sult, which is a limiting case derived within the frame 
work of boundary-layer theory for large Mach Num 
bers, is that the induced pressure due to boundary-layer 


growth can be expressed to first order as 
Po = po(l + ydox (6 


where /,,. is the actual pressure at the plate surface and 
p.. is the pressure at the plate surface, as determined 
from the inviscid flow equations, at a point far down 
stream where the flow deflection is the same as the 
surface deflection. The term x. is the interaction 


parameter defined by 
X M.2V C/V V.x/v 


where v kinematic viscosity. The subscript 

denotes values outside the boundary layer, which aré 
obtained from inviscid flow relations. The term C 
which is the constant in the linear viscosity-tempera 


ture relationship, is defined by 


em VE (EES) : 
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NORMAL FORCES 


Here 7, = surface temperature, °Rankine 
7... = temperature of fluid stream outside 
boundary layer, °Rankine 
S constant in Sutherland viscosity-tem- 


perature relations, (216° Rankine 

The term ¢, of Eq. 6 is, for Prandtl Number ¢ = 0.725, 

d 0.968 [(7./T.)/M.?] + 0.145(y7 — 1 (9 

and for the laminar boundary layer on an insulated 
plate 

Fo/s 2|M..*~V/¢ (10 


The normal force is then given by 


v/ 
fae5) Owe 
J 0 si 
p >| 1 -+ yd X dx 1] 
J0 
\ | p 2V L by X 


[he inviscid conditions on the upper and lower sur 
faces, indicated by subscripts « and /, are evaluated 
either from linearized theory or, more exactly, the 
oblique shock Prandtl-Meyer 
small angles of attack the two methods yield essentially 


and relations. For 
the same results. 

Dividing through Eq. 12 by the product of dynamic 
pressure and the area gives for the normal force coef 


ficient 


(2 qo) (P 


rhe first term in the above expression is the result given 
by inviscid theory. The term in the parenthesis is 
positive for a positive angle of attack, and thus it is 
seen that the weak interaction theory predicts an in- 
Calcu- 
lations made on the basis of the above theory showed 


crease in the normal force over inviscid theory. 


that the normal force varied linearly with angle of at- 
tack for angles up to 5 degrees. 

It will be recalled that the weak interaction theory 
is meant to apply far from the leading edge. Here, 
however, it has been applied to the entire plate surface, 
what the errors involved 


and it is natural to ask 


in this procedure are. One consequence is easily 
The zeroth order skin friction, in the 
Weak 


interaction theory predicts a correction to the skin 
and, thus, the 


demonstrated. 
absence of interaction effects, varies as 1/+/ x. 


Iriction which likewise varies as 1/7 x, 
first order interaction-corrected skin friction is nonin 
The theory, if 


nght up to the leading edge, gives infinite forces in all 


tegrable. weak interaction carried 


directions except the normal one, since only the induced 
pressure is integrable. The same holds true for Mas- 


len's results. The difficulty stems from the nonuni 


(2) 
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formity in the expansion at Q. It is possible that 
the integration of the induced pressure right up to the 
leading edge, as was done here, does not introduce a 
serious error because the contribution to the normal 
force from the region near x = 0 is small, but it is clear 
that the procedure cannot be justified rigorously 

As yet, a satisfactory theory for the region near the 
The 
theory,'® which was developed for 


leading edge is not available. Lees-Probstein 
‘strong interaction’ 
large values of /* + Re, does give finite values for all 
the forces at x 0. However, the strong interaction 
theory does not join on to the weak interaction theory, 
and it does not seem warranted to force the two theories 
into agreement by an arbitrary fairing procedure in 
pressure variation which 


order to obtain an induced 


will be associated with an integrable skin friction 


Previous Experimental Work 


The inviscid linearized theory was checked by 
Ferri. He measured the lift on various airfoils at 
Re = 10°, and found that linearized theory slightly 


overestimated the lift curve slope. 

‘he first attempts at checking the weak interaction 
theory were made by Bogdonoff and Hammitt’ in the 
Princeton helium hypersonic wind tunnel at Mach 
Numbers up to 15. Their initial conclusion was that 
the experimental data on induced pressure correlated 


well with the parameter .1/*, y/ Re, but their measured 


pressures were considerably greater (by 2-3 times 


than those predicted by weak interaction theory. 
However, it has since been discovered that this discrep 
ancy was due mainly to leading-edge effects. It is now 
clear that the inviscid pressure field generated at a 
blunt leading edge in hypersunic flow strongly influences 
the pressures on the surface of a plate for a consider 


able Kendall 


the induced pressures on plates with sharp leading 


distance downstream.'* measured 
edges at W/ = 5.8 in air, and found that these pressures 
were linear in A/*/+/ Re, and only 15 per cent higher 
than predicted by weak interaction theory, up to a 
value of M?/~/ Re = 3.2. 

The effect of leading-edge bluntness is now the sub 
ject of considerable investigation. Bogdonoff and 
Hammitt!* have reported more helium data for blunt 
edged flat plates which they correlated with the semi 


empirical equation for induced pressure 


App 0.0161V t/x}1 + 40V C’Re,} 


where Ap = Pw — P is the pressure increase over the 
inviscid value, ¢ the leading-edge thickness, and Re 
the Reynolds Number based on leading-edge thickness. 
rhe correlation is formed by linear superposition of the 
weak interaction induced pressure result and an em 
pirical equation for the inviscid leading-edge pressures 
At a Mach Number of 12 in helium, which is typical of 
the data reported in reference 5, the inviscid induced 
pressure equals the viscous induced pressure at a value 
Re, = 


larger. 


400; for greater Re, the inviscid pressure is the 
The Hammitt-Bogdonoff experiments of ref 


erence 5 corresponded to values of Re 1.000 or 
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greater; whereas Kendall's data were mainly in the 


range of Re, less than 50; it is easy to see why Hammitt 
and Bogdonoff found induced pressures considerably 
larger than predicted by weak interaction theory. 

For the experiments reported in the present paper, 
the values of Re, were of the order of 1-4. Leading 
edge thickness effects may therefore be expected to be 
negligible. 


DESCRIPTION OF EXPERIMENT 


Test Facilities 


The experiments were performed in the No. 4 Low 
Density Wind Tunnel of the University of California 
Low Pressures Project at the Engineering Field Station 
in Richmond, Calif. It is an open jet, continuous flow, 
non-return type wind tunnel driven by a five-stage 
steam ejector vacuum pump system. 

An axially symmetric nozzle, designated No. 8, was 
used for the present experiments. The design of this 
nozzle is presented in reference 12, and calibrations are 
given in references 12 and 4. The nozzle has an exit 
diameter of 4.31 in. At flow rates from 5 to 26 Ibs. of 
air per hour, it has a Mach number range from 3.7 to 
1.1 and a Reynolds Number range of about 920 to 
3,600 per in. The maximum nonuniformity of the 
Mach 
cent. 
of the nozzle boundary-layer thickness with Reynolds 
Number and, hence, Mach and Reynolds Number can 
Table | 


Number in the test section is about +0.5 per 


The Mach Number variation is due to the change 


not be varied independently in the nozzle. 
presents the characteristics of this nozzle. 


Models 


The models used in the experiments were designed 
to serve two distinct purposes. One set was designed 
to provide normal force data over a range of Reynolds 
Numbers and are designated Models No. 1 and No. 2. 
The other set of models, four in number, were designed 
to provide data on the effect of aspect ratio variations, 
and will be referred to as ‘‘aspect ratio’’ models; they 
are designated models A, B, C, and D. The aspect 
ratio models also provided normal force data over a 


limited range of Reynolds Numbers. The model di- 
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TABLE 1 
Nozzle Characteristics 
No. 8 Nozzle (Mach 4 
Static 
Flow Pressure, Isentropj 
Rate, Microns i Reynolds Core Dia 
Lb./Hr Mach No Hg R No per in in 
5.2 3.70 50 143 920 L.2 
10.3 3.91 74 132 1,620 1.9 
14 3.98 90 128 2,130 2.3 
20 $4.05 112 125 2 800 24 
26 4.13 134 121 3,560 ie 
TABLE 2 
Model Dimensions 
Model Span, in Chord, in \spect Rati 
l 0.400 0.253 1 58 
2 1.105 0.690 160 
\ 0.501 0.499 | 
B 0.604 0.490 1 23 
Cc O S04 QO 500 1.6] 
2) 1.002 0.501 » 1) 


inensions are given in Table 2. It may be remarked 


that despite the low aspect ratios used, the Mach 
Number was sufficiently high so that the Mach cones 
from the opposite wing tips did not intersect upstream 
of the trailing edge. 

The plates were made of phosphor-bronze stock 
0.009 in. thick, and the leading edges were sharpened 
symmetrically to a 20° included angle. One surface of 
each plate was highly polished to provide a reflecting 
surface for the light beam used with the optical system 
for measuring the angle of attack. The models were 
sting-supported at their trailing edges, as shown in Fig 
1, and were positioned in a vertical plane in the air 
stream. The relation between the model and _ nozzle 
is best seen in Fig. 2. Fig. 2 also shows the sting 
shield used with the models. The vertical portion of 
the sting passed freely through the sting shield and 
balance housed inside the 


transmitted forces to the 


aluminum enclosure shown. 


Balance Apparatus 

The balance used for measuring the extremely small 
forces was a specially designed microbalance developed 
by R. N. 


scribed in detail in reference 15, 


Latz. The balance and its design are de 
and only its general 


features will be described here. 
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[he balance in its present form is a single component, 
beam-type balance with all operations remotely con 
trolled \ schematic diagram of the balance is shown 
in Fig A flexure pivot is used as the pivot unit 
of the b ilance. The pivot consists of two jaws joined 


by means of two pairs of crossed wires. The upper 
‘aw is fastened rigidly to the base of the balance, and 
the lower jaw is free to rotate with respect to the upper 
iaW \ shaft is mounted on the lower jaw and passes 
through clearance holes in the base and upper jaw to 
receive the model. A the 
lower pivot jaw to provide means for force and dis- 
placement measurements. 


beam is also fastened to 
The flexure pivot has prac- 
tically a zero restoring moment and eliminates the hys- 
teresis and wear problems associated with knife-edge 
pivots 

Helical quartz springs are used as the force measur- 
ing element on the balance. They are relatively in- 
sensitive to temperature variations and can be selected 
to provide the desired sensitivity and force range. 

The spring extension can be measured on a mechani- 
cal counter mounted on the balance itself or by a 
counter outside the tunnel operated through a selsyn 
linkage system. The least count on the counters cor- 
responds to 0.001 in. spring extension. 

The null indicating system of the balance is a linear 
variable differential transformer (L1’7D7) which pro- 
vides an electrical output whenever the balance is not 
at the null, or the 
schematic sketch the null-indicating elements are in- 
A counterclock- 


force-measuring, position. In 
dicated on the left arm of the beam. 
wise rotation of the beam causes the core of the L17DT 
to move out of the transformer shell. This displace- 
ment of the core induces currents in the circuit of the 
unit which The 
sensitivity of the unit enables the balance arm to be 


are indicated on a galvanometer. 
positioned easily to within 0.001° of its null position. 
The balance, together with its associated equipment, 
is shown in Fig. 4. The first element of the group 
reading from left to right), is the vacuum tube volt- 
meter which measures the output of the audio oscilla 
The output of the 


The third element 


tor, the second element pictured. 
audio oscillator is fed to the LIDT. 
is the microammeter which measures the output of the 
LVDT. The the control 


unit for the balance which includes switches to actuate 


fourth element shown is 
the spring extension mechanism and the angle of attack 
system, together with mechanical computers to register 


spring extensions and angles of attack. 


Angle of Attack Mechanism 


mechanism is also shown in 


The 
Fig. 4 


mercial speed reduction unit with a 400:1 reduction 


attack 
It consisted of a turntable mounted on a com 


angle of 


ratio. Through gearing and a selsyn linkage system the 
angular position of the balance could be read remotely 
on a mechanical counter) to within 0.09 


lo provide an independent measurement of the 


angle of attack of the plates, and to check for deflec 


tions from a set angle due to aerodynamic forces, an 


ORC EB? 
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Fic. 3. Schematic diagram of balance 


optical system was installed in the tunnel. It may be 
consisted of a 
The 


beam of light from the source was reflected by the 


seen in the background of Fig. 2. It 
light source and a curved galvanometer scale. 


polished surface of the plate to form an image on the 
galvanometer With this the 
could be measured to within 0.03 


scale. system angles 


Description of Tests 


Prior to any set of tests on a particular plate it was 
necessary to establish zero incidence. Since the plates 
were symmetrical, the condition of zero incidence was 
established by adjusting the incidence until zero normal 
force was obtained under flow conditions. 

The desired flow conditions were established roughly 
by a rotameter reading. Finer adjustments were made 
by reading the stagnation manometer. The supersonic 
jet was then balanced, to free it from expansion waves 
or shock waves at the exit plane of the nozzle, by in 
troducing secondary air in the downstream manifold 
until the nozzle wall tap pressure agreed with the 
chamber pressure as measured on a pressure ring lo 


cated around the nozzle exit 


Tare Corrections 

Since the plates were supported by stings of finite 
size, a portion of the normal force measured each time 
was due to a tare force contribution from the sting 
To the 
models which were intended to be identical in every 


evaluate tare force contribution, tare force 


respect to the lift and aspect ratio models were fabri 
cated. These tare force models, however, were pro 
vided with a support rod which attached to the tip 
rather than the trailing edge as in the case of the gross 
force models. Details of the differences between tart 
force model support rods and gross force model support 
rods are presented in Fig. 1. 

Fig. 2 shows the arrangements used for the gross 
force and tare force systems respectively as well as the 


relationship of the models to the nozzle. It may be 
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noted that the dummy sting in the tare force installation 
does not experience the same flow that the real sting 
does during force measurement because the 
presence of the model support rod in the tare force 


unknown 


gross 


installation alters the flow to an degree. 
However, the magnitudes of the gross forces and tare 
forces are such that a 10 per cent error in the tare force 
causes only | per cent error in the net force, so that small 
inaccuracies in the tare force measurements may be 


tolerated. 


EXPERIMENTAL ACCURACY 


On the basis of the variations in the various nozzle 
calibrations and the 
could be measured, the maximum possible errors in the 
Mach Numbers and Reynolds Numbers are estimated 
The 


per 


accuracy with which pressures 


to be +1 per cent and +5 per cent, respectively. 
values of J/°/+/ Re are accurate to within +34 
cent. 

The estimated maximum and minimum errors in the 
gross normal forces and tare forces for the various 
models are given in Table 3. 

The maximum errors correspond to conditions for 
which the flow rate and angle of attack are a minimum 
Minimum errors correspond to conditions of maximum 
flow rate and maximum angle. 

The normal force coefficients depend on the difference 
between the gross forces and the tare forces, the free 
static the Mach 


Accounting for the functional relationship 


stream and free-stream 


pressure 
Number. 
that exists between these variables and the probablk 
error in the found that the 


minimum probable error in the normal force coeffi 


each of variables, it is 


cients is +3 per cent. The maximum probable errors 


for models No. 1 and No. 2 are +5.6 per cent and 


t4.0 per cent, respectively. For model A the maxi 
mum error in the normal force coefficient is about +4 
per cent. 

A least squares analysis of the normal force coefficient 
variation with angle of attack leads to the following 
estimates for the maximum probable errors in the 
values of the normal force curve slope, assuming that 


the normal force varies linearly with a: 


Model Max. Error, 
per cent 
l 1.7 
2 $1.9 
A > 0 


RESULTS OF EXPERIMENTS 


Examples of the basic data are shown in Fig. 5, 
which gives the variation of normal force coefficient 
with angle of attack and Reynolds Number, and in 
Fig. 6, which shows the tare force measurements on 
model No. 2. 

The major results of the experimental work are pre 
sented in Figs. 7, 8, 9, 11, and in Table 4. The com 
parisons of the experimental results to linearized theory 
(Fig. 7) use the theory in the form 


TIC 


LL 


SCIENCES DECEMBER, 1956 
TABLE 3 
Gross Forces fare Forces 
Max. Error Min. Error Max. Error Min. Err 
Model Pet cent Per cent Per ‘cent Per cent 
l 1.4 0.2 $0) J 
2 Pr 4 0.3 z.] 0.3 
\ + ] 0.5 2.a ) 
B o.2 0.4 2.0 18 
& De 0.5 2.5 i 
I) 1.8 0.2 2.5 0.6 
dCy o ba RSS 
7 (¢ Na) LT | - 
\ da /i1 V M | RV AM 
0.5 
(Cya)i7 } 
RV M / 


where the subscripts are defined as 
LT 
LI 


linearized theory for finite aspect ratio 


linearized theory for infinite aspect ratio 


The comparison of the experimental results with th 
weak interaction theory (see Fig. 8) also uses linearized 


theory as a standard of reference. For this comparisor 


the ratios 


(Cro) expr! (CNa)L7 and (Cnva)wr (CNna)z 


were used. In the second ratio the subscript I17/ 


the 
case according to weak interaction theory. 


denotes two-dimensional or infinite aspect rati 
The values of Cy, for all the models and flow cond 


tions are tabulated in Table 4. 


DISCUSSION OF RESULTS 


Normal Force Data 


The most interesting result of the experiments was 
the fact that the experimentally determined norma 


TABLE 4 
Reduced Data 





Model Re VW W3/V R Cva 
l 228 3.09 3.de 0.0317 
380 3.88 2.30 0 _O280 
531 3.98 2.73 0.0261 
683 $ O4 2 53 0. OP49 
S60 $1] aa 0.0245 
2 633 3.70 2 Ol 0 O60 
1,100 3. 90) 1.79 00 
1,405 3.97 1.66 0 
1,900 $05 1.52 0 
2,412 :.12 1.43 0 
\ 147 3.69 2.27 0 
803 3.90 2.09 0. 0244 
1,390 $05 1.78 0 O24 
1,751 4.12 1.67 0.0222 
B 17 1 B.d2 ? 36 0.0805 
784 3.90 2.%2 0.0257 
1,355 $05 1.80 0 O24 
1,705 4.12 1.70 0 0238 
Cc 144 >. 67 2.35 0.0302 
797 > 90 2 10 0.0267 
1 336 $ 04 1 S80 Q 0249 
1,747 4.12 1.67 0.0238 
D 178 3.41 2.0 0. O309 
797 3.90 2 10 0.0271 
1,389 t.05 1.78 0.0252 
1,765 1.12 1.67 0.02 
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force curve slopes were substantially greater than pre 


dicted by linearized theory. Fig. 7 shows the variation 


with Reynolds Number of the ratio between experi 


mental and theoretical normal force curve slopes, 


(Cva)exer, (Cna)ir- 
The ratio between experimental and theoretical nor 


mal force curve slopes decreases from 1.77 at a Reyn 


olds Number of 228 to about 1.33 at a Reynolds 
Number of 2,412. It should be noted that on Fig. 7 


the Mach Numbers for models with different chord 
lengths are somewhat different for a given Reynolds 
Number. However, the between 


the Mach Numbers at a given Reynolds Number is 


greatest difference 
only about 2 per cent. 

Solid curves are plotted in Fig. 7 for the data of 
models No. 1, No. 2, and C, which had aspect ratios of 
1.6. A single dashed curve is also plotted for the data 
of model A, which had an aspect ratio of 1.0. The fact 
that the data and C do not coincide 
shows that the inviscid aspect ratio correction does 
not hold accurately for the viscous flow situation, 
since if it did the ratio (Cy,)expr’ (Cne)tr would cor- 
relate data for models with different aspect ratios. 


for models A 


This problem is discussed in greater detail in the sec- 
tion on Aspect Ratio Data. 

It appears from Fig. 7 that there was an influence 
present which affected the relative order of the data. 
The relations that the curves for models No. 1, No. 2, 
and C not with 
changes in model dimensions; 
pect the individual curves for models with the same 


bear to each other are consistent 


in fact, one might ex- 


aspect ratio to coincide in their overlapping portions 
and form a single continuous curve. The relative 
order of the curves is difficult to explain in terms of 
errors since the data for any one model are consistent 
among themselves and were reproducible during the 
course of the experiments. The deviations of the curves 
are also greater than can be accounted for by experi 
mental accuracy and the slight Mach Number dif- 
ferences. 

One possible explanation for the deviation between 
the curves is that the sting effects varied with each 
model. 
metrically similar and, for model No. | particularly, it 
is possible that the sting had an appreciable influence 
Further experiments are needed to 


The model-sting combinations were not geo 


upon the forces. 
settle this question. 

The linearity of the Cy vs. a curves of Fig. 5 indicates 
that flow separation did not occur near the trailing 
edge. The pressure ratio Ap’/p necessary to separate 
a laminar boundary layer at Reynolds Numbers of the 
xX 10° is about 0.20. 
Ap/p varies as (Re,) 
0.75. 
shock with flow turning angle of between 10 and 11 


order of 5 The pressure ratio 
2,900, Ap/p 


an oblique 


so at Re, 


This pressure ratio corresponds to 


degrees, which is greater than the maximum angle of 


attack used in the experiments. The assumption that 


flow separation did not occur thus appears justified. 
Fig. S presents a comparison between the experi 


mental results and weak interaction theory. The 


rt 


\ 
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agreement is seen to be quite good. It will be recall 
that Kendall found induced pressures which were soy 
15 per than 
whereas here the theory is, for one modei, higher tha 
the data. 
normal force measurements 1s a “‘trailing-edge effect 


cent higher weak interaction  theor 


One factor which may have influenced th 
Expansion waves are formed at the lower boundan 
layer seam, and compression waves at the upper bound 
ary-layer seam, as shown in Fig. 10. The compressio 
waves feed high pressures upstream in the subsonj 
portion of the boundary layer on the upper surface 
and the expansion waves will propagate low pressures 


surface. Both effects 


upstream on the lower will 
tend to decrease the normal force experienced by th 
plate, with the effect diminishing with increasing 


Reynolds Number. It seems quite likely that this 
trailing edge effect might account for part of the dif 
ferences between the Kendall's 
data. A decrease in normal force of the order of 10 per 
cent due to trailing edge influences would be sutflicient 


present results and 


since, as will be seen in the discussion which follows 
the 
theory which was used to adjust the weak interaction 
The 
use of the experimentally deduced aspect ratio cor 


aspect ratio correction derived from linearized 


results for finite aspect ratio is in fact too small. 


rection would move the theoretical curves down, on th 
average, about 5 per cent. 

The experimental data can also be compared with 
the theoretical calculations of Mack. He derived ex 
pressions for the normal and tangential force and the 
moment coefficients for a flat plate at an angle of at 
tack, using an Oseen slip boundary-layer correction 
This boundary layer presents a sharp leading edge t 
the flow; the displacement thickness slope d6* dx at 
y = 0 in the Oseen slip theory has the value 2 3 J/ 
Mack found for the normal force coefficient the result 


Cy = (4a/V M? — 1)X 

f1 + (2/3MV AP — 1) exp (J?) erfe (J 14 
where / (2'3)(V Re AJ). At an undisturbed 
stream Mach Number |/ { and plate Reynolds 


Number Re 100, which are representative of the 
present experimental data, Mack's correction to C 
is of the order of magnitude of | per cent. It seems 
reasonable to conclude that in the present tests the el 
fect of slip was negligible, and that the boundary layer 


presented an essentially blunt leading edge to the flow 


Aspect Ratio Data 


The Mach Number and aspect ratio range used 
the experiment are such that the maximum difference 
in normal force coefficients between two models is & 
pected to be about 7 per cent, on the basis of the invisct 
correction, for models having aspect ratios of 1.0 an 
2.0 at a Mach Number of about 3.70. 

Fig. 9 presents a plot of the normal force curve slope 
for various aspect ratio plates plotted against the hy 
The data 


personic similarity parameter, 1/*/+/ Re. 
1.6 enable 


for models with aspect ratios of 1.0 and 
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reasonably good curves to be drawn. The data for 
models with aspect ratios of 1.2 and 2.0 appear to be 
somewhat less reliable and curves were faired through 
these data so that they were consistent with the curves 
based on the better data. 

The curves show that increasing the aspect ratio 
has the effect of increasing the normal force curve 
slope. The difference in the normal force curve slopes 


jor models with aspect ratios of 1.0 and 2.0 at J = 
3.70, (118 V Re = 2.30) is approximately 13 per cent, 
whereas the theoretical difference predicted for these 


conditions is approximately 7 per cent, as mentioned 


ibove. Thus it is seen that the inviscid aspect ratio 
correction term underestimates the effect of aspect 
ratio when applied to a viscous flow situation. This 


is as expected, since the theory confines the influence 
of the end of the plate to a region lying within the 
Mach cone originating from the leading-edge wing tip. 
[he flow outside these affected areas is everywhere 
supersonic according to the theory, so that it remains 
uninfluenced by the flow 
field. In the real case, however, the thick boundary 


laver over the entire surface of the plate provides a sub 


disturbances elsewhere in 


flow through which disturbances and 
this, 


sonic laver of 


pressure fields may be transmitted. Because of 
portions of the plate surface well within the region be 
tween the tip Mach cones may be affected by pressures 
it the plate ends, and the aspect ratio influence will be 
greater than predicted by inviscid theory 

lo enable a better comparison of the experimental 
ind theoretical aspect ratio effects, it is convenient to 
Eq. 5) by in 


generalize the aspect ratio correction 


troducing a function of Re, 


RV AP 1}} 


1D 


C la/~V M 1h 


f(Re 


has the value 0.5 for Re > Now form 
Ce, 


denote particular aspect ratios, we have 


where f( Re 


ing the ratio, Cy where the subscripts | and 2 


a 


Cy RV AP 1 | 


RV MM? 1] 


| f( Re 
Re 


Va lf 
By rearrangement of the terms there results 

RV M? — 1[(Cx,/Cve.) — 1] 

"Na, hi/Cn,. Me | | 17 


Now assume that for the range of Reynolds Numbers 


involved 


f( Re constant K 


rhis is almost a necessity at this point, since the scatter 
in the data make it impossible to obtain a check on a 
more complex, but more accurate expression for f( Re). 

From the faired data of Fig. 9 the terms of Eq. (17 
were calculated and plotted. The results are shown on 
For any given 
Each 


Fig. 11, which also shows the theory. 
pair of aspect ratios, four points are plotted. 
point represents a particular Mach Number and Reyn- 


» 


© 


f( Re) introduced in Eq. (15 
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olds Number, but the Reynolds Number variation far 
the Mach The Mach 
Number varies about +5 per cent around mean value, 


exceeds Number variation 
whereas the Reynolds Number increases fourfold as 
noted on Fig. 11. 

A straight a slope A 
through the data. The 
given pair of aspect ratios scatter about this line show 


line with 0.75 was drawn 


calculated points for any 
ing that in fact a Reynolds Number dependency does 


exist. However, for these experiments with their 


limited range of variables, f( Re 0.75 (approximately 
and the aspect ratio correction which multiplies the 
infinite aspect ratio normal force coeflicient to give the 
normal force coefficient for a finite aspect ratio plate, 


mav be written 
| 0.75 RV M | 


The above correction should not be taken as applying 
accurately for Mach Numbers, aspect ratios, and Reyn 
olds numbers outside the range of these tests 

It is expected that when more aspect ratio data be 
come available, a better aspect ratio correction will 
have to devised the effect of 
Mach Number in the Possibly, the ef 
fective width of the tip Mach cones at the trailing edge 


be which will include 


correlation. 


could be related to the boundary-layer displacement 


thickness 6,* there. Instead of the inviscid value 
L/V M 1, one could writ« 
& - 
ké,* 
VM | 
| L.73k(1 0.27.M 
Ld 


V Re, 


where & is an unknown “influence factor. With this 


formulation, we obtain 
} j | 
RV M | 


: L.73kR0 + 0.27172)V M ly 
4 (1S 
2 2V Re; J 


The quantity in the brackets replaces the quantity 
rhe limited range of the 
variables and the scatter in the data make it difficult 
to reach any conclusions concerning the value of 


in particular, whether it is actually a constant, and 


therefore Eq. (18) does not offer any advantage over 
the method used here for determining the aspect ratio 


correction 
CONCLUSIONS 


Measured normal forces on rectangular flat plates 
f 200-2,400 


and a 


over a Reynolds Number range « 
1.1 were from 33 to 77 per 


The 


experimental normal forces decreased with increasing 


Mach Number range of 3.7 
cent greater than predicted by linearized theory 


Reynolds Numbers. First-order weak interaction 


theory gave a good approximation for both the magni 
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tude of the increase in normal force and the trend of the 


normal force with Reynolds Numbers. No slip ef 
fects were discernible. 

The decrease in normal force due to a decrease in as 
pect ratio was found to be from 3 to 7 per cent greater 
than that predicted by linearized theory. Over the 
Reynolds Number range of 450-1,700 and Mach 
Number range of 3.7-4.1 the aspect ratio correction 
which multiplies the infinite aspect ratio normal force 
coefficient to give the normal force coefficient for a 
finite aspect ratio plate can be approximated as [1 — 
(0.75/.RV IP? — 1)]. 
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the actual maximum change in frequency squared to 
that based on elementary theory is 
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k? = 4fe(d) (A) (G/E)(b/ 
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PRESSURE, P 


SMALL-DISTURBANCE RELATION: 
P—Poo = Ro Iee W(t) 
SECOND-—ORDER RELATION: 
Ti) 


P— Post Poo sola) + (Gee 
LIGHTHILL’S THIRD-ORDER RELATION 
2 3 
P—Pon*Pael(at) +744 (at )+ BG) | 


Formulas for the pressure on a piston in a one-dimen 
sional channel 


2, 


Fic. | 


a, &B uncoupled natural frequencies of oscillation 
of spring-restrained typical section in 
bending, torsion, and flap rotation, re 


spectively 


INTRODUCTION 


5 ges TERM ‘“‘piston theory,’ as used in this paper, re 
fers to any method for calculating the aerodynamic 


loads on aircraft in which the local pressure generated 
by the body's motion is related to the local normal com 


ponent of fluid velocity in the same way that these 


quantities are related at the face of a piston moving in 
a one-dimensional channel. 
which such an analogy is valid are now well understood 
and are outlined below. In general, piston theory may 
be employed for large flight Mach Numbers or high re- 
duced frequencies of unsteady motion, whenever the 
surface involved is nearly plane and not inclined too 
sharply to the direction of the free stream. The fore- 
going shape conditions are fulfilled by all but the im- 
leading-edge) regions of 


mediate tip (and possibily 


supersonic wings. In most aeronautical applications 
the normal component of fluid velocity is the given 
quantity and the surface pressure is the unknown to be 
determined, so that a point-function relationship be- 
tween the two is a great convenience. These observa- 
tions, coupled with the fact that arbitrary small de 
formations and arbitrarily time-dependent unsteady 
motions can be taken into account, recommend piston 
theory to the aeroelastician as a powerful tool for 
analyzing high-speed problems. 

The authors’ main purpose is to focus attention on 
the remarkable simplicity and utility of the method by 
presenting a selection of new applications and compari- 
which were 


sons with previously published results 


originally obtained at the expense of much more 


laborious computation. The number of different ex 
amples which can be given in a single paper is naturally 
limited, but it is hoped that readers engaged either in 
the field of aeroelasticity or elsewhere will be encouraged 
to seek other uses of piston theory on problems of cur 
rent interest. 

Fig. | depicts a piston moving with velocity w(¢) in 
the end of a channel containing perfect gas, whose un- 
disturbed pressure, density, and speed of sound are p., 


RONAUTIC 


The circumstances under 


\ 
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Provided that the piston generates only 
ie, 


p., and a 


simple waves and produces no entropy changes, the ex 
act expression for the instantaneous pressure /(/) on its 
face is 
j ‘ (2 
PP ' \Y 1)/2|(w/a { 


Depending on the magnitude of the ratio w/a., Eq 
may be approximated by the linear relation, 
p p Pod W 


} 


by its second-order binomial expansion, 


p Pp p..a Z wea r ily + 1 ti(w a } 
or by its third-order expansion, 
fy p 
f 
pa »(@ a T Y l 1} ( w/a 
- ‘ i ‘ 
} 1)/12\(z 
Eq. (3) resembles Busemann’s quadratic formula for 


steady supersonic flow past an airfoil. Reasoning from 
a suggestion of Hayes,! Lighthill’ pointed out that Eq 

t) can be used with excellent accuracy even under non 
isentropic conditions to calculate the pressure on an air 
foil in steady or unsteady motion whenever the flight 
Mach Number has such an order of magnitude that 
M?> 1. 
cannot be too large, where for practical purposes 6 repre 
sents the larger of the thickness ratio of the airfoil or the 
ratio of maximum amplitude of unsteady motion to air 
foil chord length. 

It is of historical interest that a theory equivalent t 
applying the linear relation Eq. (2) to simple harmon 
oscillations of a flat plate with hinged flap in two-di 
mensional flow was given by Collar’ in 1944. Collar’s 
work consists of an extension by physical reasoning ol 
for steady 


formula supersonic 


Ackeret’s well-known 
flight. It was provided a mathematical foundation by 
Temple and Jahn,‘ who reduced the same results from a 
series expansion of the exact solution when 1/ — 

independent of the frequency of oscillation. (The con 
dition of large Mach Number also implies that the fac 
tor V Al’ l 
Collar’s formulation into exact 
In 1946 Jordan’ presented the high 


can be replaced by ./, which brings 
agreement with linear 
piston theory. 
Mach Number expansion for two-dimensional sinus 
oidal motion and noted the point-function relation be 
tween p and w which is the characteristic feature of pis 
ton theory. 

More recently Hjelte® and Landahl, Mollo-Christen 
sen, and Ashley’ have shown that linearized piston 
theory can be used for arbitrary small motions of thin 
two-dimensional airfoils whenever any one of the con- 
ditions M7? > 1, kAZ? > 1 or k? MM? > 1 is met, k being an 
appropriate measure of the ‘‘unsteadiness’’ of the flow 
The last of these three conditions suggests that certaitl 
subsonic flows may be covered, but comparison with 
more that satisfactory ac 


curacy will be attained only in very high-frequenc) 


exact solutions indicates 


problems such as local skin vibrations. 


An additional limitation is that the product 1 
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Lan l's systematic presentation’ of the circum 
stance nder which linear and second-order piston 
theor\ n be deduced from the general theory of un 
stead\ tential flow brings together the various lines ot 
pproach outlined above and furnishes a rigorous esti 
mate the magnitude of errors to be expected when 


Eg. (2) or Eq. (3 is adapted to practical problems 
He tre ul upper wing surface close to the X- Y-plane 
na uniform stream of velocity LU’ parallel to the 
lirectiol rhe instantaneous surface position is given 

4 bZ0(X, Y, f 5 


the parameter 6 being factored out so that Z) becomes a 
jantity of order unity when made dimensionless by 
division with a typical length-—1.e., the midspan chord 
Phe normal fluid velocity produced by the wing's 


motion 1S 
Ow’ (O Ot) + L(O OX) |6Z 8) 


By means of an iterative procedure Landahl derives a 


general solution for the pressure coefficient 


on the surface. When this is expanded in powers of 6 


ind 1, .\/*, a series develops in the following form 
( 6 MIC + (6 MW*)C + (6/M°*)C, T 

OW AM 6°C (6°/M?)C + O(67/ AT ‘ 
The various coefficients C,,,") and C are of order 


nity whenever the solution 1s uniformly convergent, so 


1 


that the 6 and .\/ factors indicate the magnitude of the 
rious terms. Those of importance to the present dis 


ussion are listed below 


( 2(2 [ Ya 
( y +t 2](wo l Yb 
fo O 5 
£) X 
re) Ov P 
wylé, v, (Ut 2b 
g dé Ve 
l 


\ll space coordinates in Eq. (9c) are divided by the 


typical length 20; x,(y 


gives the projection of the wing 
leading edge on the x-y-plane. 

The 6 J and 6° terms in Eqs. (S) and (9) are obv1 
ously identical to the linear and second-order portions 
ol Eq 5) 
contains absolute errors no greater than the orders of 


7 and 6 Mf? 
neither of the products 6.7 or ké\J can be too large, 


Hence, in general, the piston approximation 
There is also a general limitation that 


which means physically that the surface can nowhere 
be inclined too sharply to the direction of L’ and that 
normal velocities produced by the unsteady motion 
must be limited in comparison with the speed of sound. 
Because the mathematical solution is not uniformly 
convergent in the neighborhood of a cut-off wing tip, as 
might be expected from the failure of the one-dimen 
sional flow analogy when fluid can spill laterally around 


HEORY 








——)— 2 bT(x) 
U 


ELASTIC AXIS 


e 2 bx, + 
Fic. 2. Geometric properties of an airfoil section with flap 
the tip, the error estimates are incorrect there rhis 


does not forbid the use of piston theory to account for 
other three-dimensional effects, however, and the region 
where inaccuracies may occur is clearly bounded by the 
inboard Mach line from the tip leading edge. Finally, 
the third-order term in Eq. (4) cannot be predicted by 
Landahl’s analysis, which is based on a flow assumed 
free of strong shocks and entropy changes. The justi 
fication for including third-order contributions, when 
they are appreciable, is regarded as excellent but 


nevertheless empirical 


Two-DIMENSIONAL FLUTTER PROBLEMS 


Phe aerodynamic derivatives used for predicting the 


flutter characteristics of a two-dimensional airfoil o1 


“typical section’ are calculated by considering sinu 


soidal motion in the three degrees of freedom /: (vertical 


translation), @ (pitching), and 9 (flap rotation) illus 


trated in Fig. 2. The axis of pitch and the hingeline of 


the aerodynamically located at 


unbalanced flap are 
distances 2bx) and 2bx,, respectively, behind the leading 
Up to terms of second order, this motion gives 


rise to normal components of fluid velocity 


edge. 


Here the negative sign refers to the upper surface and 
the positive to the lower, since w must always be taken 
in the outward direction from the surface as used 1n 
piston formulas like Eqs. (2)—-(4 

Although arbitrary time dependence can be handled 
without additional difficulty, only the case of simpk 
harmonic motion is treated here, because it permits the 
introduction of the familiar notation of Garrick and 


Rubinow.’ For example, the lift and pitching moment 


about x vy per unit span, due to 
h he lla 
a ae llb 
y) Be Li¢ 


can be written in general 
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Note that & now represents the reduced frequency wh, l 
of the oscillation. 

When second-order piston theory, according to Eq. 
(3), is employed to approximate these dimensionless 
coeflicients, the following relatively simple expressions 


are found 


0 l5a 
L kM l5b, 
/ 1/k?M (15¢ 
/ | kAl y+ 1) R](A,,/80 (15d 
l v1) RM (7 1) Rk?) [r(x1)/4] (15¢ 
lie | \ RAI | y+ 1) /R\(Ap/S8d (151 
WM,’ 0 log 
Wo’ = (1, RM y + 1)/kR](A,/8b 15h 
V/ 1 kM y + 1)/k?](A,,/8b (151 
WV 1 BRM (4 1) R\(M, 46%) (15; 

Vs ate V1 2vo(1 v1) | RM 
7 1) /R7] y (x vo) [r(xy), 2 (A »/S8b?) 1 15k 

M, 
(4/3)( — Nv) D(x Vo) C1 vy bvgv (1 V1) 
kA 
y + | 

_Q | m - “ , a | (150) 


In each of Eqs. (15) the first term is the linear part, 
which arises even in the absence of thickness; each con 
tains the Mach Number as a factor /~—', in agreement 
with Eq. (8). The thickness terms, when not identically 
zero, are linearly dependent on the thickness distribu 
tion function 7{v). Most of them involve simple geo- 
metric properties of the profile, such as A, and M,, 
which are the area and first moment of the area about 
the leading edge. Similarly the flap derivatives are de- 
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Fic. 3. Elastic and inertial properties of a two-dir i01 
typical section 


pendent on flap cross-sectional area A » and first moment 
J, about the hingeline. Certain of these geometrical 
interpretations have been pointed out by Miles"’ in the 
discussion of second-order effects in his comprehensive 
monograph on unsteady flow. 

It is of considerable importance to the flutter analyst 
that, although Eqs. (15) are derived from a nonlinear 
theory, symmetries between the upper and lower sur 
face motions cause the lift and moments to relate lin 
early to the oscillations producing them. At the higher 
Mach Numbers, aerodynamic loads computed from Eqs. 

12)-(15) are known to agree closely with calculations 
by much more complicated second-order theories, such 
as that of van Dyke.!! One instructive comparison can 
be made in the case of single degree of freedom rota- 
tional instability of the whole airfoil or the flap: both 
theories predict the impossibility of such an instability 
and show the same influence of thickness on the damp 
ing moment associated with the rotation. Taking air 
foil pitching as an example, stability is dependent on the 
sign of the coefficient .1/;, being positive when \/, 1s 


positive. From Eqs. (15), the piston theory value 1s 


VU, = (1/RM)[(4/3 hyo + 4x 
(y + 1)/k V/ 77) 16 


where J/,,, denotes the first moment of the cross-sec- 
tional area about . vo. Eq. (16) is positive for all 
values of x» and all airfoils thin enough to fall within 
the scope of the theory. If the section has fore-and-aft 
symmetry (double-wedge, circular-arc, etc.), thickness 1s 
seen to make an unstable contribution to .1/, when the 
axis of rotation 1s ahead of midchord. 

For flutter computations involving more than one de 
gree of freedom, the airfoil is assumed to have the 
elastic and inertia properties shown in Fig. 3. The 
flutter equations are not reproduced here, since thet 
derivation and solution have been extensively dis 
cussed in the literature-—e.g., Garrick and Rubinow 
rhe significant point to be made is that, when the aero 
dy namic coefficients have forms as simple as Eqs D); 
implicit or explicit algebraic expressions can be worked 


ds 


out for the critical flutter speed LU’, and frequency 


functions of 7 and the other parameters defining the 
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PISTON 


system rhus, the trial-and-error process which was 


formerly characteristic even of two-dimensional flutter 
inalvses is eliminated. Flutter prediction becomes 
more elementary at high speeds than any other similar 


incompressible flow included-—in the aero 


elasticia The 
profile shape and initial angle of attack* can be intro 


calculation 


n's experience. influences of thickness, 


duced, for the first time, in a rigorous and routine 


fashion 


\s a first illustration, Fig. 4 shows exact and ap 
proximate curves of dimensionless bending-torsion (1 — 
a) flutter speed L’y bw, vs. M for a typical section with 
i low frequency ratio, elastic axis at mid-chord and cen- 
ter of gravity at 60 per cent chord. The notation is 
that of Garrick and Rubinow® and is reviewed in the list 
of symbols. Fig. 4 is taken from a previous paper by 
Heller and the present authors,!? in which the applica 
tion of piston theory to two-dimensional bending-torsion 
flutter is extensively elaborated. A conclusion stated 
there, to which no exceptions have yet been found, is 
that the simplified approach yields Uy bw, values within 
Q per cent of their exact counterparts for .\/ 2.0. 
Incidentally, third-order terms from Eq. (4) were car 
ried through the original calculations but have since 


where x 


So 
ve) 
~ 
al 
T 
—_ 
~ 
St 
~ 


lhickness enters only through the constants A and B; 
its influence may be strongly destabilizing, however, as 
suggested by the data for a 6 per cent double-wedge air- 
foil given in Fig. 4. Other interesting implications of 
Eq. (17) and various reduced forms thereof have been 
discussed previously. '° 

Figs. 5 and 6 compare piston theory results with exact 
calculations by Woolston and Huckel'*® for two-dimen- 
sional bending-aileron (i — 8) flutter of an airfoil with 
zero thickness and a 20 per cent chord unbalanced flap. 
In each case the ordinate, abscissa, and section parame- 
ters are chosen to correspond with those given by ref 
erence 13. Accordingly, Fig. 5 presents the dimension- 
less bending stiffness ba,/a., required to produce neu 
tral stability, plotted vs. static unbalance parameter x, 
fora system with unrestrained aileron (ws = 0) and the 
‘ther parameters stated. Fig. 6 presents the same or- 
dinate, plotted vs. frequency ratio ws w,, for two values 
of static unbalance. 

Although the Mach Number is only 2 


alue included in the exact calculations 


the highest 
the agree- 
ent with piston theory seems quite satisfactory. As 
shown by a few representative points, changes in the 
critical condition caused by replacing the flat plate with 


* Camber and initial angle are accounted for only by Eq. (4 
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Fic. 4 3ending-torsion flutter of typical section 
been discovered to have negligible influence within the 
ranges of variables treated. 

All the piston theory curves on Fig. 4, together with 
similar curves for any other combinations of the bend 
ing-torsion parameters, are calculable from the single 
closed-form solution 

x — lire(x — 1 7 
+ (1 — 2x) + A)? — [(43 bey + da B) 
$3) — 4xyp + 4p? + B 
= lNa 
Ivy + 4x0? + B 
a 6 per cent double-wedge profile are small in this 


No explicit formulas, comparable to Eqs. 
can be written for the critical conditions of 


situation. 
17)—(18), 
torsion-aileron, or bending-torsion 


From the practical standpoint, how- 


bending-aileron, 
aileron flutter. 
ever, the implicit formulas that must be solved demand 
little more computational effort than the bending-tor 


s10n case. 
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Jending-torsion flutter of typical section 


In connection with certain of the foregoing conclusions 
and also with flutter analyses on finite wings in the 
Mach Number range between 2 and 3, it must be pointed 
out that a theoretical inconsistency may arise from in 
cluding thickness effects. Examination of Eq. (S) re 
veals that, if the 6° term is retained while the 6 .\/* term 
is omitted, it is automatically implied that 


6>1/M 19 
This inequality is not fulfilled, for example, by a 6 per 
cent thick airfoil at A/ 2. 
influence of the 6 J/* term on flutter eigenvalues has 


A preliminary study of the 


demonstrated that its omission often produces com 
pensating errors in the elements of the flutter de 
terminant. This would be a fortunate conclusion, if 
generally valid, because this term is not of the piston 
theory type. 

Following a suggestion by Landahl, it is possible to 
introduce the higher terms of Eq. (8) into an iterative 
flutter calculation. For a given J/, a first estimate of 
flutter reduced frequency ky is taken from piston 
theory. Corrections based on this ky are then inserted 
into the flutter determinant, and a second approxima 
tion is computed with relative ease. It appears that a 
single iteration will normally be sufficient for most pur- 
poses, and this involves no especial difficulty. The 
6/\/* corrections to the aerodynamic coefficients for 


bending-torsion flutter are as follows: 


AL, sin 2k/2k\ (20a 
AL, = cos 2k/2kM' (20b 
AL,’ = sin 2k/4k* A (20c 
AL,’ (sin? k/2k® M8 (20d 
AM,’ = (sin 2k/kA*) + (cos 2k/2k2.M*) 
(sin 2k/4k*M*) (20e 
AM,’ (cos 2k/kM*) — (sin 2k/2k?M?*) + 


(sin? k/2k*M*) (20f 
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AM;' sin 2k/ 2k? .M sin? k/2k4*M 20g 


AM,’ (cos 2k/2k?M (sin 2k 4k1M 20h 


Fig. 7 illustrates the improvement in predicted bend 
ing-torsion flutter speed that is attained when piston 
theory is modified by a first iteration in the manner just 
described. The directions and magnitudes of — the 
changes in L’,/bw, thus found can be expected to turn 
out about the same, whether or not thickness effects are 


also introduced into the calculation. 


THREE-DIMENSIONAL FLUTTER PROBLEMS 


When adapted to flutter calculations for finite lifting 
surfaces, piston theory is automatically “‘strip theory, 
in the sense that the pressure at any point, as well as the 
loading along any chordwise section, is independent of 
\s dis 


cussed above, the loads are predicted with the least 


the motions of other portions of the system. 


accuracy near a cut-off tip, but this same criticism can 
be leveled at any of the other strip theories which have 
been used so extensively in the past. If the flutter 
problem is set up in terms of assumed modes and fre 
quencies of vibration, the principal novelty is that the 
generalized forces acting on the various modes can 
usually be written down in simple algebraic forms for 
arbitrary Mach Number and reference reduced fre 
quency. Certain unwieldy trial-and-error procedures 
common to more conventional methods are thus avoided 
or minimized, as with the typical section in the fore 
going discussion. 

The first example presented in this section is selected 
because unclassified experimental data are available and 
because it offers startling evidence of the importance of 


thickness in supersonic flutter. The case is one oi 
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Fic. 8. Bending-torsion flutter of supersonic wing -comparison 
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bending-torsion flutter of a straight-tapered wing of 
moderate aspect ratio,* whose other properties are 
carefully chosen to be representative of lifting surfaces 
i aircraft designed for appreciable supersonic speeds. 
[he actual properties cannot be revealed for reasons of 
security 

Fig. S presents theoretical curves of the dimensionless 
flutter characteristics vs. Mach Number. These were 
btained from the standard Rayleigh-Ritz equations, 
including one uncoupled mode of cantilever bending 
ind one The solid 


curve results from using the two-dimensional linearized 


of cantilever torsional vibration. 


aerodynamic theory of Garrick and Rubinow:’ at higher 
Mach Numbers it is seen to merge smoothly with the 
curve based on zero-thickness piston theory, as would be 
expected. The dot-dashed curve differs only in that the 
influences of profile shape and thickness have been in 
troduced through Eqs. (15). There is little doubt about 
which calculation agrees best with the measured data, 
even though .1/ is less than 2. The linearized prediction 
is dangerously unconservative, despite the fact that pre 
vious experience with lower Mach Numbers suggests 
that neglecting tip effects should lead to an error in the 
terms of the order 


opposite direction. Introducing 


6/\M/*, as outlined at the end of the preceding section, 
might be likely to displace the dot-dashed L’¢, bw, 
curve downward by roughly the distance between the 
two zero-thickness curves, thus giving an even better 
comparison with experiment. 

Quite a different sort of application of piston theory to 
three-dimensional wings concerns the possibility of find 
ing exact solutions for the damping and frequency of the 
aeroelastic modes, when the inertia and elastic proper 
ties are uniform along the span. A laborious calcula- 
tion of this type was once made, using incompressible 
strip techniques, by Goland and Luke.'* By contrast, 
the aerodynamic derivatives for arbitrary motion ac 
cording to piston theory are such that little more dif- 
ficulty is encountered than when analyzing modes of 
free vibration in a vacuum. 

One simple example which serves to illustrate the 
method is uncoupled pure torsion in a supersonic stream 
of a wing having constant chord 2), constant torsional 
stiffness GJ, and constant running moment of inertia J, 
ibout the axis of twist. The differential equation gov 
erning the spanwise distribution of twist @(Y, ¢) is 
L°b?/M) X 


GJ(0°6 OY" (Sp 


[,(0°6 Of 


(2b U7)(0@ Ot) [(1/3) — xo + x0?] + O[C1 2 vo] | 


In Eq. (21) the aerodynamic moment acting about the 


xis 2bxy behind the leading edge was found directly 


from Eq. (2), with 


W) - {U6 + 2b(06/0t)(x ‘ 22 


* This wing was flutter tested at M.I.T. as a portion of experi 


ental research conducted for Dynamics Branch, Aircraft Lab 
Wright Air Development Center, under Contract No 


AF33(038 


rator 


22955 


\ 


r 


HEORY 111 


























24 a 
06 , 
n=3 
——————— 
0 4 ; 30 I i 
05 ° er | 
A 2) 
e 3|3 
wh a 
04 26 
‘ ~~ 
oO a 
= a 
a O32 r22 
a F — 
© n= ~ S 
2 02 2 18 
a a 
= Co 
a uJ 
Qa e) x 14 + 
uw n=! ee ee 
—tT 
Ct) 10 
2 3 4 e 2 3 4 2 
M M 
Fic. 9. Frequency and damping ratios of two torsional aero 
elastic modes of uniform wing, ; 1.443, k 
025/M 


Thickness effects are omitted, and the minus and plus 
signs in Eq. (22) refer to the upper and lower wing sur 
faces, respectively. 

Eq. (21 
viously solved by products of exponential functions of } 
and f. 
root and zero torque at the tip force the } 


is linear with constant coefficients and is ob 


The boundary conditions of zero twist at the 
functions to 
be simple sine curves, as in the case of free vibration 
These conditions also yield a characteristic equation, 
from which can be computed the frequency and_ the 
ratio of actual damping to critical damping for the 
mode of each order. For the mth mode, the damping 


ratio and the frequency ratio prove to bs 


1/2ur. 1/3 iv, +- 4y 
S 23a 
V n2( MP, + [(] Dy ur,?|M 
A Wa (1 Wk \ 7 Wh, + ] Dy YT, V x 
V i ( nN e  & im 2b 


where k; is the reduced frequency w,/ lL’ corresponding 
to the fundamental mode of vibration in vacuum—+.e., 
21)V GJ I, 23c) 


Ww (7 


a 


All other quantities are defined in the list of symbols 
¢, is always positive, confirming the previously men 
tioned stability of pure rotational motions at high Mach 
Numbers, but, curiously, it turns out to be inversely 
proportional to » for the modes of larger order. Fig. 9 
shows curves of damping ratio and ratio of damped 
natural frequency in air to the fundamental frequency 
w,, both plotted vs. .\/ for the two lowest orders 1 
and 3. Numerical properties of the particular wing are 
listed on the Figure 


PANEL FLUTTER 


A problem which has caused analytical difficulties out 
of all proportion to its apparent complexity is the flutter 
of a flat skin panel exposed on one side to a supersonic 
air stream. Since most of the trouble stems from in 
ability to construct exact solutions based on available 
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Fic. 10 Pwo-dimensional panel adjacent to a supersonic air 


stream 


aerodynamic theories, it is reasonable to expect that 
piston theory can yield valuable information within the 
ranges of \/ and k for which it is valid. This is definitely 
the case, as will be illustrated here by studying the ele 
mentary example of the unbuckled two-dimensional! 


panel. 

The panel, pictured in Fig. 10, is embedded in an 
otherwise rigid, plane surface parallel to the X- 
axis. The two-dimensional displacement Z(X, ¢ 
occurs between supports at X = 0 and X = 2b. m 


represents the mass per unit distance perpendicular to 
the direction of the stream U; EI is the effective bending 
stiffness; and there is a uniform initial tension 7° per 
unit span which lies somewhere between the yield 
strength and the (negative) value which would cause 
buckling. The differential equation of motion, includ- 
ing air loads by linear piston theory, reads 


EI(0!Z/O0X*) — T(0?Z/O0X") + (p.U?/M)(0Z/0X) +4 
(m/2b)(0?Z/Ol?) + (p.U/M)(0Z/ot) = O (24 


Pressure variations on the back of the panel are omitted 
from Eq. (24) because they can usually be expected to 
increase the stability by removing energy from vibra 
tions. If the panel is simply supported, Z and 0°2/0..° 
must vanish at Y = 0 and 2d. 

The procedure for determining the stability of solu 
tions of Eq. (24) is toassume that Z(Y, f) is of the form 
| Z(X )e""|, and to try to find whether, for some com 
binations of the physical parameter /:/, 7, etc., it is 
possible for 2 to have a zero or negative imaginary part. 
If such is not the case, one is led to the conclusion that 
all small motions are damped. It will now be proved 
that for the membrane problem, in which the fourth 
derivative term of Eq. (24) is assumed to be negligible, 
Q can have only positive imaginary part, implying 
stability. 

The equation governing the complex dimensionless 


amplitude function 
B(x) = (1/2b6)Z(X/2b) 
of the membrane can be written in symbolic form 
[(d?/dx?) + a(d/dx) + (ce + id)|z = 0 (25 


where the constants a, c, and d are real numbers. The 


boundary conditions are 
(0) = 2(1) = O (26 
If 2 is separated into real and imaginary parts 


a(x) = 2,(x) + 22,(x) (27) 
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an elimination of z; leads to the following real differen 
tial equation and boundary conditions for 2 


[(d?/dx*) 4- a(d/dx) + c]*z, + dz, = 0 d ~ (0) (28 
s. = F ss 
) 9 ior Y Q, 1 (29 
|(d*/dx7) + a(d/dx) |z 04 
If the same procedure is used to eliminate 2,, it de- 


velops that the differential equation and boundary con 
ditions for 2; are identical with Eqs. (28) and (29 
However, it can be shown that Eqs. (28) and (29) have 
no eigenvalues if d # 0. If d = O, one obtains a very 


simple formula for the imaginary part of Q: 
Q; = podob/m 


This shows that the membrane panel is always stable, 
the degree of stability depending on the ratio of the air 
density to the running mass of the plate. <A few calcu- 
lations of damping ratios have been compared with 
numerical results from the exact solution of Goland and 
Luke at M = V8: satisfactory agreement is found 
in all cases. 

The foregoing proof holds true only in the case of a 
membrane.* When the fourth-derivative bending term 
is included one can demonstrate that, even though the 
separated differential equations for 2, and 2; are identical 
with the same boundary conditions, d does not neces 
sarily have to be zero, and the resulting equations will 
have two linearly independent eigensolutions for a given 
set of eigenvalues (c, d # 0). In this case, the sign of 
2; cannot be determined readily. 

Recently doubts have been cast on the validity of the 
Rayleigh-Ritz method as applied to supersonic panel 
flutter. 
exact solutions is easy to work out at Mach Numbers 


The relationship between Rayleigh-Ritz and 


where piston theory applies. As an illustration, a 
Rayleigh-Ritz calculation was made on a panel with 
membrane stiffness only, having (mA//p..b? 10 (this 
single parameter fixes the dimensionless stability prop 
erties). According to the exact analysis, the panel ts 
absolutely stable, having (2,b/ U’ 0.025. However, 
the assumed-mode procedure yields the dimensionless 
flutter frequencies and speeds, based on the fundamental 
frequencies w, of free vibration, which are tabulated 


below: 


Number of Normal Modes 
Used in the Analysis of a 
Membrane 


bs 1.58 +S] 
3 2.48 4. 82 
4 3.45 +. S4 


It is seen that taking a few modes results in a finite 
flutter speed, inconsistent with the exact solution. The 


* The authors are indebted to J. M. Hedgepeth of NACA and 
Professor Budiansky of Harvard University for pointing out this 
fact. The proof given in the original preprint for the plate-typt 
panel is incorrect. (For further discussion, cf. Readers’ Forum 
note by Hedgepeth.?') 
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PISTON 


values of l'y bw, might even seem to be converging 
toward a real limit, except that reference to the sequence 
of we #, Ss and careful examination of the flutter de- 
terminant for higher numbers of modes shows that some 
sort of divergence toward an infinite or imaginary LU’; 
ultimately will take place. Unfortunately, no such 
careful attention to convergence is observed in routine 
flutter calculations, so this difficulty could be com- 
pletely overlooked in the present case. It is interesting 
to note that if only odd or only even modes are con- 
sidered, the Rayleigh-Ritz solution will always indicate 
a damped motion as in the exact solution. 

rhe reader can imagine many other more complicated 
cases of high-speed panel flutter whose prediction will be 
expedited by the use of piston theory. Almost-plane 
plates with finite dimensions normal to the stream 
direction, and, in fact, with arbitrary shape, fall rigor 
jusly within its scope, since they do not have cut-off 
tips Nonlinear aerodynamic effects due to initial 
bucklings are readily included and fit well into an 
inalysis which also involves nonlinear structural proper 
ties. No further discussion is presented here on these 
ispects of the problem, because they are currently under 
study by Sechler and Fung at California Institute of 


Technology 
OTHER AEROELASTIC PROBLEMS 


Equations like (3) and (S) hold for motions with 
irbitrary time dependence, so that it is natural to con- 
sider employing them in unsteady aeroelastic calcula- 
tions where the displacements are not sinusoidal. 
lypical cases of this sort are the response of aircraft to 
gusts or blasts, and high-speed loads due to rapid 
maneuvers. Certain steady-state analyses may also be 
simplified because of the point-function character of 
piston theory, or rendered more accurate through the 
introduction of thickness effects. A good example of the 
latter kind is Biot’s section on nonlinearities in his study 
of the chordwise divergence of supersonic wings.'® To 
the authors’ knowledge, his is the first publication where 
Lighthill's theory is applied to a practical aeroelastic 
problem. 

Steady-state supersonic calculations with two-dimen 
sional aerodynamics are not facilitated by piston theory, 
because Ackeret’s formula for the supersonic airfoil al- 
ready exhibits the point-function relationship between 
pressure and normal velocity. A considerable saving is 
achieved, on the other hand, when three-dimensional 
ur loads are needed but when some inaccuracy in the 
treatment of the immediate tip region can be over- 
Weatherill’s recent thesis’ on the high-speed 
divergence and loss of aileron control of elastic wings 


looked. 


confirms this point very well; his use of piston theory 
came as a brief afterthought following an elaborate 
alysis based on three-dimensional aerodynamic in- 
fluence coefficients. 

To demonstrate the utility of piston theory as a tool 
lor dealing with unsteady transients, consider the gust 


response of a two-dimensional airfoil which can bend or 
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Fic. 11 Bending response of typical section to sharp-edged gust 


properties given in text 
twist but undergoes no elastic chordwise deformation 
Bisplinghoff, Ashley, and Halfman" 
such problems can be solved by appropriate superposi 


discuss how all 


tion of six indicial response functions, three lifts and 


three pitching moments. For example, the piston 
theory lift and moment following a step change / in the 
airfoil’s vertical velocity at / 0 assume the following 


constant values instantaneously 
Pr tp.bUhy M 30 


Av tp.b2Uhy MC 2) 31 


a 


Thickness effects are omitted from Eqs. (50) and (31 
but they can be introduced by means of a single multi 
plicative factor. Similarly, the lift and moment de 
veloped due to an encounter at ¢ = 0 with a uniform 


sharp-edged gust of vertical velocity wy, are 


eu ) —(4dp.. Ubw, A) (Ut 2b 0=< 07/6.< 2 G2) 
(—(4p.. Ubw,/ M 2< ltd 
\—(4e Ub*w,/M)((1/4) X 

M, = (Ut/b)? — (Ut/b)xo]) O< Lt b < 2 (33) 


(4p. Ub’w,/M) [1 —2xo] 2 < Utd 


At Mach Numbers above 2, Eqs. (30)—(33) check very 
well with the more exact supersonic indicial functions 
given by Chang!’ and others. To illustrate the use of 
these functions, the bending response to a sharp-edged 
gust has been calculated for the airfoil shown in Fig. 2 


with infinite torsional and control-surface stiffnesses 


(w, = ws = ©). The various constants entering the 
computation are VW = 3, uw = 33! 5, w, 100 rad. ‘sec., 
6b = 3 ft., U = 3,000 ft./sec. Fig. 11 presents this re 


sponse in the form of the ratio of the instantaneous 
bending displacement /(¢) to its final steady-state value. 
The curve is actually just the transient of a second 
order, linear system with constant coefficients. The 
slow rate of decay indicates that the aerodynamic 
damping at high Mach Numbers is quite small, so that a 
large number of cycles of vibration will go by before a 
significant amplitude reduction occurs. 

It should be evident that computations like the above 
can be performed on three-dimensional wings with 


other elastic and rigid-body degrees of freedom in 
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cluded in the equations of motion. The only increased 
difficulties are organizational and numerical; the gain in 
simplicity relative to more exact aerodynamic theories 
is proportional to the size of the problem. 


CONCLUDING REMARKS 


As was convincingly demonstrated by Bisplinghoff in 
his 1955 Wright Brothers Lecture,” high-temperature 
structural problems, boundary-layer heat transfer, and 
a host of related effects are certain to complicate the 
work of the aeroelastician in the not-too-distant future. 
It is, therefore, a pleasure to be able to report the prac 
tical success of one new element in the picture which 
represents a significant reduction in labor as compared 
with corresponding theories appropriate to flight at 
slower speeds. A Mach Number zone can be defined, 
lying roughly between 2.5 and the lower limit of hyper- 
sonic flow (a boundary determined by the size of the 
product 6.1/); it is characterized by the fact that thermal 
effects are usually quite important but that the dis- 
tinction can still be made between a boundary layer and 
an external potential flow. Within this zone it seems 
clear that piston theory is the logical tool to use 
for estimating aerodynamic forces when analyzing a 
large majority of all aeroelastic and aerothermoelastic 
phenomena. 
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Remarks on the Convergence of an Iterative 
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Remarks on the Convergence of an Iterative Procedure in 


I A RECENT research note,! Hua Lin examined the rate of con 
vergence of an iterative procedure in solving dynamic respor 


problems by plotting graphically the successive approximatior 





Solving Dynamic Response Problems. . Tse-Sun Cuo 1119 I 
of the following equation 
Axial Compression Buckling of a Pressurized Cylinder 19 , ; 
with a Thermally Induced Ring Compression . bat, ‘ ; 
LEONARD A. Harris” 1120 where & is the nondimensional displacement, £' its derivative 
with respect to the nondimensional time 7, and the damping 
ie es natasiséc - 2%e are = 
A Study of the Applicability of Rabotnov’s Creep Parame vette tee taltio’ conditions bein | £’(() 0 The Kes 
ter for Aluminum Alloy ’ ‘ 
FreD H. TURNER AND Kurt E. BLomguist 1121 itive procedure is defined by 
An Approximate Solution to Supersonic Conical Flow t(()) = ( < 
TENRY E. Fettis 1122 . " 
ind 
Erratum —‘‘Applied Aerodynamics and Flight Mechan- ‘ = 
ics’’ W. BaILEy OswaLp 1123 é r& = 1 2c | 
‘ t)} { 
Nitric Oxide Formation in Hypersonic Flow n=1,2.3. \ 
G. LOGAN, JR 1125 
where £, denotes the th approximation Lin made the observa 
On the Initial Vertical Descent of a Helicopter Following tion that the rate of convergence deteriorates rapidly with high 
Power Failure B. W. McCormick, JR 1125 damping ratio ¢and particularly at large values of time, 7 It i 
§ Stiff i matter of interest to examine further the validity of such 
Comment on the Reliability of the Strength of Stiffened procedure In this note we shall prove that, for this problem, 
Cylinders in Bending and or Compression 
BERTRAM KLEIN 1126 the successive approximations will converge to the actual solu 
tion for all values of ¢ and 7 It appears, therefore, for larg 
The Asymptotic Turbulent Boundary Layer values of ¢ and 7, high order of approximations are required 
R.A. Dutton 112% 3y treating the right-hand side of the differential equation of 


Effects of Free-Steam Vorticity on the Behavior of a Vis- 


cous Boundary Layer 


A Shock-Free Airfoil 


Eq. (3) as known functions of 7, the solution of Eq. (3) can be 


Fenn Wels 1190 put in the following form 


sf  ¢ m@.. B COS T 
H.C. Levey 1129 , 12 /du®) 2 qt fy 9 


i d-/d7 Li, | =(é 


ind, by applying the boundary conditions 


Eq. (4) is reduced to 





Phu 
Hence 
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Ho toto) ga f fh 


eae) Maal tel HOT WATER 
h) — golth)| dt OUTLET 


It is easy to find from Eqs. (2) and (3) that INLET fo 
e,  abliiaais ) ditt ciel 


CIRCULAR | HOT 
) SHELL —= WATER | -HOT WATER 
Ct 3 | | (+200°F// FEEDER 
custiely RING {pi 
)) into Eq. (8) we have 
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so that 





since r > U0. Substituting Eq. (1 


Sn(7 $n- 7 LOAD CELL — 














(25)" cf (r “)"-2 (1 + th) dh , —— : = 
1 2)! JO CLAMP —eo- hn UU NITROGEN 
pee (SPF) 
lias Pe 11) 

= Jay + nf r —t)"~*dt — Fic. 1. Schematic view of test specimen 
n 2)! 1 0 


~ 
i T t:)"—" dts length of 24 in., and a thickness of 0.0087 in rhe length of tl 
cylinder was perpendicular to the direction of rolling 
An aluminum ring was machined to provide a room temper 
Thus, the series ature force fit on the outside of the cylinder \s indicated i1 
Fig. 1, the ring, which was very stiff compared to the thickness 
S* fe,(r) E,-1(7)j the shell, was not attached to the shell but bore against the shel 
n=1 at mid-height with a moderate force fit pressure. The ends of 
the cylinders were clamped by straps to the 3-in. thick, circular 


converges uniformly and absolutely for all given values of ¢ and . : 
iluminum head plates of the test jig. When placed in the jig 


rand lim &,(7) = &7) is the solution required of Eq. (1). 
u—a the weight of the ring was supported by three wires which wer 
connected to the frame of the test jig 
REFERENCE The evlinder was heated internally by circulating hot water 
Lin, Hua, On the Convergence of an Iterative Procedure in Solving Dynami inside the tank at a temperature of 200°F \ pipe directed 
Response Problems, Readers’ Forum, Journal of the Aeronautical Sciences constant flow of hot water against the skin which was in contact 


Vol. 23, No. 4, p. 391, April, 1956 _ . P ° , 
aie . : with the ring at the mid-height of the cylinder When a te! 


perature of approximately 200°F. was attained, liquid nitroger 


+ ° 

Was poured into the trough around the external ring. Befor 

‘ - ‘ applving the axial compressive loading, the cylinder and ring 
Axial Compression Buckling of a Pressurized des Sap, poet re somal = 
. ° were allowed to reach equilibrium temperature At equilibriun 
Cylinder with a Thermally Induced Ring , aati ; : ae vo 
Cc . the skin temperature 1/4 in. from the point of contact of th 
ompression ring and the shell was about IS0°F., and the ring temperatur 


averaged about —200°F. across the section. The thermal de 
Leonard A. Harris * 
Senior Structures Engineer, Missile Development Division, 
North American Aviation, Inc., Downey, Calif 
July 12, 195¢ 


ee: AERODYNAMICALLY HEATED, supersonic aircraft 
structures must resist thermally induced stresses and de 
formations in addition to those stresses and deformations caused 
by external loads. The effect of the thermal stresses superim 
posed on other stress conditions is not always readily apparent 

The experimental investigation reported herein was performed 
to establish the effect of one particular thermal condition on the 
axial compressive buckling strength of a ring-stiffened, pressurized 
eylinder. The thermal condition investigated is one in which a 
circumferential line load is induced at the mid-height of the 
cylinder by refrigerating the ring and heating the skin. The 
ipplied load is a uniform axial compression. This combination 
of thermally and externally applied stresses would simulate a 
flight condition in which the skin is aerodynamically heated and 
the ring is kept cool by immersion in fuel 

\xial compressive buckling stresses were obtained from two 
model cylinders which were rolled from half-hard 18-8 stainless 
steel (E 27.0 * 108 psi The cylinders were fabricated with 


i single longitudinal seam weld and had a radius of 8.75 in., a 


* The author gratefully acknowledges the assistance of his colleagues of 





the Stress Design Data and Methods Unit and of the personnel of the Missile 


Fic. 2 3uckle pattern on tested cylinder 


Test Laboratory 
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iluminum ring with respect to the steel cvlinder 
oP vas great n if the ring and cylinder had been made of th 


because of the greater thermal coefticient of con 
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cti uminum as compared to steel s 
\fter attaining temperature equilibrium, the internal pressure 
URE vas inc! d to 14 psi, and a uniform axial compressive loading 
b central internal strut until the cvlinder buckled Z 
| rit buckling stresses for the two evlinders were 11,500 i 
50 psi. Fig. 2 shows one of the cylinders after re , 4 
from the test jig and after the ring had been pressed off ° 
[wo distinct buckle regions are visible rhe larger buckles, } 
me : ‘ . ‘ "3 
ATE hove the mid-height of the cylinder, are the axial compression tf + —- 

ER pes é : . } 
ER buckles which had wave lengths of about 1 in. in the longitudinal % 
E direction and about 4'/. in. in the circumferential direction 

Before the load was removed, these buckles were of the diamond LY ara 
4 shape typical for compression buckles of unstiffened cylinders 
' ‘ , . . 
The smaller buckles, which run circumferentially around th 
. tank at mid-height, were caused by the thermal expansion of I] 
4 
EN the skin and contraction of the ring It is believed that these 
F) smaller buckles were present before the axial compression loading / 
id been applied. The circumferential buckles were about / 
1'/,in. from crest to crest and about 1 in. wide 
Unstiffened cylinders of the same geometry had previously been 'A | 
ngth of the tested at room temperature under axial compression with internal E - 
; ‘ ; 6£<(,<t,<< G S7TRA/N 
pressure Comparison of the present tests with the unstiffened o . * 

mm temper evlinder tests indicated that the ring-stiffened cylinders failed at Fic. 1] Family of isochronic stress-strain curves and mastet 
ndicated ir stresses Within the lower portion of the scatter band obtained curve 
thickness of from the tests of the unstiffened cylinders. Therefore, this 
ist the shel thermally induced circumferential line load did not appreciably 
The ends of lower the axial compressive buckling strength of these cylinders 

; ; 2 — ’ : > 40 
ick, circular \ full-scale structure of similar geometrical configuration might pias 
| in the jig likewise be capable of resisting the effects of this type of thermal ~~ | 
é , 2) 
which wer condition superimposed on an axial compressive load » 
v 
. © aa 
. 
r hot water 9 ‘ 
- directed a + 4 J4 
s in contact 30 “a 
hen a tem 
id nitroge A Study of the Applicability of Rabotnov’s 
1g. Belor Creep Parameter for Aluminum Alloy 20 
‘r and ring ; 
quilibriun f T . / : 5 
quilibriun Fred H. Turner and Kurt E. Blomauist 
tact of th tre fhice AABA rcratt Company a nkoping, Sweden 0 4 
emperatur ¥19.195 
hermal de 
EFERENCE 1, which is an NACA translation of a paper by o 
R™ 7 . PS eee oe o Oz 04 | O« 08 - , ) 
\ N. Rabotnov, has apparently not been observed by Os S/TAIM (fe 
creep specialists in Western aeronautical circles. This note is . on 
Ve/ve vsed for determining rhe conston?r fo] 
intended to show that at least one expression due to Rabotnovy 
can give useful results for transforming creep curves from one Fic. 2. Stress-strain curve for 245-T3 Alclad Sheet 450°! 
stress level to another. This method is a useful complement to 
the well-known Larson-Miller parameter which permits trans 
formation of creep data from one temperature level to another raps 
31 
It has been shown by Rabotnov! that it is possible to trans 248-T3 Alclad Sheet 450 


form a family of isochronic creep curves to a single master curve 
by using a suitable parameter. This parameter is Dete from WAC TR 


53-24. Part / 








ge =o(l+as-t 7 = const 
Time fo OS % Predicted 
is a constant which is characteristic for the material and Ofress [psi Fetal strane vo/ves 
temperature in question, o is stress, € is strain, ¢ is time, and TJ is (Aewrs) (Aevurs) 
temperature Thus, for every given value of the strain e, the a ] T 
quation shows the relation between the stresses o; and the cor ‘ 5000 86.0 
responding times ¢ When / = 0, ¢ = o—e.g., the master curve 


/7000 37.0 4/4 


uncides with the tensile stress-strain curve at high strain rat 


see Fig. 1 20000 16.0 /44.0 


. 22000 6.0 66 


\) Given a tensile stress-strain curve and one point on a 
25000 RS /9 


eep curve at a certain stress, Rabotnov's parameter provide 








i possibility for predicting creep curves at other stresses. Cree] 30000 O15 ai 


results in reference 2 have been used in order to illustrate this 
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point. The first and second columns in Table 1 have been taken 
from reference 2, and the third column has been tabulated a 
follows. The first line in Table 1 gives the time 86 hours for 0.5 


er cent total strain when the stress is 15,000 psi The stress 


34,000 psi gives 0.5 per cent 
is the 


] 
strain curve in Fig 
an ordinary 


2 shows that ¢ 


tensile test Using this value of o % 


n the equation enables us 


strain in 
corresponding g-value 1 to determine 


is follow 








34,000 15,000-(1 + a-86 a = 0.29 (hours 
[he times for the same total strain at other stresses have been 
calculated, and a comparison with actual results is shown in 
Table 1 

4 lA } fy Ve v 
x 
Ld é 
" 
- O HOURE 
} / ~ 
- 20 HOURS 
9, i / ra 
a) x / / 
/ / 
/ 
_ 4 
STRA/N 
Fic. 3 Evaluation of constant a and construction of master 
curve for 248-T3 at 400°F 
98 
wy 
WK 40 HOURS 
e ye 
K J 
aa 20 E 
STRA/IN 
preaicled curve 
CX DeCl//I7C/ITO/ CUS VE 
a 
Vv) 
v) 
PS ” 80 HOURS 
iS J 
4 y 
20 ; 
if Os fo) 5 
STRAIN 
1 $, Comparison of extrapolated values with actual test 
results for 248-T3 at 400°F 
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TABLE 2 
248-T3 400°F. 


\Data from AST/V Vo/ 49 i 
-— : EE a a cre DS 
| z , 
| /o%o/ Strarry me O7ress vo/ve $ 
o) ) 
|e (ours ( psd) (Ps? 
O54 /OO 20000 | 20400 
0.30 /00O /5 000 1/4600 
/.¢ 200 20000 20800 
0.35 200 ‘5000 4400 
O 4% 400 5000 /5000 
OS 600 5000 “4700 
€ 3SOC 45000 4700 
B) A tensile stress-strain curve may not always be obtaii 
ible In this case the master curve can be determined from t 
isochronic stress-strain curves, where the times can be rather 
short \t any given value of the strain, corresponding point 
on the two curves make it possible to calculate the constant 
enables the master curve to be constructed (Fig 


which 


Some slight variation in the expected at vari 


value of a may be 


value has been used for the 





ous values of the strain \ mean 


whole of the master curve r 
It is now possible to predict stress-strain curves for other 
rable 2 together 


in Fig. 3 
his is shown in Fig. 4+ and 


longer times 
The master curve of Fig 


the actual test results 3 is used in this 
calculation 


These investigations cannot be claimed to prove the universal 
applicability of Rabotnov'’s parameter, but the results must b 


considered satisfactory for aluminum alloy of 24S-T type within 


the temperature-time-stress region of interest in current aircraft 
design 
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353, April, 1953 
Guarnieri, G. J Inte j t Stre Heati lest 1 
Structu Vetals, WADC Technical Report 53-24, May, 1954 
Dorn, J. E and Tietz, T. E., Cre and Stre Rupture I 
ASTM Proceeding Vol. 49, 1949 
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An Approximate Solution to Supersonic Conical 
Flow 


Henry E. Fettis 
Aeronautical Research Laborator Wright Air Developmer 
Center, ARDC, Wright-Patter AFR. Ohio 
July 23, 195¢ 
5 ee EQUATIONS which determine the supersonic flow field 
in adiabatic gas about an unyvawed cone may be writte! 
d*u/dé?) + u = lau 4 cot 6 } 
z= |( 1)/2] (c? { v2), v = du/dé\ 
radial 


In the above equations, u is the velocity component in t 
normal direction, @ the ang 


component in the 
local speed 


radial direction, a the 


direction, the 


between the cone axis and the 
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1 constant which, for any given problem, may be 
(See Thus all 


and a may be regarded as referred to c as 


ler s known further, reference 1 


qu tities u, UV, 


The standard technique for handling the system of Eqs. (1) is 


nt te numerically starting with the conditions 


ul @. = Uy! 


4 = () \ 


the semiapex angle and u is some number between 0 


it The 
tion ee reference | 


btained for the Mach Number of the free stream 


1 


solution is continued until the shock condi 
ire Satisfied. At this point, a corresponding 
This 
procedure has been carried out in reference 1 for numerous values 
rhe results show that, for small apex angles and 


Ma Numbers above 2, the shock conditions are met before the 





has become appreciably large In such cases, it 


be seen that if the first equation of Eqs. (1) is rewritten in the 


10 tL cot O(du/dé ~- 2y = 
a~ 1 + peot é 3 


id that if is small compared with a? (this is certainly true 
near the cone wall), then a reasonable approximation can be ex 
pected if Eq. (5 is replaced by 

d2u /dé?) + cot 9du/dé) + 2u = 0 j 


ubject to the conditions 


The general solution of Eq. (4) 1s 


( = Acosé + B}1 + cos 6 ln [sin 6/(1 + cos 6 6 
ind, by differentiation, 
+ 1 sin @ + Bicot 6 sin @ In {sin 6/(1 cos 6 7 
Phe initial conditions of Eqs. (5) give 
1 = uy (cos 6 sin? @ In [sin @)/(1 + cos 4)]| {1 
' (8 
b = My sins 6 \ 
Eq. (4) is the same as would have been obtained if the fluid had 


been assumed incompressible (a? = Nevertheless, the solu 


tions of Eq. (4) do not deviate appreciably from Kopal’s values 


even when v itself is no longer small. In fact, by comparing the 


solutions of Eq. (4) with the more exact solutions, a maximum 
error of 1.5 per cent of the total velocity in the #-component and 
5.3 per cent in the v-component was found to exist over the com- 
plete range of the table rhe greatest errors were found when 
the free-stream Mach Number was close to unity and occurred 


it the largest value of 6 in each case For Mach Numbers above 





maximum error was found to be about 2 per cent 
rhe approximations given by Eqs. (6) through (8) represent a 


substantial improvement over that given by the isothermal 


solution 


TABLE | 


Isotherm 
Solution 


Numerical 


\pproximation of 
Solution! j 7 


Eqs. (6 y 


Gg = 5D uy = 0.4, uz 10] 
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9 , a li = Uu 
5 0.40000 0 0. 40000 0 0.40000 0 

20 0.38144 0.12973 0.38149 0.12939 0.38637 0.10353 
3 0.33419 0.22849 0.33437 0.22766 0.34641 0.20000 
”) 0.26333 0.30961 0.26380 0.30821 0.28284 0.28284 
65 0 95 0.36935 0.17480 0.36711 0.20000 0.34641 
o 1.07206 0.40520 0.07379 0.40073 0.10353 0.38637 

% = 20°, uw = 0.8, M= 4.15 

ye) O SO0000 () 0. SO000 ) 0. SOOO0 0 

2 1.7 0) 0.79975 0.02729 0.79988 0.01396 
22 7 0 0.79904 0.05: 0.79951 0.02792 
y ) 7 0 0.79789 0.07859 0.79899 0.04187 
24 796! 0.10504 0.79632 0.10270 0.79805 0.05581 
25 79422 0.13050 0.79432 0.12608 0.79695 0.06973 


FORUM 


= My Sin (6 4 


rhe improvement is greatest in the v-component which, according 
to the isothermal solution, may be in error by as much as a fact 
j solutions f 


of two lable 1 gives a comparison of the three 


two typical sets of initial values 
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should be « 


O> PAGE 473* the paragraph including Eq. (4 \ 
led il and meat f 


panded to include both the loc n mean coefficients of 





turbulent skin friction as shown below 
shown approximately in Fig. 2 and given by the relations 


for turbulent skin-friction coefficients as follow 


Loc il ( 


] - Ld 4.15 log Rye } 
Mean, C 
0.242/Cy** = logy (Ry kK. E. Schoenherr $+] 
An increment 
Journal of the Aeronautical Science Vol. 23, N Py 169-484 
May, 1956 
t See also errata givenin the Journal of the Aeronautical Science \ 


23, No. 9, p. 892, September, 195¢ 


Nitric Oxide Formation in Hypersonic Flow* 


J. G. Logan, 
Head—Theoretical Flow Physics Section, Aerodynamic Researct 
Department, Cornel! Aeronautical Laboratory, Inc., Buffa 





July 23, 195¢ 


g ee HIGH TEMPERATURES resulting from high-speed flight wit! 
in the earth's atmosphere are sufficient to cause partial dissoci 
nitrogen and oxygen molecules as 
Between flight Mac! 
Numbers of 12 and 20, the high temperatures achieved behind 
normal shock waves may lead to the formation of NO, NOs, NeO, 
and NO 
ciated with the process of dissociation 
and N 
Wood.? 
the degree of 1onization 
Apart from the 
ind rotational states of the air constituents at high temperatures 


ation and ionization of the 
well as the formation of new components 


itomic species usually asso 
namely, O, N, O-, O 


Bethe and Teller! and 
wccording 1 


in addition to O.~ and the 


which have been studied by 
Also, of course, free electrons are present 
thermal excitation of the higher vibrational 


any variation in chemical composition must obviously influence 


the aerodynamic characteristics of the flowing medium throug] 





changes in molecular weight and specific heat.*»* In assessing 


the effect of these chemical changes, one must know the equilib 


rium concentration of the various constituents However 


it is equally important to have knowledge of the rate at whicl 


the real gas equilibrium is achieved since understanding of ac 


* The investigations of the NO reaction described were sponsored by tt 
Air Research and Development Command under Contract AF 18(600)-1332 
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dynamic heating effects and boundary-layer phenomena depend 
to a large extent upon the rates of dissociation and recombination 
Relatively small concentrations of components may have large 
effects, as for example, heating produced by recombination of 
oxygen atoms at a surface.® 

From the analyses of Bethe and Teller! and Wood,? one would 
conclude that the rates of dissociation of O. and Ne are such that 
the presence of atomic oxygen may be neglected below 3,000°K 
(5,400°R.) and atomic nitrogen below 6,000°K. (10,400°R 
below flight Mach Numbers of 10 and 16, respectively, because 
of the long relaxation path length required. At these temper 
atures, free oxygen and nitrogen are present in approximately 
at pressures 


1 per cent concentration in thermal equilibrium 
that, 


Bethe and Teller assume 


than 10 atmospheres 
the final equilibrium concentration of 


of the order of 1 per cent 


greater 
when the dissociated 
constituent is small the time re 


quired to reach the equilibrium state is large since two molecules 


must collide with a relative energy at least equal to the disso 
ciation energy When the final concentration of the dissociated 
constituent is small, the number of molecules possessing energ) 
would be small 


the dissociation also 


Consequently 


sufficient to bring about 
and effective collisions infrequent 
cluded that the relaxation time or corresponding relaxation path 
at least of the order of the projectile or 


they con 


length would be large, 
missile length 

Recent experiments at Cornell Aeronautical Laboratory, Inc., 
suggest that the relaxation time for the formation of free oxygen 
mav be theo 


much shorter than one would conclude from these 
retical investigations, and, consequently, it may not be possible 


AL 
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temperature i the 


to neglect the formation of free oxygen at 
h 


even if the conclu 





order of 3,000°K Nevertheless, 
Bethe and Teller and Wood are correct in regard to f lati 
of free oxygen and nitrogen, additional thermal reacti CL 
which may have significant aerodynamic effects at 3,000°K. |] 
particular, the reaction No + O, s 2 NO occurs at temperature 
below 3,000°K., and 1 per cent concentrations can be obtaing 
it temperatures of about 2,000°K Investigations of the rea¢ 
tion rates for the formation and decomposition of NO are being 
carried out at Cornell Aeronautical Laboratory These ¢ 
periments are being conducted in a shock tube designed 

i single controlled pressure and temperature pulse of é 
milliseconds’ duration Fig. | When the reactant g Ir, 1s 
processed by a pulse of this type at the end of the constant te: 
perature heating period, the gas is cooled rapidly by an « j 
fan at the rate of about 500,000°K. per sec 

This cooling rate is so rapid that the reverse reacti ,O— 
O N»2 cannot proceed fast enough to achieve the equilibriu 
concentration of NO at each instant during cooling Thi 
temperature is reached with a gas which retains mu ft 
NO corresponding to the high temperature Because t 
verse NO reaction (as well as the forward reaction 

it temperatures below about 1,500°K., the high-t eratur 
NO concentration is, in effect, partially frozen and can bi 
jected to leisurely chemical analysis at room temperature 

Fig. 2 shows the result of some recent measurements of nitri 

xide formation from ait Equilibrium amounts of nitric i 
were not achieved in the low-temperature region of t Irv 
because the reaction rate and heating time 1 At 
higher temperatures, the reaction rates were large, d, althoug 
equilibrium amounts of nitric oxide were appare f ik 
the cooling rates wert inadequate to freeze the igh-temperatu 
equilibrium. From studies of this type, it is possible to det 
mine the reaction rate and the equivalent reaction path leng 
At 2,000°K. (3,600°R the reaction time was calculated t 
ipproximately 70 millisec. decreasing to approximatel 
millisec. at 2,800°K. (5,100°R At approximately 3,000°K 

5,400°R the reaction time becomes sufficiently short-—less 
than one half a millisecond—so that the full equilibrium con 


centration can be achieved behind a shock wave ahead of a blunt 


body 


Since the air just outside the boundary layer of a blunt body 


traveling at hypersonic speeds may be processed by a tempera 


ture pulse similar to that employed in the chemical shock tube, 
it is possible for the full equilibrium concentration of NO attained 


behind the shock to be frozen as the air particles accelerate away 


from the stagnation point. Fig. 3 shows the temperature pulse 


encountered outside the boundary layer of a 1-ft. diameter sphere 
at a flight Mach Number of 10. The calculated cooling rate ts 
large enough to freeze the equilibrium concentration of NO at 
tained behind the shock. Since the NO reaction is exothermic, 
it is possible for boundary-layer heating to occur as a result of the 
breakdown of this nonequilibrium concentration as it enters int 
and diffuses within the boundary layer. The ‘‘freezing’’ process 
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strictly to free oxygen and nitrogen atoms since 


ul t stable at room temperatures. Information is not 


is to the recombination rates for these atoms, and 


it the relaxation time ts 





h comparable to the 
hypersonic flight problems 
ilso of 


erest 1n 


I tri xide concentration 1s importance becaus 


first component to ionize, and, at relatively low tem 
free electrons may be formed to cause the 


t shock 
At temperatures of tbout 5,000°K 


ufficient 


become highly 


GQ OO00PR 


Waves to conductive 


in thermal equi 


ri ere are approximately 10! free electrons per cu.cm 
product lost entirely by thermal ionization of NO. Caleu 
i cate that t flight Mach Numbers of 15, the con 
ductivit f the air behind a normal shock wave becomes equiva 
nt to that of sea water, 5 mhos per mete Knowledge of the 
ctr ncentration is also important because electrons may 
become efficient carriers of heat if sufficient concentrations of 
ctr re obtained 

At 1 tely the same temperature conditions at whicl 
NO f i ccurs, small quantities of NO, are formed 


Knowledge of this component may be necessary in Investigations 


radiation p} en 

In summar ve conclude that, although the dissociation reac 
tion for nitrogen may not have a significant aerodynamic effect 
I flight Mach Numbers of approximately 16, the breakdow1 


nitrogen t form other 


tent at temperatures as low as 3,000°K., and the re 


compounds occurs i significant 


ultant 


ier 





c effects 1 not be negligible even at ht Mach Num 
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On the Initial Vertical Descent of a Helicopter 
Following Power Failure 


x ~ 
Aer I gineer, Verto! Aircraft Corporation, Morton, Pa 
Y 2 1Y5¢ 
—— ESTIMATES of the rotor-speed decay and descent v« 
locity of a helicopter following power failure can be made 
which, for the first 2 or 3 sec., agree closely with more exact 


leulations 


Consider a single-rotor helicopter rhe equation of motion in 


il direction would bx 


Vert 


FORUM 


1 
i 


= rotor thrust 


= weight of helicopter 


vertical direction positive 


downward 


The equation of motion of the rotor would be 


where 


Providing the 


locit\ 


not ¢ 


pressed 


hanged, the 


rotor more 


nt 


Tinertia 


shaft torque tending to increase 


= acrodynam 


ingular \ 





pproximatel 


veight of the helicoy 
shaft torqu 

\fter power fail 

1 (6) into Eq 
loc il rot ) ‘ 
} 

i ) = 
greater the ratio of 1 
the tu f pow 
speed an l 


locity of ro 


thrust 


torque 
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Fic. 3. Predicted descent velocity as compared with reference 2 
(J = 1,000 
Eqs. (7) and (8) are compared with the calculated results of 


reference 1 in Fig. 1 and with the results of reference 2 in Figs. 


2and 3. In the case of reference 2, the information given in the 
reference was not sufficient to determine the constant Q,/w 
This was, therefore, chosen to satisfy the w/w) curve for J = 
1,000 slug ft.2, and then the rate of descent and w/w» curves for 
the other rotor moments of inertia were calculated using this 
same constant 

The simple approximations of Eqs. (7) and (8) and the more 
exact calculations of the references are seen to be in substantial 


agreement. 
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Comment on the Reliability of the Strength of 
Stiffened Cylinders in Bending and/or 
Compression 


Bertram Klein 
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I A RECENT PAPER,' some doubt was expressed concerning the 
reliability of tests (and presumably calculations) made on 
stiffened cylinders under axial or bending loading. The conten 
tion is that stiffened-cvlinder test (and possibly calculation) data 
exhibit the same large random scatter as unstiffened-cvlinder 
data. It is the purpose of this note to show that in many prac 
tical cases this belief is not justified 

For purposes of discussion, stiffened cylinders are divided 
into two broad types: 

1) Skin-Buckling 
cases, the skin is made thin and allowed to buckle well below 


Stiffened Cylinders In many practical 


the collapse load of the stiffened cylinder. The structure then 
is a gridwork of stringers and frames connected by a buckled 
medium still capable of resisting shearing action. This tvpe is 
called skin-buckling stiffened cylinders. Let us consider th 
test data given in Table | for such thin-walled skin-buckling 
cylinders.2, The specimens were of the same nominal dimensions 

Next, consider the results of calculations for thin-walled skin 
buckling evlinders loaded in pure bending*® shown in Table 2 

\ll of the cylinders discussed failed in the general instability 
mode. However, similar degrees of accuracy can be shown for 
stiffened cylinders that fail in a panel instability mod 

2) Orthotropic Cylinders When the skin gage is large enough 
or the stiffening elements close enough together) so that no ski 
buckling occurs prior to failure, the stiffened cylinder may be 
considered to fall into the orthotropic-cylinder category In 
many such practical cases, the test scatter band is still within 
reasonable limits. For example, consider the test data of 
Notice 


that the deeper the stringer compared to the sheet gage /, the less 


Table 3 for stringer-stiffened orthotropic cylinders.* 


the scatter 


TABLE | 


Maximum 


Test Com- Maximum % Deviation 
Cylinder Sheet Type pression Bending in Bending 
No Gage Test Load Moment Moment 

PIBAL 5 0.02 B 0 308 , 000 
6 0.02 B 0 359 , 000 ld 
1 0.02 B&C 13,000 296 , O00 
S 0.02 B&C 14,600 302 , 000 2 
9 0.012 B&C 13,900 227 , 000 
10 O.O12 B&C 14,200 249,000 10 
12 0.02 B&C 7.400 336 , 000 
3 0.02 B&C 7.200 362,000 7 
TABLE 2 
Maximum 
Compressive 
Strain X 10 
Cylinder theory test % Deviation 
PIBAL 1] 22.59 20.86 8.3 
PIBAL 54 18.13 16 13.3 
GALCIT i 14.0 22.9 
GALCI 


T 30 34.14 29.5 15.8 
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TABLE 3 
Stiffener Maximum a 
Sneci Spac Compressive Devi 
I ae , 
el * t ing Width Dept Load ition 
\ 
1 j 0.050 3.56 O.2d55 Qan 5.600 
1) ( 0.050 3.54 O OO) 1 460) 23 
1 t 0.052 3.56 O 186 6.650 
) 65 0.051 a.4e © 1902 5,400 14 
17 16 0.053 3.56 0.258 POU 8,370 
) 64 0.052 3.51 0.248 PO4 7. 360 Ss 
M | 5, the load is practically independent of 
vlit igt 


In conclusion, let it be stated that this brief note is not in 
tended to cover all aspects of the problem of collapse of stiffened 
cylinders loaded 


iim is to illustrate that scatter in the strength of such structures 


in bending and/or axial compression. The 


js, in many instances, entirely within the bounds of engineering 


wccuracy 
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The Asymptotic Turbulent Boundary Layer 


ClENTITIC /tricer Royal Aircraft Estab hment, Farnborough 


LAMINAR BOUNDARY LAYER on a flat plate with uniform 


suction develops asymptotically toward the 


form 
, 
ufU=1-— 
ere the boundary-layer thickness and the velocity profile are 
independent of 4 sv analogy, it has been suggested® * that 
the turbulent boundary laver with uniform suction should also 
od 


briefly describes experiments in which the 
laver with uniform 


equilibrium condition, with both thickness and velocity 


1 


levelop to a common asymptotic form. The present note 
turbulent boundary 
suction was, in fact, found to achieve an 
profile 
ndependent of However, the experiments suggest that the 

isured velocity profile corresponding to this asymptotic con 
lition is not universal, and, in other respects also, the behavior 
f the turbulent boundary laver with uniform suction evidently 
liffers from that of the corresponding laminar layer 


Experiments were carried out on a flat perforated metal sheet 


ft. long which replaced the floor of the working section of a 
ind tunnel This perforated surface was preceded by a 4-in 
lid entry length which formed the upper lip of a bleed through 
which the tunnel boundary layer was removed \ transition 


wire Was fitted on the front of this solid surface so that the bounc 


sequent devel 


iry laver at the start of suction was fully turbulent. The sub 


pment of the boundary laver was measured for 


veral different suction velocities at a tunnel speed of 60 ft. per 
For the first series of measurements the perforated surface 
s left exposed; for the second it was covered with calendered 


fabric to simulate a continuously porous surface 


* Formerly, Research Student, Engineering Department, University of 


imbridge, England 
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turbulent boundar 


Fic. 1. The development of an initially 
l 1 perforated surface 


iver for various rates of suction through 


The development of the boundary layer on the perforated sur 
face at several different rates of suction is shown in Fig. 1 It 
will be seen that, for a particular rate of suction | l 0.0073 
thickness substantially constant 


the boundary-laver remains 


over Virtually the full length of the plate; for greater suction vi 

locities the turbulent boundary laver decays toward the laminar 
asymptotic form, while for smaller suction velocities it continues 
nvlon-covered surface 


to grow The results obtained with the 


were qualitatively similar, but the suction velocity required t 


boundary-laver thickness was 


With both surfaces, the 


maintain a constant approx! 


mately velocity pr 
files of the 


unchanged following 


4) per cent less 


constant-thickness layer remained substantially 
in initial period of adjustment from tl 

solid-boundary profile at entry. Moreover, in spite of the ver) 
different suction velocities required in the two cases, the velocity 


profiles of the constant-thickness layers were identical except 


very close to the surface (see Fig. 2 
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Fic. 3. Rate of loss of mean-flow energy in turbulent boundary 


laver with and without distributed suction 


With increased lengths of solid entry and conditions otherwise 


unchanged, it was possible to maintain a turbulent boundary 


layer which was substantially constant in momentum thickness, 
but the 


velocity profiles now changed continuously 


plate, tending away 


along the 
from rather than toward the equilibrium 
profiles measured earlier 

It is of interest to note that, when conditions are independent 
of x, the distribution of turbulent shear stress through the bound 


ary layer can be simply obtained from measurements of the 


mean velocity profile; the boundary-layer equation reduces to 


Ou/Oy) = (1/p) (Or/Ov 
which upon integration gives 
Fre on : ‘ a . ' 
pu'v’ = pv,U;[1 — (u/l — (v/2, | Ou/oy); 
since r = plOu/Oy) — pu'v’, andv = —% const. fro con 
tinuity 


Distributions through the boundary layer of 
production and direct dissipation terms have been obtained from 


the mean-flow 


the turbulence 


velocity profiles and are plotted for the nylon 
surface in Fig. 3, which also shows results obtained by Klebanoff! 
for the flat plate. It 


markedly different. 


solid will be seen that the results are 


Plotting the asymptotic profile in the form u/ Ul, Uyy/v 
gives a straight-line relation which extends to the edge of the 
boundary layer, thus suggesting that intermittency is absent 
It is believed that the asymptotic turbulent boundary layer 
exhibits sufficient 
The 
analogy 


features of interest to repay further study 
term asymptotic has been used partly because of the 


with the laminar case and partly to differentiate this 
equilibrium layer from the type studied by Clauser® and others, 
where only the velocity defect of the outer part of the layer ex 
hibits longitudinal similarity 
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’ | NHE HYPERSONIC BOUNDARY LAYER on a flat plate is associat 
with a region of inviscid rotational flow behind cury 
shock wave from the plate leading edge lo provide an unde 


standing of the effects of free-stream vorticity o 


£ a viscous boundary layer, the writer 
example of a viscous boundary layer in a free str 
vorticity The reformulation of 
problem wherein the pressure gradient induced by 
ment effect of the boundary laver is taken 


highlights of the analysis can be most simply pres 


1 ; 
aiscusseda i sin 


im of const 


present note presents 


the displa 


into account | 


‘ase of steady incompressible two-dimensional flow I 


Navier-Stokes equations can be written in the usual notation 


follows 


u(Ow/Ox) + v(Ow/Oy) = v[(O%w/Ox?) + (d%w/d 
w = (Ou/OY Ov/ OX 

The continuity equation is 
Ol / OX r (Qv/Oy) = U 


By examining the order of magnitude of the var 
these equations, in the 
famous boundary-layer approximation, it can be shown that in 
fluid of small viscosity the variables w, u, 7 


must vary in accor 


ince with the following boundary-layer equations: 


( 


10us terms 1n 
same manner that Prandtl originated his 


1 


u(Qw/Ox) + v(Qw/dy) = v(0%w/dy? la 
w = Ou/OyY 2a 

Ou /Ox dv/dy) = 0 3 

ina thin layerO0 < y < 6,6/x < 1 This thin boundary layer 
adjacent to the surface is a region of strong vorticity In the 


siinple problem of a uniform shear flow past a flat plate, the free 
stream condition is 


In reference 1, it was assumed that 7 = ()—1.e., the 


When the 
laver displacement effect is taken into account, this condition 1s 


boundary 


layer displacement effect was neglected boundar 


expressed more properly as 


th) t 
rhe usual no-slip conditions on the surface are 
u“uwz=v=—U0 5 
Eq. (la) can be integrated once with respect to y. Thus 
u( OU, OX + v(Ou/Oy) = v(O°u/Oy* + f(x 
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n the behavior 
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onditions of Eqs } d (4b Eq. (6 s > I r. s e She / ease ‘ ‘ 
~ J P| Journal of the Aeronautical Scienc Reader Forun 
3 22, No. 10, pp. 724, 725, October 19 
w 7 Schlichting, H Bou l 1 Or MeGra H Book 
Company, Inc., 1955 
} ts the duce¢ ressure g i due to the dis 
! ce ffect of the boundary laye ntroduct mw the 
C m f tion ¥ such that s 
‘ Ov °) Oy O S) 
( lation Of Mot1o d its boundary conditions are 
A Shock-F ree Airfoil 
»y /O vy /Ox0) Oy O O*y 04 
0"'y 9 Jy 9 3 H. C. Levey 
Ov Ox = OY/O 10 Aeronautical Research Laboratori epartme 
Australia 
—> Ory O = w 1] y 30, 195¢ 
ence I, let 
INTRODUCTION 
v\ V3 Blasius variable 12 
TT EQUATIONS which govern the steady tw imensional 
é r= «& V ex i vorticity number 13 flow of an inviscid ideal gas admit solutions whi have the 
d express y in terms of an ascending power series in & with property of transonic isentropic recompression In particular 
functions of n as coefficient it is possible to find profiles which theoretically give rise to shock 
free supersonic regions which placed in an appropriate subsoni 
sociated ¥ V vu re Tt 14 stream rhe fact that this phenomenon does not seem to be 
curved Formal substitution of Eq. (14) into Eqs. (9), (10), and (11 observed experimentally is a matter of some controvers\ 

n unde ields the following sequence of equations though evidence is accumulating to show that the reason ma\ 
behavior r be “the non-existence of neighbouring flows Before any ex 
simpl i ’ woes 1S perimental investigations can be made, it is necessary to hav: 
constant ) = ) 16 profile for which the exact solution exhibits a comparativel 

| of the _ , 17 large shock-free supersonic region, and such profile is give! 
lisplac this note 
it Phe - ro = m 18 The hodograph method is used, in which the profile shape can 
d im the ) = — 19) only be specified in advance in fairly general terms It was 
" rh : judged that typical conditions would be met by yvymmetrical 
ition #s =] ad profile with a thickness ratio of about 10 per cent in a free stream 
he unction defined by Eqs. (15 17) is the Blasius function of Mach Number 0.8. Since the method of « insion in eigen 
The factor (f fo’ )n— aw =—1.73 Thus, the equation for f functions and analytical continuation throughout the hodograph 
n be solved with the boundary conditions as specified plane®: * leads to slowly convergent series, it was decided to ust 
2 It can be stated that, to the order of terms, the induced pres more complex solutions developed by Cherry vecause of the 
sure effects as well as the free-stream vorticity effects are ac rapid convergence of the series involved hese solutions lead 
ounted for in the present formulation \s it is demonstrated to a profile pointed at both ends (despite a number of attempts at 
bove, to the zeroth approximation, neither the free-stream vor modification ), but this was not thought to be a serious limitation 
ermis in ticity nor the boundary-layer displacement effect contributes since the solutions are only valid for zero incidence 
ited his iny change on the viscous boundary layer. But the free-stream It was found that two solutions of the type given by Cherry 
hat ina vorticity effects may become significant as soon as the bound were necessary, the first (with subscript 4 ) to give the basic shape 
accord- iry-layer displacement effects become important. In_ low and the second (with subscript 2) to modify the thickness ratio 
peed boundary-layer calculations these effects are both small In addition, a particular elementary solution was added to avoid 
In some hypersonic boundary-layer investigations, however, the occurrence of a limit line in the flow field 
both effects are expected to be significant As a matter of 
a fact, the induced pressure effects on a hypersonic boundary FORMULAS 
2 layer cannot be accurately studied without a full knowledge of At the point x, y, let g, 6 be the polar velocity components 
Laver the free ee “Varres Conners Thus, the study of the with g normalized so that g = 1 is the vacuum speed, and denot 
tn the inviscid flow held around a body at hypersonic speeds becomes a the sonic speed by q.[= 1 1 The positive con 
sata necessary condition for successfully carrying out hypersonic stants a, 8, \ are now defined by 
viscous-flow calculations. In general, it can be argued that one sa ; 
1 the effects of free-stream vorticity is to introduce a modified Zod} 4 7: : a l 
ta pressure field. Consequently, the skin friction, heat transfer, iat 
ndary boundary-layer separation, and laminar stability characteristics ind a transformation is made to the variables g, ¢ through 
seine re all affected These effects are important as soon as the dis ee eae Ss aevtan sin o/(1 ” 2 
tion is placement effects become significant ‘ ; 
where « is defined in terms of the free-stream speed b 
$b) t seapanaiees «x = arctanh } [1 dao?/gs?)|'/? (1 72) 1/4} 
Small \ ty, Journal of the Aeronautical Science Reader Forum arctanh ; “as ws 
E <2, No. 9, pp. 651, 652, September, 1955 he solutions used are then given (with real parts understood 
by 
as t required that the inviscid free stream have constant vorticity V4 = s(] 2 7X (> 
f net considered here, this condition is, strictly speaking X44 = (g0Q24/0 1 [%ag sin ¢ Q | Do « P } 
ory /o ary /d ’ Y,=2 
hrough t 0*-y/0 y2y 0 0 R < 
the present calculations, Ee 1] consistent with the boundary *Acknowledgment is made to the Chief Scientist, Department of Sup] 
blish this irticle 


er appre mation Australia, for permission to pu 





1150 FOURNAL OF THE AERONAUTIC 


M-0-755 
M0536 


M-t3 Mel 


M-0-667 





°o 





0 0 20 


» 


= 
50 60 70 80 90 100 
% 


Fic. 1. Airfoil with shock-free supersonic region at 1/.. = 0.8; 
upper surface and sonic line 


= are? be bag?” - 1+ -—s , & r-2& 


V sin 6, y= X siné@ + JV ecos@ (6 


tL SCIEZNCES DECEMBER, 956 


Of), 3) = * l ] 
.? > 
re sa" 
l l l 
( : ) f , . , 
2 } 3 | 
The coefficients ¢ —) . are tabulated by Cherry ind 
Gnr, On, Cnp AYE privately tabulated by the same author I 
tary solutions are given by 
vw" = (n 1 )g sin néF 
X”) = n(1 g?)-P « cos 29 } F,( 4? 
24 F D 1)} 
) = (1 ] ( sin 29 {nF,(g 
29? F , 1)! 
where the F/,(qg?) are the standard hypergeometric functions 
occurring in hodograph theory 
DISCUSSION 
rhe final solution used was of the form 
YW = 4 + AW, — 2B,Y 10 
where A» and &, were adjusted so that VW = 0 gave the desire 
profile. The occurrence of a limit line was avoided by adjusting 


the profile streamline in the hodograph plane to avoid tangenc 
with the characteristics, and it was found that, within the limita 
tions of the initial requirements, the greatest supersonic Mac 
Number—circa 1.13—was attained with A, = 1.4 and 148A;; = 
107 with the other B, zero. The field points have been con 
puted to five significant figures, and the upper surface of tl 
profile and the sonic line are shown in Fig. 1 There is a flatten 
ing near the downstream sonic point, which is not discernable or 


the figure, but otherwise the profile exhibits no marked peculiari 


ties 
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